问题描述
我正在尝试在R中运行Poisson模型。我检查了许多网站来做到这一点，并且我相信下面的代码应该可以工作。但是，我收到一条错误消息，指出有一个未使用的参数。我不明白问题是什么。
library(lme4)
m1<-lmer(sum ~ Lunch + intervention_y_n + (1 | prcid),data=data_long_1,family = "poisson"(link = "log"))
lmer中的错误（总和〜午餐+干预_y_n +（1 | prcid），数据= data_long_1，： 未使用的参数（family = poisson（link =“ log”））
解决方法
您需要data raw;
call streaminit(123);
do id = 1 to 100;
do t = 0 to rand('integer',1,2);
do rep = 1 to rand('integer',20);
y = round(rand('uniform'),0.01);
output;
end;
end;
end;
run;
proc sql;
create table want as
select id,t,avg
from (
select id,mean(y) as avg format=5.2
from raw
group by id,t
) as averages
group by id
having avg = max(avg)
;
用于广义线性混合模型（例如，泊松模型）。 glmer()仅适用于线性混合模型（带有身份链接的高斯响应），因此它没有lmer()自变量。
此代码一直有效到最近（4个月前）。如果您检查NEWS file for lme4，将会看到2020年4月发布的1.1-22版本禁用了将family自动转发到lmer(...,family=...)的功能。但是，从2013年开始，此用法已被弃用（即发出警告消息，告诉您不要使用它）。