### R代码来自Vignette源'tdt-vignette.rnw'############################################################### ###代码块编号1：家庭数据########################################################################################误差############################################ mis < - 误解（data = peddata，snp.data =基因型）dim（mis）#############################################################################，na.rm = true）per.snp < - applion（mis，2，sum，na.rm = true）par（mfrow = c（1，2））hist（per.subj，main ='main ='tocture perakectam persuctive'，xlab ='主题'）hist（per.snp，main ='直方图per snp'，xlab ='snp'）#########################################################################################“ familyid”]体育r.fam < - tapply（per.subj，fam，sum）par（mfrow = c（1，1））hist（per.fam，main ='finastarge pernive per family'，xlab ='family'）################################代码块数字5：TDT测验########################################################测试<-tdt.snp（data = peddata，snp.data =基因型）cbind（p.values.1df = p.value（tests，1），p.values.2df = p = p = p = p.value（测试，2））qq.chisq（chi.squared（tests，1），df = 1）