Dear List,

 

I would like to calculate the Goodman-Kruskal gamma for the predicted
classes obtained from an ordinal regression model using lrm in the Design
package. I couldn't find a way to get gamma for predicted values in Design
so have found previous positings suggesting to use :

 

Rcorr.cens(x, S outx = TRUE)  in the Hmisc package

 

My question is, will this work for predicted vs observed factors?  I.e. x =
predicted class and S = observed class? Or is there a better way to obtain
this? I used the maximum individual probability for each observation to
determine the predicted class. 

 

Any help appreciated, 

 

Thanks

 

Colin

 

 

 

Colin Robertson

Dept of Geography

University of Victoria

 

 


	[[alternative HTML version deleted]]

