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)