Description Details Author(s) References Examples

Description: This package performs a Latent Class Analysis of phenotypic measurements in pedigrees and a model selection based on one of two methods: likelihood-based cross-validation and Bayesian Information Criterion. It computes also individual and triplet child-parents weights in a pedigree using an upward-downward algorithm. It takes into account the familial dependence defined by the pedigree structure by considering that a class of a child depends on his parents classes via triplet-transition probabilities of the classes. The package handles the case where measurements are available on all subjects and the case where measurements are available only on symptomatic (i.e. affected) subjects. Distributions for discrete (or ordinal) and continous data are currently implemented. The package can deal with missing data.

Package: | LCAextend |

Type: | Package |

Version: | 1.0 |

Date: | 2011-11-15 |

License: | GPL |

LazyLoad: | yes |

Arafat TAYEB <arafat.tayeb@crulrg.ulaval.ca>, Alexandre BUREAU <alexandre.bureau@msp.ulaval.ca> and Aurelie Labbe <aurelie.labbe@mcgill.ca>

Maintainer: Alexandre BUREAU <alexandre.bureau@msp.ulaval.ca>

TAYEB, A. LABBE, A., BUREAU, A. and MERETTE, C. (2011) Solving Genetic Heterogeneity in Extended
Families by Identifying Sub-types of Complex Diseases. Computational Statistics, 26(3): 539-560. DOI: 10.1007/s00180-010-0224-2,

LABBE, A., BUREAU, A. and MERETTE, C. (2009) Integration of Genetic Familial Dependence Structure in Latent Class Models. The International Journal of Biostatistics, 5(1): Article 6.

1 2 3 4 5 6 7 8 9 | ```
#data
data(ped.cont)
fam <- ped.cont[,1]
#the function applied for the two first families of ped.cont
res <- model.select(ped.cont[fam%in%1:2,],distribution="normal",trans.const=TRUE,
optim.indep.norm,optim.probs.indic=c(TRUE,TRUE,TRUE,TRUE),
famdep=TRUE,selec="bic",K.vec=1:3,tol=0.001,x=NULL,var.list=NULL)
plot(1:3,res$bic,type="l",col="blue",xlab="K",ylab="BIC",main="model
selection using BIC")
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.