[R] boot: how is t0 calculated?

David Winsemius dwinsemius at comcast.net
Tue Oct 23 07:15:31 CEST 2012


On Oct 22, 2012, at 3:28 PM, Thomas Hoffmann wrote:

> Dear all,
> 
> I am doing log-transformed bootstrap regression using:
> 
> x=c(0.038,0.054,1.030,1.250,10.240,52.000,228.100,240.000,758.000,1502.600,0.140,9.170,280.000,1.000,0.150,0.388,20)
> y=c(3961.5,25987.5,526557,321094.5,6433332,23592715.5,40500000,
> 228853521.1,320980392,429000000,58435.5,13525240.5,621195600,1020000,30567.0,296100.0,4800000)
> 
> xy = data.frame(x=x,y=y)
> 
> reg.ln = function(storage, indices){
>  storage = storage[indices,]
>  res.lm  = lm(log(y)~log(x), data=storage)
>  coefficients(res.lm)
> }
> 
> xy.boot = boot(xy, reg.ln, 2000)
> 
> 
> 
> Why does the Intercept given by xy.boot$t0 differs from mean(co.boot$t[,1])?

Hard to tell from this perspective:

> mean(co.boot$t[,1])
Error in mean(co.boot$t[, 1]) : 
  error in evaluating the argument 'x' in selecting a method for function 'mean': Error: object 'co.boot' not found

Are you asking why 12.95 is not equal to 12.94?

> mean(xy.boot$t[,1])
[1] 12.95764
> lm(log(y)~log(x), data=xy) 

Call:
lm(formula = log(y) ~ log(x), data = xy)

Coefficients:
(Intercept)       log(x)  
     12.942        1.055  


> How are the t0 values calculated?

t0?

> Any help is appreciated!

That is what some of us are thinking, too.

-- 

David Winsemius, MD
Alameda, CA, USA




More information about the R-help mailing list