求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

为什么我的R程序运行的时候没有显示错误，但就是没有输出结果？求教高手帮忙解决

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：输出结果 R程序 求助

Could you provide the code?

2# junuylia offset=2
setwd("/home/ubu/zhangbb/")
x=read.csv("ptpn22-500.csv",header=T)
y=x[,1]
snp1=x[,-1]
snp=snp1[,-17]
names=colnames(snp)
n=ncol(snp)
p=vector()
q=vector()
for(i in 1:(n-offset+1))
{

id=i:(i+offset-1)

data=snp[,id]

sename=names[1:(n-offset+1)]

pr=prcomp(data,scale=T)

pre=pr$x

c1=pre[,1]

c2=pre[,2]

log=lm(y~c1+c2)

p=log$stats[5]

}
result=cbind(sename,p,q)
col.names(result)=c("rs","pca")
#col.names=c("rs","pca")
path=paste("pvalue_",offset,".csv",sep="")
write.csv(result,path,col.names=T,row.names=F,quote=F)

offset没有定义

what do you expect to see?
by using log=lm(y~c1+c2), the linear summary are stored in log. You may use summary(log) to check. If you want to print the output out, use print function.