Удельное в группе сравнения для lsmeans

голоса
0

Я хотел бы ограничить ретроспективные контрасты рассчитывается с использованием пакета 'lsmeans'в R.

library(nlme) #for gls
model<-gls(time~ benchmark*gc*opt, method=REML,data=d)

lsmeans.d<-lsmeans(model, ~benchmark:gc:opt)
pairs(lsmeans.d)

Это выводит сравнение всех возможных комбинаций тестов х дс х опт.

Мой вопрос: как я могу указать в функции lsmeans только вычислить и выходные сравнения в пределах каждого уникального теста: Gc комбинаций?

Я только хочу знать, если их разница между sigifciant «и выключаться» лечения в рамках каждой комбинации тестов х дс. Я добавил красные линии, чтобы показать это.

Необработанные данные ниже:

введите

Пример из 2

d=expand.grid(obs=0:10, benchmark=c('antlr', 'bloat', 'chart', 'eclipse', 'fop', 'hsqldb', 'jython', 'luindex', 'lusearch', 'pmd', 'xalan'), gc=c('CopyMS', 'GenCopy', 'GenImmix', 'GenMS', 'Immix'), opt=c('on', 'off'), heapSize=seq(from=1.5, to=4, by=0.5))
d$time = rexp(nrow(d), 0.01)+1000
d$time = d$time + abs(d$heapSize-3)*100
d$time[d$opt=='on'] = d$time[d$opt=='on']-200

d$time[d$opt=='on' & d$benchmark=='bloat'] = d$time[d$opt=='on' & d$benchmark=='bloat'] + 190
d$time[d$opt=='on' & d$benchmark=='pmd' & d$gc=='Immix'] = d$time[d$opt=='on' & d$benchmark=='pmd' & d$gc=='Immix'] + 600

ggplot() +
facet_grid(gc~benchmark) +
geom_boxplot(data=d, mapping=aes(x=opt, y=time, color=opt))
Задан 20/10/2018 в 05:25
источник пользователем
На других языках...                            


1 ответов

голоса
0

Найдено ответов благодаря вдохновению здесь: первый ответ на РВС

Using the FIRST version of your 'fit' model, do lsm = lsmeans(fit, ~A*B|C) and then contrast(lsm, list(c = c(1,0,0,-1)) – rvl Nov 24 '16 at 21:33

Ответ:

lsm = lsmeans(model, ~benchmark*opt|gc)
pairs(lsm)
Ответил 20/10/2018 в 08:13
источник пользователем

Cookies help us deliver our services. By using our services, you agree to our use of cookies. Learn more