On this Picostat.com statistics page, you will find information about the l51 data set which pertains to An example pedigree data. The l51 data set is found in the gap R package. You can load the l51 data set in R by issuing the following command at the console data("l51"). This will load the data into a variable called l51. If R says the l51 data set is not found, you can try installing the package by issuing this command install.packages("gap") and then attempt to reload the data. If you need to download R, you can go to the R project website. You can download a CSV (comma separated values) version of the l51 R data set. The size of this file is about 966 bytes.

An example pedigree data

Description

The data contains data on 51 individuals in a pedigree. Below it is used for comparing results from various packages.

Usage

data(l51)

Format

A data frame

Source

Morgan v3.

References

Morgan v3. http://www.stat.washington.edu/thompson/Genepi/MORGAN/Morgan.shtml

Examples

## Not run: 
km <- kin.morgan(l51)
k2 <- km$kin.matrix*2# quantitative trait
library(regress)
r <- regress(qt ~ 1, ~k2, data=l51)
names(r)
r
# qualitative trait
N <- dim(l51)[1]
w <- with(l51,quantile(qt,probs=0.75,na.rm=TRUE))
ped51 <- within(l51, bt <- ifelse(qt<=w,0,1))
d <- regress(bt ~ 1, ~k2, data=ped51)
d
# for other tests not shown here
set.seed(12345)
ped51 <- within(ped51,{r <- rnorm(N); bt[is.na(bt)] <- 0})
library(foreign)
write.dta(ped51,"ped51.dta")## End(Not run)

--

Dataset imported from https://www.r-project.org.