data(pedigree)
data(pheno)
data(geno)
foo <- gls.batch.get(
phenfile=pheno,
genfile=data.frame(t(geno)),
pedifile=pedigree,
outfile="example_output.txt",
covmtxfile.in=NULL,covmtxfile.out=paste(phen,"_cov_matrix.txt",sep=""),
phen="Zscore", covars = "IsFemale",
med = "rfgls", sizeLab = "OOPP", Mz = TRUE, Bo = TRUE, Ad = TRUE, Mix = TRUE,
indobs = TRUE, col.names = TRUE, pediheader = FALSE,
pedicolname=c("FAMID","ID","PID","MID","SEX"),
sep.phe = "", sep.gen = "", sep.ped = "")
str(foo)
##Also works:
pedigree2 <- pedigree
pedigree2[,-2] <- NA ##<--Change all but column 'ID' to NA.
foo2 <- gls.batch.get(
phenfile=pheno,
genfile=data.frame(t(geno)),
pedifile=pedigree2, ##<--Note change.
outfile="example_output.txt",
covmtxfile.in=NULL,covmtxfile.out=paste(phen,"_cov_matrix.txt",sep=""),
phen="Zscore", covars = "IsFemale",
med = "rfgls", sizeLab = "OOPP", Mz = TRUE, Bo = TRUE, Ad = TRUE, Mix = TRUE,
indobs = TRUE, col.names = TRUE, pediheader = FALSE,
pedicolname=c("FAMID","ID","PID","MID","SEX"),
sep.phe = "", sep.gen = "", sep.ped = "")
Run the code above in your browser using DataLab