Вычисление невязки из матрицы в 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

Если я четко понял, что вы хотите сделать, это вычислить остатки от многокомпонентной логистической регрессии хирургии по специальности. Вот решение с использованием пакета Nnet:

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

Я транспонированная свой Matrice и создал столбец 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 

Надеюсь, это поможет.