[R] GLM FITTED VALUES TABLE
Spencer Graves
spencer.graves at pdf.com
Fri Nov 28 17:59:54 CET 2003
Did you try something like the following:
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
>______________________________________________
>R-help at stat.math.ethz.ch mailing list
>https://www.stat.math.ethz.ch/mailman/listinfo/r-help
>
>
More information about the R-help
mailing list