在R中使用加权数据进行方差分析和多重比较测试

时间:2019-02-13 16:07:57

标签: r anova weighted

这是我的问题的上下文:我有一个调查数据集,其中包含受访者,产品,喜欢的得分和“权重”列(根据他们的人口统计信息与每个受访者相关的权重)。 行是一个受访者对一种测试产品给出的喜好评分。每个受访者测试了13种产品(每个受访者有13条产品线)。共有500名消费者参加了测试。

目标:对加权数据执行ANOVA和多重比较测试以比较产品性能(在这种情况下,将执行LSD测试-即使这不是我所知的最佳方法)

在R中,我可以执行加权方差分析: model_lm = lm(liking~product+respondent, data = survey,weights=weighting)

并且根据我的加权数据可以很好地计算出系数。

但是在执行LSD测试时:

library(agricolae)
test = LSD.test(model_lm , trt='BRANDS' , alpha = 0.1,p.adj='none')

1 /每种产品计算的平均分数未加权 2 / LSD计算基于非加权数据

问题: 在R中,如何获得基于加权平均分数的多重比较测试LSD的结果?

也有一个统计程序可以根据 标准误差的权重:基于实际N的标准误差 而不是加权N。 see幻灯片10:

为什么?因为随着对数据进行加权会增加标准偏差,所以如果我不使用此过程,由于标准错误权重的存在,判别就更少了。

谢谢

摩根

0 个答案:

没有答案