Open RandiLGarcia opened 7 years ago
upp=parti+api(xmean+xsd)
lwer=parti+api(xmean-xsd)
covpint = yay$varBeta[kk,jj] varp = yay$varBeta[kk,kk] varint = yay$varBeta[jj,jj]
seupp = sqrt(varp + 2(xsd+xmean)covpint+(xsd+xmean)^2varint) selow = sqrt(varp + 2(-xsd+xmean)covpint + (-xsd+xmean)^2varint)
plow = 2-2pnorm(abs((parti+api(xmean-xsd))/selow))
See Dave's email on Jun 21. Can be used instead of re-centering.