Displaying 2 results from an estimated 2 matches for "newxi".
Did you mean:
newx
2007 May 21
1
Boostrap p-value in regression [indirectly related to R]
...model
mod <- lm(yi ~ xi)
summary(mod)
b1 <- coef(mod)[2]
t1 <- coef(mod)[2] / coef(summary(mod))[2,2]
### 1000 bootstrap replications using (X,Y)-pair resampling
t1.star <- rep(NA,1000)
for (i in 1:1000) {
ids <- sample(1:n, replace=TRUE)
newyi <- yi[ids]
newxi <- xi[ids]
mod <- lm(newyi ~ newxi)
t1.star[i] <- ( coef(mod)[2] - b1) / coef(summary(mod))[2,2]
}
### get bootstrap p-value
hist(t1.star, nclass=40)
abline(v=t1, lwd=3)
abline(v=-1*t1, lwd=3)
2 * mean( t1.star > abs(t1) )
##########################################...
2004 Dec 14
7
deploying dspam
Hi everyone,
I've been thinking of phasing out my spamassassin installation (it's
hardly doing me any good these days) in favour of dspam. Now since I'm
working on this anyway, I thought I'd configure it so when users get any
spam, it is delivered into a special folder. I figure I need the 1.0
test series for this, so I can make that folder fixed. But, I also want
to allow