Well, im trying to fit a model with 5 "classes" Controls - DiseaseA_soft - DiseaseA_hard - DiseaseB_soft - DiseaseB_hard.
Being soft and hard 2 levels of their prognosis. "Hard" is way severe than "soft". Anyway, i don think a classical "Class ~." model should be used, cause we got 2 sub-ordinal levels for each disease. ¿Is there a way i can correctly represent this "design" with nnet::multinom?
For a 5 classes (clase = 1 column with 5 levels of a factor unordered) is just
multinom(clase~variables, data= data_training, maxit=10000 )
I dont know if i should make 2 columns, 1 with a 3 lvl factor determining Controls, Dis_A, and Dis_B, and the other with the severity of the disease. (Mostly because i dont know what to do with the Controls, i dont think is correct to treat them as the "gentle" version of the disease.
¿Can i do something like this? ¿And that will be correct?
multinom(disease_type+severity~variables, data= data_training, maxit=10000 )