haplo.cc.match(SimHap)
haplo.cc.match()所属R语言包:SimHap
Haplotype analysis for matched case-control data
匹配的情况下,控制数据的单倍型分析
译者:生物统计家园网 机器人LoveR
描述----------Description----------
haplo.surv performs a series of conditional logistic regression models to matched case-control data with haplotypes using a simulation-based approach to account for uncertainty in haplotype assignment when phase is unknown.
haplo.surv进行了一系列的条件Logistic回归模型匹配的情况下,控制数据,使用模拟为基础的方法来解释阶段是未知的单倍型分配的不确定性时的单倍型。
用法----------Usage----------
haplo.cc.match(formula1, formula2, pheno, haplo, sim, effect = "add",
sub = NULL)
参数----------Arguments----------
参数:formula1
a symbolic description of the full model to be fit, including haplotype parameters. The response must be binary indicator of case-control status, and the formula must contain a variable indicating strata, or the matching sequence.
一个象征性的描述,完整的模型是合适的,包括单倍型参数。响应必须是二进制的情况下,控制状态的指标,计算公式必须包含一个变量,表示阶层,或匹配序列。
参数:formula2
a symbolic description of the nested model excluding haplotype parameters, to be compared to formula1 in a likelihood ratio test. The response must be binary indicator of case-control status, and the formula must contain a variable indicating strata, or the matching sequence.
一个象征性的嵌套模式,不包括单倍型参数的说明,相比formula1似然比检验。响应必须是二进制的情况下,控制状态的指标,计算公式必须包含一个变量,表示阶层,或匹配序列。
参数:pheno
a dataframe containing phenotype data.
一个数据框的表型数据。
参数:haplo
a haplotype object made by make.haplo.rare.
单倍型对象由make.haplo.rare。
参数:sim
the number of simulations from which to evaluate the results.
的数量的模拟来对结果进行评估。
参数:effect
the genetic effect type: "add" for additive, "dom" for dominant and "rec" for recessive. Defaults to additive. See note.
的遗传效应类型:"add"添加剂"dom"主导和"rec"为隐性。默认至添加剂。见注。
参数:sub
optional. An expression using a binary operator, representing a subset of individuals on which to perform analysis. e.g. sub=expression(sex==1).
可选的。表达式中使用一个二元操作符,代表个人的一个子集进行分析。例如sub=expression(sex==1)。
Details
详细信息----------Details----------
formula1 should be in the form:
“formula1的形式应该是:
值----------Value----------
haplo.cc.match returns an object of 'class' hapClogit.
haplo.cc.match返回一个对象的“类”hapClogit。
The summary function can be used to obtain and print a summary of the results.
summary函数可以用来获取和打印结果的摘要。
An object of class hapClogit is a list containing the following components:
一个对象的类hapClogit的是一个列表,其中包含以下组件:
参数:formula1
formula1 passed to haplo.cc.match.
formula1传递给haplo.cc.match。
参数:formula2
formula2 passed to haplo.cc.match.
formula2传递给haplo.cc.match。
参数:results
a table containing the odds ratios, confidence intervals and p-values of the parameter estimates, averaged over the n=sim models performed.
一个表,其中包含的胜算比,置信区间和参数估计的p值,平均在n = sim模型进行。
参数:empiricalResults
a list containing the odds ratios, confidence intervals and p-values calculated at each simulation
一个列表,其中包含的胜算比,在每个模拟计算的置信区间和p-值
参数:logLik
the average log-likelihood for the n=sim linear models fit using formula1.
平均对数似然为n =sim线性模型适合使用formula1。
参数:LRT
a likelihood ratio test, testing for significant improvement of the model when haplotypic parameters are included
似然比检验,显着改善的模型进行测试时,单倍型参数
参数:ANOVA
analysis of variance, comparing the two models fit with and without haplotypic parameters.
方差分析,比较这两个模型并没有合适的单倍型参数。
参数:Wald
The Wald test for overall significance of the fitted model including haplotypes.
Wald检验整体意义的拟合模型,包括单倍型。
参数:rsquared
r-squared values for models fit using formula1 and formula2.
模型的R平方值适合使用formula1和formula2。
参数:effect
the haplotypic effect modelled, "ADDITIVE", "DOMINANT" or "RECESSIVE".
单倍型效果为蓝本,“添加剂”,“显性”或“隐性”。
注意----------Note----------
To model a codominant haplotypic effect, define the desired haplotype as a factor in the formula1 argument. e.g. factor(h.AAA), and use the default option for effect.
要对一个显性单倍型的效果,定义所需的单倍型的一个因素formula1参数。例如factor(h.AAA),并使用默认选项为effect。
(作者)----------Author(s)----------
Pamela A. McCaskie
参考文献----------References----------
参见----------See Also----------
snp.cc.match, haplo.bin
snp.cc.match,haplo.bin
实例----------Examples----------
data(SNP.dat)
# convert SNP.dat to format required by infer.haplos[SNP.dat转换,格式化所需的infer.haplos]
haplo.dat <- SNP2Haplo(SNP.dat)
data(pheno.dat)
newdata <- prepare.cc(geno=haplo.dat, pheno=pheno.dat, cc.var="DISEASE")
newhaplo.dat <- newdata$geno
newpheno.dat <- newdata$pheno
# generates haplotype frequencies and haplotype design matrix[产生单倍型频率和单倍型设计矩阵]
myinfer<-infer.haplos.cc(geno=newhaplo.dat,
pheno=newpheno.dat, cc.var="DISEASE")
# prints haplotype frequencies among cases[打印单倍型频率在病例]
myinfer$hap.freq.cases
# prints haplotype frequencies among controls[打印单体型频率在控制]
myinfer$hap.freq.controls
# generate haplo object where haplotypes with a frequency [单倍型的频率产生单倍体对象]
# below min.freq are grouped as a category called "rare"[下面min.freq分为一类被称为“罕见”]
myhaplo<-make.haplo.rare(myinfer,min.freq=0.05)
mymodel <- haplo.cc.match(formula1=DISEASE~SBP+DBP+h.N1AA+strata(STRAT),
formula2=DISEASE~SBP+DBP+strata(STRAT), haplo=myhaplo,
pheno=pheno.dat, sim=10)
summary(mymodel)
# example using a subsetting variable - looking at males only[例如,使用一个子集变量 - 只在男性]
mymodel <- haplo.cc.match(formula1=DISEASE~SBP+DBP+h.N1AA+strata(STRAT),
formula2=DISEASE~SBP+DBP+strata(STRAT), haplo=myhaplo,
pheno=pheno.dat, sim=10, sub=expression(SEX==1))
summary(mymodel)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|