從 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

我想知道如何計算矩陣的殘差輸出，以便我可以推斷specialty和surgery之間關聯的性質。 謝謝

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 

希望它會有所幫助。