问题标签 [asreml]

我正在使用asreml-R包进行分析。每当我存储一个拟合模型对象时asreml-R，都会提供一些我想隐藏的额外信息。请参阅以下代码和信息：

如果您帮助隐藏这些额外信息，我将不胜感激。请记住，这 asreml-R不是开源的。谢谢

R下面给出的是使用包分析可解析的 alpha 设计（alpha lattice design）的代码asreml。

我正在尝试使用Rpackage复制上述内容lme4。

如何计算vblup和vblue（差异的平均方差）lme4相当于？predict()$avsed ^ 2asreml

有谁知道对于缺少响应变量但所有固定效应值和随机效应因子水平已知的情况，是否有一种简单的方法可以从 ASreml-R 模型生成预测。通用的“拟合”函数仅返回对具有观察响应的案例的预测。

我将用于从安装在例如 lmer 中的模型对新数据集生成预测响应的“预测”函数似乎不适用于 asreml。

我不确定 predict.asreml 函数会完全满足我的要求，但也许我误解了如何使用它。

由于 ASReml 不是开源的，我没有提供可重现的示例。

谢谢，帕特里克

我们可以从 R 中的混合模型计算 BLUP 的标准误差 (SE) 吗？我正在使用一个名为 ASReml 的自定义包，它可以根据预测值计算 SE。但是我不确定在哪里可以找到 BLUP 值的 SE。

到目前为止，我有这样的事情......

使用 asreml 来拟合混合模型。

然后我们要求预测值：

上面的这些 std.errors 是针对 predict.value 列的。有没有办法获得 BLUP 或随机效应列的标准误差？

随机效应的 BLUP 是：

我正在使用来自在线练习教程的数据集，代码可以在第 7 页底部找到（https://tomhouslay.files.wordpress.com/2017/02/indivvar_mv_tutorial_asreml.pdf）

在教程中，他们声明他们正在引入“trait”作为关键字来表示多变量模型，但是当我运行完全相同的代码时，我得到以下信息：

eval(parse(text = x), envir = data, enclos = asreml4Env) 中的错误：找不到对象“特征”。

哈吉斯练习 csv 文件可以从这里下载：https ://figshare.com/articles/Haggis_data_behavioural_syndromes/4702540

这是教程提供的代码，asreml 函数有什么变化吗？

我正在使用同父异母的血统（分配的个人、父系、母系和群体（种群））。当我运行模型时，我收到消息“列出了零 df 的项……”，当我要求对固定效应进行 wald 测试时，我的一个效应丢失了。为我的数据分配谱系如何影响可用于测试我的固定效应的自由度？下面是我的代码。“sacode”和“treatcode”是我感兴趣的固定效果。当我运行模型并生成 Wald 测试时，“sacode”不包括在内。我在 R 工作室中使用 asreml-R4。

谱系示例：

示例代码：

我尝试添加和删除不同的变量，以查看是否有任何特定变量造成问题，但无论模型中的变量组合是什么，都会出现 zerodf 消息。

对于我的博士学位，我目前正在使用以下命令在 AS Reml R 中运行双变量方差分量估计：

但每次我尝试运行它时，我都会收到以下错误消息：

错误levels <- (*tmp*, value = as.character(levels))：因子级别 [4] 重复

我的固定因子（Betriebsnummer和gestLength）halfyear都被定义为因子，我还将它们的属性设置为“唯一”。我也用我的随机术语做到了这一点，SampleID是一个因素和Gs我的遗传关系矩阵。phenW是我的数据集，它是一个数据框。Reflex并且lakt_1是我的响应变量。

有人可以帮助我吗？

我正在尝试通过 R 运行一个 .bat 文件，其中包含 ASReml 的输入文件。.bat 文件在 Windows cmd 提示符下执行时可以正常工作，但我无法通过 R 使其工作。我使用的代码以前（几个月前）可以工作，但现在它给了我错误。

.bat 文件如下所示：

我用来运行 batfile 的代码是这样的：

错误代码是这样的：

有没有人建议我做错了什么？非常感谢！