計算殘差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 

希望它會有所幫助。