s******1
发帖数: 39 | 1 请教高手，下面这一段R code 是什么意思？谢谢
getLMOutFunc = function(lmOut){
highSEL = coef(lmOut)[,"Estimate"][["(Intercept)"]] + coef(lmOut)[,"
Estimate"][["pred"]]
lowSEL = coef(lmOut)[,"Estimate"][["(Intercept)"]] - coef(lmOut)[,"
Estimate"][["pred"]]
sampleN = as.numeric(lmOut[["ngrps"]])
if(coef(lmOut)[row.names(coef(lmOut))=="pred",][[4]]<.05) {
sig = coef(lmOut)[row.names(coef(lmOut))=="pred",][[4]]
} else {
sig = "ns"
}
return(list("High SEL" = highSEL, ... 阅读全帖 |
|
S********a
发帖数: 359 | 2 有个conditional logistic regression,想看AIC,但是结果里并没有给出来，请教什么
code或者命令把AIC求出来。谢谢。
> result0<-clogit(case~hr0+strata(subjid), data=mydata)
> summary(result0)
Call:
coxph(formula = Surv(rep(1, 9704L), case) ~ hr0 + strata(subjid),
data = mydata, method = "exact")
n=8397 (1307 observations deleted due to missingness)
coef exp(coef) se(coef) z Pr(>|z|)
hr0 0 0.001158 1.001158 0.005376 0.215 0.83
exp(coef) exp(-coef) lower .95 upper .95
hr0 1.001 ... 阅读全帖 |
|
h********g
发帖数: 5 | 3 Finite Element Analysis of the Indentation-Induced
Delamination of Bi-Layer Structures
Ming Liu and Fuqian Yang∗
Materials Program, Department of Chemical and Materials Engineering,
University of Kentucky, Lexington, KY, 40506, USA
Journal of
Computational and Theoretical Nanoscience
Vol. 9, 851–858, 2012
Contact deformation can cause local damage of mechanical structures and lead
to structural failure
of mechanical devices. In this work, we use the finite element method
to analyze t... 阅读全帖 |
|
D********r
发帖数: 135 | 4 看看图吧，空难死亡从1984到2017一直在线性下降(R2=0.65)，跟1994年枪禁无关。看
在你还想学点统计的份上，从数据到code都给你了，好好研究一下。
Year Deaths incidents
2017 399 101
2016 629 102
2015 898 123
2014 1328 122
2013 459 138
2012 800 156
2011 828 154
2010 1130 162
2009 1108 163
2008 952 189
2007 981 169
2006 1298 192
2005 1463 194
2004 767 178
2003 1233 201
2002 1418 197
2001 1539 210
2000 1586 198
1999 1150 221
1998 1721 225... 阅读全帖 |
|
B********e
发帖数: 10014 | 5 f(x)=(x^{m+1}-x)/(x^{n+1}-1);
f'(x)=f1(x)/(x^{n+1}-1)^2,
where f1(x)=-(n-m)x^{n+m+1}-(m+1)x^m+nx^{n+1}+1.
want: f1(x)\leq 0 when x\geq 1.
sufficient: f1(1)=0, f1'(x)\leq 0, when x \geq 1.
f1'=x^{m-1}*f2(x), with f2=something.
then sufficient: f2(1)=0, f2'\leq 0.
which is not difficult because
f2'=-(n+m+1)(n-m)(n+1)x^n+n(n+1)(n-m+1)x^{n-m}.
homework: complete it by showing the absolute value of the coef of first
term
is larger than the coef of the second term, plus the fact n>n-m and x>1. |
|
g********r
发帖数: 8017 | 6 abline(fm$coef)
另外检查一下fm$coef是不是对的。abline这种函数不应该出问题。我怀疑是
coefficient不对，所以线画到plot范围外边去了。 |
|
q**j
发帖数: 10612 | 7 打听一下，glmnet default给前面加一行constant。但是我的X matrix里面已经有了这
个常数了，能否改变glmnet？还是我自己要改变了？
还有
set.seed(1010)
n=1000;p=100
nzc=trunc(p/10)
x=matrix(rnorm(n*p),n,p)
beta=rnorm(nzc)
fx= x[,seq(nzc)] %*% beta
eps=rnorm(n)*5
y=drop(fx+eps)
px=exp(fx)
px=px/(1+px)
ly=rbinom(n=length(px),prob=px,size=1)
set.seed(1011)
cvob1=cv.glmnet(x,y)
以后
coef(cvob1)的结果和
coef(cvob1,s=cvob1$lambda.min)不一样。
我想第二个应该是我们最后想要的东西，但是第一个是什么？
最后，cv.glmnet里面如何能够取出来那个balance L1和L2 penalty的常数呢？
... ...你好歹应该先知道Lasso和Ridge是什么吧？ |
|
a*******g
发帖数: 80 | 8 Hi, I used Dr. Harrell's rms package to make a nomogram.
Below is my code for nomogram and calculate total points and probability in
original data set used for building the nomogram. My question is how I get
the formula for calculating the survival probability for this nomogram. Then
I can use this formula to do validation by using other data set.
f1 <- cph(Surv(retime,dfs) ~ age+her2+t_stage+n_stage+er+cytcyt+Cyt_PCDK2 ,
data=data11,
surv=T, x=T, y=T, time.inc=5)
surv<- Survival(f1)
surv10 <- ... 阅读全帖 |
|
k***g
发帖数: 7244 | 9 其实政治学里作中国基层研究的挺多的，不过好像在国内没有什么影响力；
也许是因为写得太学术了使得一般的读者失去兴趣了，主要是一方面要试图
放进西方政治、社会学的理论框架，引入大量的一般人不懂的术语和专有名
词，另一方面大量的使用计量的方法，没有相关背景的人很难理解，比如给
一个 estimation 结果和一个 Hessian matrix，一般的读者不会明白如何用
Hessian inversion 的 diagonal 来解读 estimation 的 coefs;
而这个又是由于期刊的审稿制度导致的，如果你只是讲一个story，没有 th
eory 和 相应的 方法，很难发表在好的期刊上；
自然科学其实也有这样的问题。
《美国政治学评论》近十年来唯一篇关于中国国内政治的文章（印象如此）
，就是关于中国农村的 public goods 的问题：如果非民选的官员的 accou
ntability 很弱的话（因为没有选举束缚）? 那么他们如何给人民提供 pub
lic goods，他们为什会有这样的 incentive ? 社会的稳定是如何实现的？
这篇文章其实回答的是《中县干部》里问题... 阅读全帖 |
|
k***g
发帖数: 7244 | 10 其实政治学里作中国基层研究的挺多的，不过好像在国内没有什么影响力；
也许是因为写得太学术了使得一般的读者失去兴趣了，主要是一方面要试图
放进西方政治、社会学的理论框架，引入大量的一般人不懂的术语和专有名
词，另一方面大量的使用计量的方法，没有相关背景的人很难理解，比如给
一个 estimation 结果和一个 Hessian matrix，一般的读者不会明白如何用
Hessian inversion 的 diagonal 来解读 estimation 的 coefs;
而这个又是由于期刊的审稿制度导致的，如果你只是讲一个story，没有 th
eory 和 相应的 方法，很难发表在好的期刊上；
自然科学其实也有这样的问题。
《美国政治学评论》近十年来唯一篇关于中国国内政治的文章（印象如此）
，就是关于中国农村的 public goods 的问题：如果非民选的官员的 accou
ntability 很弱的话（因为没有选举束缚）? 那么他们如何给人民提供 pub
lic goods，他们为什会有这样的 incentive ? 社会的稳定是如何实现的？
这篇文章其实回答的是《中县干部》里问题... 阅读全帖 |
|
s********z
发帖数: 5411 | 11 那天跟一个伊朗朋友吃饭，他问我以后怎么办。
我说工作几年，然后会中国。 他问为啥，我说想给自己的祖国做点事情。 然后他滔滔
不绝想告诉我，不一定回国才能做事情，在国外也可以，也许更有效。
我就在想，大家在这里能不能想出来我们力所能及能做的事情呢？我先说说我能想到的
，希望能抛砖引玉。
1 买Made in China. 不管是自己用，还是给公司买东西，能买国货就买国货。比如我
们公司都买lenovo。
2 给charity捐钱。像地震这样的事情不用说，平时也可以捐。 比如，very least,在
amazon上买东西可以通过coef的链接。
3 把自己的钱投到国内股市。 我想这样不但支持国内企业发展，而且对自己的
portfolio也起到diversification的作用，况且中国毕竟是发展中国家，上升的潜力大。
4 帮中国做PR. 记得我当时去反对藏独，支持奥运的游行的时候就花了不少时间教育
美国人。 平时，对外国朋友，或者公司同事什么的可以教育教育。
还能做啥呢？ |
|
B*Z
发帖数: 7062 | 12 所以单拿质量一个来说事就是瞎说了。在一些特殊情况下，loaded truck还可以比空车
停得快。（if friction coef increases with weight, but still below some thres
hold.)
看看texas考试题。注意，我没说满载的总是停得快。
http://www.quest4.org/txcdl/tx-cdl.htm
29. Empty trucks have the best braking. True or False?
False. Empty trucks require greater stopping distances, because an empty
vehicle has less traction. It can bounce and lock up its wheels, giving muc
h poorer braking. (page 2-15) |
|
s****n
发帖数: 1237 | 13 the optimal solution may not have max coef for Dn, see my counter example
above.
, |
|
|
a***e
发帖数: 2094 | 15 Rayrose一般情况下没错， 但是这回没有考虑细致， 比如30%精神出轨算不算，或怎么
算呢？ 40%， 50%呢? 最好做个sampling,看看distribution. 如果是normal
distribution的话， 就容易判断。 先要测试是不是normal, 可以这样来test:
JB = n/6 x [S^2 + (K-3)^2/4]
30%精神出轨一个月，100%精神出轨但只有两天，哪个更严重？比一下duration.
Macaulay duration = modified duration x (1 + yield/2）
如果只有出轨(down)和不出轨(up)两种状态，这样来看probability：
P(up）= (e^rt - D)/(U - D）
如果有出轨的风险，最好用F来hedge一下。
optimal hedge ratio = corr. coef.(S,F）x [st.d(S)/st.d(F)]
我看书看昏头了，嗬嗬 ^_^ 还是让pro statistician 来发话。 |
|
w********c
发帖数: 2632 | 16 ☆─────────────────────────────────────☆
jetchen (飞机) 于 (Fri Nov 2 14:44:48 2007) 提到:
比如已知有n个 Gaussian mixture, 其中已经有m个gaussian的mean和variance知道了
， 要求剩下的(n-m)的mean, variance， 和mixing coef.
是不是在EM update的时候加一些限制条件就行了？
☆─────────────────────────────────────☆
PaulPierce (Paul) 于 (Fri Nov 2 18:11:17 2007) 提到:
E的时候，不重新计算知道的m个m & var
☆─────────────────────────────────────☆
sonyisme (偶静感乖类 :)) 于 (Fri Nov 2 18:11:56 2007) 提到:
re,hehe
☆─────────────────────────────────────☆
PaulPierce |
|
b*****d
发帖数: 7166 | 17 下面这个函数，输入后说else 是invalid syntax。
def coef(m):
p=1
if (m%2)==0:
for i in range(1,m,2):
p=p*i
return p*sqrt(pi)/2^(m/2+1)
else:
return factorial((m-1)/2)/2
为啥python不能用if。。。else？ |
|
y*******g
发帖数: 6599 | 18 我试了没这个问题
只是
Traceback (most recent call last):
File "", line 1, in
File "", line 6, in coef
TypeError: unsupported operand type(s) for ^: 'float' and 'int' |
|
|
a**n
发帖数: 3801 | 20 the coef estimates cannot be taken seriously due to
big st dev.
coefficients |
|
c****o
发帖数: 49 | 21 what happens is ur filter is operating on high sample rate, so coefs will
not converge... filter taps = (80/22)*fs/f_transition. conclusion is, boy u
need to do sth smart |
|
H*******s
发帖数: 382 | 22 【 以下文字转载自 Mathematics 讨论区 】
【 原文由 HappyDays 所发表 】
俺做了一个2^3 factorial design,统计每年度香港/大陆最流行歌曲的时间长度。
三个factor分别是Area(HK, Mainland), Award(Gold, Silver),Year(1980,1990),
用Minitab算的结果是：
Term Effect Coef
Constant 3.9263
Year 0.2525 0.1262
Area -0.3125 -0.1562
Award -0.0675 -0.0338
Year*Area 0.6075 0.3038
Year*Award -0.3675 -0.1838
Area*Award 0.0675 0.0337
Year*Area*Award 0.1 |
|
h********g
发帖数: 5 | 23 Ming Liu and Fuqian Yang, Finite element simulation of the effect of
electric boundary conditions on the spherical indentation of transversely
isotropic piezoelectric films, Smart Materials and Structures 21 (2012)
105020 (10pp)
Finite element simulation was used to analyze the effect of electric
boundary conditions on the
indentation deformation of a transversely isotropic piezoelectric film
with the contact radius
much larger than the film thickness. Six different combinations of... 阅读全帖 |
|
h********g
发帖数: 5 | 24 Finite element analysis of the spherical indentation of transversely
isotropic piezoelectric materials
Modelling Simul. Mater. Sci. Eng. 20 (2012) 045019 (15pp) doi:10.1088/0965-
0393/20/4/045019
Abstract
Finite element analysis was used to analyze the indentation deformation of
a transversely isotropic piezoelectric material (PZT-4) by a rigid spherical
indenter. Three cases were considered in the analysis, which included (a)
the indentation by an insulating indenter, (b) the indentation by a
c... 阅读全帖 |
|
o****i
发帖数: 19 | 25 我觉得和阶段是相关的吧。
博士生，博士后，还有有位置的人不能一概而论。
我 postdoc
读一些关于QCD 3-loop 的 splitting function and wilson coef. 计算文章感觉还行
，虽然得下点功夫。但是最后还都放到程序里了。 |
|
V****n
发帖数: 651 | 26 I think you made a couple of mistakes. The unit of thermo conductivity is W/
m-K, your equation fails to balance the unit.
The transient heat conduction equation is
k*(d2T/dx2+d2T/dy2+d2T/dz2)+Q' = rho*c*dT/dt (figive my non-rigorous
notations here)
where k is the heat conduction coef., Q' is the internal heat generated per
unit volume per unit time,
rho is density and c is specific heat (for solid Cv and Cp are about the
same)
When a person dies Q'=0. You've already calculated the right hand si |
|
|
s******d
发帖数: 2730 | 28 assume you have x, y, state three variables
x, y are the data with the five states pooled together
and state is the label (say, 1, 2, 3, 4, 5, etc)
You only need to use
> plot(x, y, col=state)
> abline(lsfit(x,y)$coef) |
|
d*******1
发帖数: 854 | 29 好吧，我把Ripley书里的例子和我自己的数据side by side run 了一下。他的好好的
我的就是死活不行（abline(fm)和abline(fm$coef)应该都行的）
而且我的散点数据看上去很好的， 难道一个简单的regression 还会fit 错了？ |
|
l*********s
发帖数: 5409 | 30 How can coef change at all? OLS estimate shall be unique.
but
X1,
of |
|
p***l
发帖数: 1775 | 31 multicolinearity解决不了是不是只好转向ridge regresion了?
coef出来之后不知道要做哪些diagonastic |
|
S********a
发帖数: 359 | 32 > unadj.lrm
结果如下
Coef S.E. Wald Z P
Intercept -3.60925 0.132424 -27.26 0.0000
pmavg 0.02205 0.009483 2.33 0.0201
我想取pmavg 的standard error （i.e.,0.009483）
怎么样写R code呢？
> unadj.lrm$SE 好像不对呀 |
|
hs
发帖数: 1549 | 33 summary(xxxx)$coef[2,2]
xxxx是 unadj.lrm |
|
S********a
发帖数: 359 | 34 完整的结果是这样的， 如何取SE呢？谢谢大家
> unadj.lrm
Logistic Regression Model
lrm(formula = combpreecl ~ pmavg, data = mydata2)
Frequencies of Responses
0 1
118691 4360
Obs Max Deriv Model L.R. d.f. P C 123051
7e-10 5.41 1 0.02 0.504
Dxy Gamma Tau-a R2 Brier
0.008 0.015 0.001 0 0.034
Coef S.E. Wald Z P
Intercept -3.60925 0.132424 -27.26 0.0000
pmavg 0.02205 0.009483 ... 阅读全帖 |
|
S********a
发帖数: 359 | 35 > sqrt(diag(adj.lrm$var))[12]
pmavg
0.009879157
> sqrt(diag(adj.glm$var))[12]
[1] NA
> summary(adj.glm)$coef[12,2]
[1] 0.009879123
对lrm用你的这个能得出结果，对glm用summary能得出结果。
谢谢！ |
|
u*****3
发帖数: 796 | 36 不想看这些info直接fit回车，或者coef(fit) |
|
i**********d
发帖数: 26 | 37 简单讲一下，原理，怎么设定lamba
疑问是，每次coef小的var的beta被shrink到0以后，这些var就不在model里面了吗？还
是在模型里面，但是假定beta = lamba |
|
c***z
发帖数: 6348 | 38 ## build data frame
work <- c(12, 14, 4, 16, 12, 20, 25, 8, 24, 28, 4, 15)
edu <- c(6,3,8,8,4,4,1,3,12,9,11,4)
income <- c(34.7, 17.9, 22.7, 63.1, 33.0, 41.4, 20.7, 14.6, 97.3, 72.1, 49.1
, 52.0)
studay.df <- data.frame(cbind(work, edu, income))
## linear model
model_3 <- lm(income ~ ., data = studay.df) # OLS
summary_table <- data.frame(summary(model_3)$coefficients)
colnames(summary_table) <- c("coef", "std.error", "t_value", "p_value")
summary_table$regressor <- row.names(summary_table)
s... 阅读全帖 |
|
P*****r
发帖数: 554 | 39 mulitple linear regression 解读coef的时候,一般会说其他量不变的情况下
你这还做过transformation,应该是transform后1 unit 增加会增加.192 |
|
w******e
发帖数: 142 | 40 还有一种方法是不要选最低cross-validation MSE的lambda来变量选择，而用1个SD以
外的变量更少的lambda来选变量再重新OLS. glmnet package默认的coef(model)出来
的就是1 sd away的系数。忘了在讲义还是书上看见过原作者这样提出来过。 |
|
c********h
发帖数: 330 | 41 没用过sas，不过你先试下linear model里fixed effect的R2是不是也很小
coef significance不一定等价于R2很大 |
|
w****u
发帖数: 69 | 42 需要写个r code对有2个strata的data，画生存曲线图，
在网上找到code如下，程序的结果是对的，
老板问了个问题答不上来，投文章今天是deadline了，求救啊！！！！！！
其中ad5是我们的strata，treatment （1，0）表示实验组/非实验组，
贴了和第一个strata： Y 有关的部分，
先算每个每个strata的人数，对每个strata算baseline的hazard h0，然后S=exp(-H0)
×exp（beta* X1),
问题是最后算x1的时候，为什么要减去那个c(mean(treatment)呢？？？？
D=read.table('D:/survial function.txt',header=T)
attach(D)
S=Surv(time,infection)
m1=coxph(S~strata(ad5)+treatment)
d.phm=coxph.detail(m1)
#strata are Y and Z; number individuals in Y, Z:
nY=d.phm$strata[1]
#h0 for Y
timesY... 阅读全帖 |
|
R*****0
发帖数: 146 | 43 没有上下文，不过感觉fit_res$coeffs被默认设成了数组，下标[n]只能用于数组（
矩阵，向量）的运算和读写。试试这个：
fit_res <- list(coeffs = vector("list",100))
for (n in 1:100){
this_fit <- lm(...)
fit_res$coeffs[[n]] <- this_fit$coef
}
fit_res$coeffs[[50]][1:6]
以后直接调用双括号下标[[n]] |
|
K*****2
发帖数: 9308 | 44 如果model没有变，x1还是4个level，应该用0.3
如果用0，应该把x1的c和d两个level合并，重新计算全部的coef |
|
s******1
发帖数: 39 | 45 请问高手， 这样计算cohen's D， effect size对吗？
modOut = lm(GPA ~ IEP + FRL + ELL + pred, data=reg2data)
out = summary(modOut)
cohensd = coef(out)[,"Estimate"][["pred"]]/ sd(predict(modOut))
看起来cohen's D= "'pred'的regression coefficient" 除以 "回归模型的预测值"的
standard deviation.
请问这样计算cohen's D 对吗？ |
|
b*****d
发帖数: 7166 | 46 【 以下文字转载自 Programming 讨论区 】
发信人: biokold (kold), 信区: Programming
标 题: python 问题
发信站: BBS 未名空间站 (Sun May 6 22:28:13 2012, 美东)
下面这个函数，输入后说else 是invalid syntax。
def coef(m):
p=1
if (m%2)==0:
for i in range(1,m,2):
p=p*i
return p*sqrt(pi)/2^(m/2+1)
else:
return factorial((m-1)/2)/2
为啥python不能用if。。。else？ |
|
s******d
发帖数: 2730 | 47 你只能比较两个model的fit。
test coefficient不make sense因为你x换个scale，coef就变了。
你只能看prediction performance.
看r^2阿，anova F test. |
|
