Displaying 1 result from an estimated 1 matches for "obscounts".
2007 May 23
0
Replicated LR goodness-of-fit tests, heterogeneity G, with loglm?
...hetg = function() {
#creating some fake data
tran=gl(10,10) #tran
are the random replicates
cov= factor(rep(1:10,10))
#10 levels of this factor for each level of tran
tleng = c(219, 312, 178, 322, 311, 242, 235, 235, 266, 193)
#weighting for
#each transect (of
each level of tran)
obscounts = rpois(100, 50)
#these are the observed response count data
expprop = rep(c(.1, .1, .1, .05, .15, .1, .1, .1, .1, .1),
10) #expected
# proportion of
counts within 10 levels within each level of tran
vv = tapply(obscounts, tran, sum) #get marginal sum of
counts
rr = rep(vv, each=10) #...