我正在为列联表拟合多项式模型。我想提取7x8表格的预期单元格计数。我不明白如何从R中的multinom()
包中的nnet
函数中获得预期的单元格计数?
这是我的代码:
counts <- c(5,10,25,32,45,201,37,29,110,144,147,152,39,190,28,78,38,20,54,78,54,189,175,30,26,89,74,30,48,56,71,82,34,75,49,103,156,174,19,60,31,145,175,182,56,82,90,39,50,38,54,87,99,99,67,82)
outcome <- gl(8,1,56)
treatment <- gl(7,8)
mod <- multinom(counts~outcome+treatment)
我需要从模型中获得预期的细胞计数。我不了解从模型摘要中获得的拟合值。我不知道拟合值的哪一栏代表预期的细胞计数。