Kinship model with standard statistics

以前面建立的 "利用自定义链接建立亲属关系模型 "为例,增加几个统计数据,即:

小学生
成年男性
成年女性和老年人的数量
以及: 一个家庭的平均子女数

1、打开上次建立的亲属关系模型(Anylogic 创建一个简单的基于代理的人口模型

2、设置代理展示的位置

打开 Main 编辑页面,添加一个 rectangle (直接拖拽过来),名字设置为 people_presentation
并设置一下相关数值
fill color: new Color(221, 160, 221, 32)
width: 520
height: 520

在这里插入图片描述

3、为 the replicated object 创建统计项目

点击一下 people 图标,在右边出现的属性页面中打开 statistics 页面

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
在R语言中,可以使用`rrBLUP`包进行基因组选择的预测。假设我们有一个包含3个品种、每个品种包含8个SNP的基因型数据,格式为Snp(agct),其中1表示该位点为变异位点,0表示该位点为野生位点。现在需要将基因型数据进行独热编码,并使用`rrBLUP`包进行预测。以下是一个示例代码: ``` # 加载rrBLUP包 library(rrBLUP) # 创建基因型数据 genotype <- matrix(c( 1, 0, 1, 0, 1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 0, 0, 1, 1, 0, 0 ), nrow = 3, byrow = TRUE) # 将基因型数据进行独热编码 genotype_one_hot <- apply(genotype, 2, function(x) ifelse(x == 1, "A", "B")) genotype_one_hot <- data.frame(genotype_one_hot) genotype_one_hot <- model.matrix(~., data = genotype_one_hot)[,-1] # 创建表型数据 phenotype <- c(1.2, 0.8, 1.5) # 进行基因组选择的预测 model <- kin.blup(phenotype ~ genotype_one_hot, K = kinship(genotype_one_hot)) prediction <- predict(model, genotype_one_hot) ``` 以上代码中,首先需要加载`rrBLUP`包。然后,创建一个包含基因型数据的矩阵`genotype`,其中1表示该位点为变异位点,0表示该位点为野生位点。接着,将基因型数据进行独热编码,并保存在`genotype_one_hot`变量中,同时创建一个包含表型数据的向量`phenotype`。最后,使用`kin.blup()`函数建立基因组选择模型,并使用`predict()`函数对基因型数据进行预测,预测结果保存在`prediction`变量中。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值