require(GenABEL.data)
data(srdta)
class(srdta)
x <- srdta@gtdata
class(x)
nids(x)
nsnps(x)
idnames(x)[1:12]
male(x)[1:12]
male(x)[c("p1","p2","p3","p4")]
snpnames(x)[1:4]
chromosome(x)[1:4]
map(x)[1:4]
n4 <- c("rs18","rs655")
n4
map(x)[n4]
n4 <- c("rs18","rs65")
n4
map(x)[n4]
chromosome(x)[n4]
x[1:12,1:4]
summary(x[,1:10])
as.numeric(x[1:12,1:4])
as.numeric(x[c("p1","p3","p4"),c("rs18","rs65")])
as.character(x[c("p1","p3","p4"),c("rs18","rs65")])
as.genotype(x[c("p1","p3","p4"),c("rs18","rs65")])
as.hsgeno(x[c("p1","p3","p4"),c("rs18","rs65")])
Run the code above in your browser using DataLab