Skip to content

Commit

Permalink
DeBUG taxa index and GTindex
Browse files Browse the repository at this point in the history
  • Loading branch information
jiabowang committed Apr 24, 2024
1 parent f4b20a7 commit 6fd7265
Showing 1 changed file with 5 additions and 5 deletions.
10 changes: 5 additions & 5 deletions R/GAPIT.EMMAxP3D.R
View file
Original file line number Diff line number Diff line change
Expand Up @@ -822,13 +822,13 @@
}
#CV.Extragenetic specified
# beta.Extragenetic=beta
XCVI=XCV[,-c(1:(1+CV.Extragenetic)),drop=FALSE]
if(ncol(XCV)>1)XCVI=XCV[,-c(1:(1+CV.Extragenetic)),drop=FALSE]
XCVN=XCV[,c(1:(1+CV.Extragenetic)),drop=FALSE]
beta.I=as.numeric(beta)[-c(1:(1+CV.Extragenetic))]
if(ncol(XCV)>1)beta.I=as.numeric(beta)[-c(1:(1+CV.Extragenetic))]
beta.N=as.numeric(beta)[c(1:(1+CV.Extragenetic))]
BLUE.N=XCVN%*%beta.N
BLUE.I=rep(0,length(BLUE.N))
if(CV.Extragenetic!=0)BLUE.I=XCVI%*%beta.I
if(if(ncol(XCV)>1))BLUE.I=XCVI%*%beta.I
#Interception only
# if(length(beta)==1)XCV=X
BLUE=cbind(BLUE.N,BLUE.I)
Expand Down Expand Up @@ -871,9 +871,9 @@
# beta.Extragenetic=beta
if(!is.null(CV.Extragenetic))
{
XCVI=XCV[,-c(1:(1+CV.Extragenetic)),drop=FALSE]
if(ncol(XCV)>1)XCVI=XCV[,-c(1:(1+CV.Extragenetic)),drop=FALSE]
XCVN=XCV[,c(1:(1+CV.Extragenetic)),drop=FALSE]
beta.I=as.numeric(beta)[-c(1:(1+CV.Extragenetic))]
if(ncol(XCV)>1)beta.I=as.numeric(beta)[-c(1:(1+CV.Extragenetic))]
beta.N=as.numeric(beta)[c(1:(1+CV.Extragenetic))]
# print(is.null(beta.I))
BLUE.I=rep(0,nrow(XCVI))
Expand Down

0 comments on commit 6fd7265

Please sign in to comment.