我试图绘制病毒载量和试验臂的简单箱线图,但不断收到此错误“x[floor(d)] + x[ceiling(d)] 中的错误:二元运算符的非数字参数”。
数据集
PatientID trial_arm viral_load
1 club 19
2 clinic 19
3 club 19
4 club 64
5 clinic 32
6 clinic 29
7 clinic 92
8 club 19
9 clinic 19
10 club 27
Run Code Online (Sandbox Code Playgroud)
当我运行代码时:
> boxplot(dataset,viral_load~trial_arm)
Error in x[floor(d)] + x[ceiling(d)] :
non-numeric argument to binary operator
Run Code Online (Sandbox Code Playgroud)
Trial_arm 变量设置为一个因子,而 Viral_load 是数字。
如何在 R 中拟合序数(3 个级别)逻辑混合效应模型?我想它就像一个 glmer,除了三个结果级别之外。
数据结构
patientid Viral_load Adherence audit_score visit
1520 0 optimal nonhazardous 1
1520 0 optimal nonhazardous 2
1520 0 optimal hazardous 3
1526 1 suboptimal hazardous 1
1526 0 optimal hazardous 2
1526 0 optimal hazardous 3
1568 2 suboptimal hazardous 1
1568 2 suboptimal nonhazardous 2
1568 2 suboptimal nonhazardous 3
Run Code Online (Sandbox Code Playgroud)
其中病毒载量(感兴趣的结果)由三个级别(0、1、2)组成,依从性 - 最佳/次优、审核评分无危险/危险以及 3 次就诊。
这是使用广义混合效果模型代码的模型外观的示例。
library (lme4)
test <- glmer(viral_load ~ audit_score + adherence + (1|patientid) + (1|visit), data = df,family = binomial) …Run Code Online (Sandbox Code Playgroud)