# NOT RUN {
## LOAD DATA ##
data(HorseFlu)
## EXAMINE CONTENT ##
summary(HorseFlu)
## individual info
head(HorseFlu@individuals)
## DNA sequences
HorseFlu@dna
## records info
lapply(HorseFlu@records, head)
## How many individuals and sequences?
get.nindividuals(HorseFlu)
get.nsequences(HorseFlu)
## How many sequences per individual?
ind <- table(get.data(HorseFlu, "individualID", where="dna"))
ind
barplot(sort(ind), horiz=TRUE, las=1,
xlab="number of samples", cex.names=.8)
## How many sequences for this individual?
ind.42 <- subset(HorseFlu, individualID="42")
get.nsequences(ind.42)
## How many samples?
length(unique(get.data(HorseFlu, "sampleID", where="dna")))
## How many sequences per sample?
table(get.data(HorseFlu, "sampleID", where="dna"))
# }
Run the code above in your browser using DataLab