Re: [R] question for crosstable

From: eugen pircalabelu <eugen_pircalabelu_at_yahoo.com>
Date: Sun, 09 Mar 2008 12:07:52 -0700 (PDT)


Yes, Mark! Thank you!
I was repeatedly using the function with two vectors and i have forgotten that it could be used also with a matrix. Thank you once again!

Marc Schwartz <marc_schwartz_at_comcast.net> wrote: eugen pircalabelu wrote:
> Good evening R-users!
>
> I have the following problem: i want to get a weighted crosstable along with the adjusted standardized residuals test
>
> Example:
> a<-c(1,1,1,1,1,2,2,2,3,3)
> b<-c(10,10,23,24,33,45,46,70,21,66)
> c<-c(3,3,2,3,4,1,1,1,3,3)
> d<-c("a", "b","b","c","a","a","a", "b", "c","c")
> z<-data.frame(a, b,c,d)
> library(gmodels)
>
> Using the above example i want to get this :
>
> xtabs(c ~ a+d, z) , where i have
> sum(xtabs(c ~ a+d, z) ) as my total number of cases, but with the adjusted standardized residuals test.
>
> So i want something like the CrossTable function from gmodels package, but taking into account that i have weights which sum to 24 and not 10, as the CrossTable prints from this:
>
> CrossTable(z$a, z$d, format="SPSS", asresid=T)
>
> Thank you and have a great day ahead!

Is this what you want?

 > CrossTable(xtabs(c ~ a + d, z), format = "SPSS", asresid = TRUE)

    Cell Contents

|-------------------------|
|                   Count |
| Chi-square contribution |
|             Row Percent |
|          Column Percent |
|           Total Percent |
|           Adj Std Resid |
|-------------------------|

Total Observations in Table: 24

              | d
            a |        a  |        b  |        c  | Row Total |
-------------|-----------|-----------|-----------|-----------|
            1 |        7  |        5  |        3  |       15  |
              |    0.336  |    0.417  |    1.225  |           |
              |   46.667% |   33.333% |   20.000% |   62.500% |
              |   77.778% |   83.333% |   33.333% |           |
              |   29.167% |   20.833% |   12.500% |           |
              |    1.198  |    1.217  |   -2.286  |           |
-------------|-----------|-----------|-----------|-----------|
            2 |        2  |        1  |        0  |        3  |
              |    0.681  |    0.083  |    1.125  |           |
              |   66.667% |   33.333% |    0.000% |   12.500% |
              |   22.222% |   16.667% |    0.000% |           |
              |    8.333% |    4.167% |    0.000% |           |
              |    1.116  |    0.356  |   -1.434  |           |
-------------|-----------|-----------|-----------|-----------|
            3 |        0  |        0  |        6  |        6  |
              |    2.250  |    1.500  |    6.250  |           |
              |    0.000% |    0.000% |  100.000% |   25.000% |
              |    0.000% |    0.000% |   66.667% |           |
              |    0.000% |    0.000% |   25.000% |           |
              |   -2.191  |   -1.633  |    3.651  |           |
-------------|-----------|-----------|-----------|-----------|
Column Total |        9  |        6  |        9  |       24  |
              |   37.500% |   25.000% |   37.500% |           |
-------------|-----------|-----------|-----------|-----------|

Note that the first argument to CrossTable(), 'x', can be either a vector or a matrix, the latter if 'y' is unspecified.

HTH, Marc Schwartz        


        [[alternative HTML version deleted]]



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Sun 09 Mar 2008 - 19:10:02 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sun 09 Mar 2008 - 19:30:26 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive