From: eugen pircalabelu <eugen_pircalabelu_at_yahoo.com>

Date: Sun, 09 Mar 2008 12:07:52 -0700 (PDT)

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

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.
*