之前做两列数据学会了用glm的binomial 拟合,然后用multcomp的glht做多组比较显著性,现在有一组数据3列,glm没有多项的模型,查了一下用glmnet似乎可以,但是之后如何比较显著性呢?如下面的数据,用glmnet, cbind('A','B','C')~crop multinomial拟合,如何比较wheat和rice这两组的显著性?另外glmnet出来的model也没有变量对线性的显著性?谢谢!
| crop | pop | A | B | C |
| wheat | Sichuan | 75 | 8 | 1 |
| wheat | Hubei | 4 | 50 | 3 |
| wheat | Anhui | 10 | 81 | 1 |
| wheat | Jiangsu | 5 | 60 | 0 |
| wheat | Fujian | 147 | 23 | 41 |
| rice | Sichuan | 178 | 0 | 0 |
| rice | Hubei | 21 | 36 | 3 |
| rice | Anhui | 31 | 78 | 0 |
| rice | Jiangsu | 28 | 108 | 4 |
| rice | Fujian | 61 | 0 | 53 |