计算残差R中

Question

这是我对我的矩阵代码：

>physician<-matrix(c(6,22,42,23,61,127,2,3,54,1,12,43,1,12,31),byrow=T,ncol=3) 
>dimnames(physician)<-list(specialty=c("Internal","Surgery","Radiotherapy","Oncology","Gynecology"),surgery=c("R","CR","C")) 
>physician 

我想知道如何计算残差输出为我的基质，所以我可以推断出“专业”之间的关联的性质和“手术”。谢谢

Answer 1

如果我清楚地理解，你想要做的是计算专业手术的多项Logistic回归的残差。下面是使用包装NNET的解决方案：

install.packages("nnet") 
    library(nnet) 

我调换你的矩阵的计算方法，并创建了一列code_surgery：

phy=data.frame(code_surgery=as.factor(c("R","CR","C")),data.frame(t(physician))) 

在你的数据帧一起来看看：

str(phy) 

然后计算回归：

reg <- multinom(code_surgery ~ Internal + Surgery + Radiotherapy + Oncology + Gynecology, data=phy) 

您可以用看看这个回归（函数把C作为code_surgery的参照组）：

summary(reg) 

我让你分析和继续。提示：你可以得到合适的值与

reg$fitted.values 

希望它会有所帮助。