allergy$fitted <- fitted(allergy.fit.main.2int)
The documentation "?fitted.glm" says "Value: Fitted values extracted from the object 'x'." I didn't try this with your data. However, from this documentation and previous experience similar situations, it would appear that this command would give you the fitted values as a column of your original data.frame.
hope this helps. spencer graves
allan clark wrote:
Hi all
I have the following generalized linear problem.
In a study of allergic responses, patients arriving at a clinic in Groote Schuur hospital were tested for sensitivity to a number of substances. Three of these were moulds: Cladosporium (C), Alternaria (F) and Aspergillius (T). Their level of sensitivity was measured on the Rast Scale as: 0: not allergic 1: mildly allergic 2 or more: allergic
The data is supplied below. The analysis is fairly straight forward and I understand how R solves the problem. I've supplied a copy of the code in order to perform the analysis.
My question is: HOW DOES ONE CONVERT THE OUTPUT (FITTED VALUES) SUPPLIED BY R AND DISPLAY THEM IN A CONTINGENCY TABLE?
allergy<-read.table("c:/a.dat",header=T) attach(allergy) allergy.fit.main.2int<-glm(y~ .^2,family=poisson, data=allergy) fitted(allergy.fit.main.2int)
The fitted values are: 1 2 3 4 5 6 7 8 666.198041 25.858478 14.943481 60.568127 5.350317 5.081557 28.233832 3.791205 9 10 11 12 13 14 15 16 10.974962 32.170233 7.068334 1.761433 14.077680 7.039304 2.883016 7.752087 17 18 19 20 21 22 23 24 5.892362 7.355551 35.631726 14.073188 14.295086 7.354193 6.610380 11.035427 25 26 27 9.014081 12.316433 62.669487
The data for those interested is as follows:
> allergy t f c y 1 none none none 671 2 none none one 23 3 none none two 13 4 none one none 60 5 none one one 8 6 none one two 3 7 none two none 24 8 none two one 4 9 none two two 15 10 one none none 31 11 one none one 9 12 one none two 1 13 one one none 14 14 one one one 6 15 one one two 4 16 one two none 9 17 one two one 5 18 one two two 7 19 two none none 32 20 two none one 15 21 two none two 17 22 two one none 8 23 two one one 5 24 two one two 12 25 two two none 12 26 two two one 13 27 two two two 59
Regards
Allan
______________________________________________
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
