Displaying 2 results from an estimated 2 matches for "expcounts".
Did you mean:
expcount
2019 Oct 07
6
[nbdkit PATCH 0/5] More retry fixes
I think this is my last round of patches for issues I identified with
the retry filter. With this in place, it should be safe to interject
another filter in between retry and the plugin.
Eric Blake (5):
retry: Don't call into closed plugin
tests: Refactor test-retry-reopen-fail.sh
tests: Enhance retry test to cover failed reopen
server: Move prepare/finalize/close recursion to
2007 May 23
0
Replicated LR goodness-of-fit tests, heterogeneity G, with loglm?
...d 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) #rep the marginal sums across each
cell
expcounts = expprop*rr #the vector of expected counts
#Now calculate the G for each level of tran (likelihood
ratios)
G = vector()
for (i in levels(tran)) {
obsi =
obscounts[tran==i]
expi =
expcounts[tran==i]
G[i] =
2*sum((obsi*log(obsi/expi)))
}
dfs = rep(max(as.integer(leve...