create.data(varbvs)
create.data()所属R语言包:varbvs
Simulate SNP data
模拟SNP数据
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Generates samples of the genotypes and quantitative trait (i.e. a continuous outcome Y) according to the specified SNP minor
生成样本的基因型和数量性状的(即连续的结果Y),根据指定的SNP小调
用法----------Usage----------
create.data(maf, beta, sigma, n)
参数----------Arguments----------
参数:maf
vector giving minor allele frequencies of the SNPs.
矢量给未成年人的SNP等位基因频率。
参数:beta
vector of the same length giving the additive effects of the SNPs.
给予的单核苷酸多态性的累加效应相同的长度的矢量。
参数:sigma
positive scalar specifying the variance of residual.
正标量指定剩余的方差。
参数:n
number of samples to generate.
要生成的样本数。
Details
详细信息----------Details----------
Genotypes are generated from a binomial distribution with success rates given by the minor allele frequencies. Observations about the quantitative trait are generated according to y = X*beta + e, where the residual e is normal with mean zero and
基因型产生的二项式分布的次要等位基因频率的成功率。有关的数量性状的观测产生的根据y = X*beta + e的,其中的剩余e是正常的均值为零,
值----------Value----------
Returns a list containing two components:
返回一个列表,其中包含两个组成部分:
参数:X
matrix of genotype data with centered columns so that the mean of each column is zero.
与中心的列,以使每一列的平均值为零的基因型数据矩阵。
参数:y
vector of quantitative trait data, centered so that the mean of the vector entries is zero.
向量的数量性状数据为中心,使得向量项的平均值是零。
(作者)----------Author(s)----------
Peter Carbonetto
参见----------See Also----------
create.snps
create.snps
实例----------Examples----------
## Randomly generate genotypes and quantitative trait measurements for[#随机生成基因型和数量性状的测量]
## 1000 SNPs and 200 individuals, in which the variance of the residual[#1000 SNPs和200个个体,其中所述的残余方差]
## error is 4. Of these SNPs, 10 have nonzero additive effects on the[#误差为4。这些SNPs中,有非零的加性效应]
## trait.[#特征。]
se <- 4
snps <- create.snps(1000,10)
data <- create.data(snps$maf,snps$beta,se,200)
## Calculate the proportion of variance in the quantitative trait[#计算方差的比例在数量性状]
## explained by the SNPs.[#解释的单核苷酸多态性。]
sz <- var(c(data$X %*% snps$beta))
sz/(sz + se)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|