Hi Ping,
You can use zz$coefficients[,4] to get the p values for each estimated
coefficients in your context.
Wuming
On 9/11/05, Ping Yao [EMAIL PROTECTED] wrote:
Wuming:
Thanks for your help.
I use the fuction:
call(fstatistic,zz)
call(p-value,zz)
I can get each variable P-values,but I can't get P-value of the model.
How can I do ?
one of the results is following :
Call:
lm(formula = gen.fat ~ snp_marker)
Residuals:
Min 1Q Median 3Q Max
-10.5455 -3.0481 0.4545 3.9519 6.9519
Coefficients:
Estimate Std. Error t value Pr(|t|)
(Intercept)13.0481 0.4518 28.881 2e-16 ***
snp_markerallele2 0.5107 0.9102 0.561 0.5753
snp_markerBoth 1.4974 0.6927 2.162 0.0318 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 4.607 on 212 degrees of freedom
Multiple R-Squared: 0.02166,Adjusted R-squared: 0.01244
F-statistic: 2.347 on 2 and 212 DF, p-value: 0.0981
I use the code :
zz-summary(lm.D9 - lm(gen.fat~snp_marker))
coe-coef(lm.D9)# the bare coefficients
if (coe[2]=.05||coe[3]=.05||coe[4]=.05||coe[5]=.05) {
cat(phenotype is = ,x[j] , \n)
cat(snp marker is = ,x[i] , \n)
sign-call(fstatistic,zz)
call(p-value,zz)
#print(coe)
print(zz)
}
On 9/10/05, Wuming Gong [EMAIL PROTECTED] wrote:
?summary.lm and check the Value section.
Wuming
On 9/10/05, Ping Yao [EMAIL PROTECTED] wrote:
Hi:
I use lm (linear model) to analyze 47 variables , 8 responses
So I use loop to finish it .
I want the program to show the results that P-value is less than 0.05.
How can I cite the P-valus from lm result ?
Ping
The code:
#using LM to model general fati
for (j in 48:52) {
for (i in 3:46){
gen.fat-y_x[,j]
gen.fat-as.numeric(gen.fat)
snp_marker-y_x[,i]
x-colnames(y_x)
#snp_marker-as.matrix(snp_marker)
#mode(snp_marker)
cat(phenotype is = ,x[j] , \n)
cat(snp marker is = ,x[i] , \n)
zz-summary( lm.D9 - lm(gen.fat~snp_marker))
print(zz)
return
}
}
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html