-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathshelltox.R
83 lines (65 loc) · 3.11 KB
/
shelltox.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
#######################################################################
####### Dose group control and High ####################################
rm(list=ls())
library(CorrBin)
library(vcd)
data(shelltox)
shelltox=shelltox
head(shelltox)
#install.packages("CorrBin")
summary(shelltox)
ftable(xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox))
#####################################
##### 0 response for c,l,m,h#########
rbind(xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,1][1,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,2][1,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,3][1,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,4][1,])
#####################################
##### 1 response for c,l,m,h#########
rbind(xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,1][2,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,2][2,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,3][2,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,4][2,])
#####################################
##### 2 response for c,l,m,h#########
rbind(xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,1][3,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,2][3,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,3][3,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,4][3,])
#####################################
##### 3 response for c,l,m,h#########
rbind(xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,1][4,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,2][4,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,3][4,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,4][4,])
#####################################
##### 4 response for c,l,m,h#########
rbind(xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,1][5,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,2][5,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,3][5,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,4][5,])
#####################################
##### 5 response for c,l,m,h#########
rbind(xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,1][6,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,2][6,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,3][6,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,4][6,])
#####################################
##### 6 response for c,l,m,h#########
rbind(xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,1][7,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,2][7,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,3][7,],
xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,4][7,])
#####################################
##### structable performs the task above for c,l,m,h#########
structable(xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox))
#####################################
##### selects a 0 response for c,l,m,h#########
structable(xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox))[1,]
#####################################
##### selecting response for c,h #########
data=structable(xtabs(Freq~NResp+ClusterSize+Trt,data=shelltox)[,,c(1,4)])
dim(data)[1]
dim(data)[2]
dim(data)[3]