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@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.
______________________________________________
R-help@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.
