library(gmodels)
?CrossTable # by Marc Schwartz
> a<-67; b<-10; c<-79; d<-67;
> dft<- with(as.data.frame.table(as.table(c(a=77,b=10,c=79,d=67))),
data.frame(Var1 = rep(Var1,Freq)) )
There are probably more efficient ways to replicate entries. I am just
drawing a blank right now.
> dft$Rw <- with(dft, ifelse((Var1 == "a" | Var1 ==
"c"), 1, 2))
> dft$Cl <- with(dft, ifelse((Var1 == "a" | Var1 ==
"b"), 1, 2))
> table(dft$Rw, dft$Cl)
1 2
1 77 10
2 79 67
> CrossTable(dft$Rw,dft$Cl, prop.chisq=FALSE)
Cell Contents
|-------------------------|
| N |
| N / Row Total |
| N / Col Total |
| N / Table Total |
|-------------------------|
Total Observations in Table: 233
| dft$Cl
dft$Rw | 1 | 2 | Row Total |
-------------|-----------|-----------|-----------|
1 | 77 | 10 | 87 |
| 0.885 | 0.115 | 0.373 |
| 0.494 | 0.130 | |
| 0.330 | 0.043 | |
-------------|-----------|-----------|-----------|
2 | 79 | 67 | 146 |
| 0.541 | 0.459 | 0.627 |
| 0.506 | 0.870 | |
| 0.339 | 0.288 | |
-------------|-----------|-----------|-----------|
Column Total | 156 | 77 | 233 |
| 0.670 | 0.330 | |
-------------|-----------|-----------|-----------|
--
David Winsemius
On Jan 10, 2009, at 11:32 AM, Bhargab Chattopadhyay wrote:
> Hi,
>
> I want to set a make a contingency table which will look like this..
> The problem is that I can't set the table like the following.
>
>
> col1 col2 Total
> row1 a b n10
> rp1 rp2 100
> cp1 cp2 rtp1
> ---------------------------------------------
> row2 c d n20
> rp3 rp4 100
> cp3 cp4 rtp2
> ---------------------------------------------
> Total n01 n02 n
> ctp1 ctp2 100
> 100 100 100
>
>
> Suppose
>
> a<-67; b<-10; c<-79; d<-67;
> n<-a+b+c+d;
> n10<-a+b;
> n01<-a+c;
> n20<-n-n10;
> n02<-n-n01;
> rp1<-(a/n10)*100;
> rp2<-100-rp1;
> rp3<-(c/n20)*100;
> rp4<-100-rp3;
> cp1<-(a/n01)*100;
> cp2<-100-cp1;
> cp3<-(c/n02)*100;
> cp4<-100-cp3;
> rtp1<-(n10/n)*100;
> rtp2<-100-rtp1;
> ctp1<-(n01/n)*100;
> ctp2<-100-ctp1;
>
> After this I can't procced in a meaningful way. Can anyone help me
> out?
> Thanks in advance.
>
> Bhargab
>
>
>
> [[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.