Hi R-users,
I've found the error
In these rows:
# Model (a)
testtemp <- lm.bp(doctorco~sex+age+income, prescrib~sex+age+income,
data=bootdata)
betafound <- c(testtemp$beta,testtemp$beta3)
results[i,] <- betafound
betafound must to be equal to:
betafound <- c(testtemp$beta1,testtemp$beta
Thanks for the answer.
The R code is here:
library("bivpois")
data("ex3.health")
# Bivariate Poisson models
ex3.model.a <- lm.bp(doctorco~sex+age+income, prescrib~sex+age+income,
data=ex3.health)
Bootstrap standard errors can be obtained easily using the following script
for model (a):
n <- lengt
Can anyone help me please?
--
View this message in context:
http://r.789695.n4.nabble.com/Bootstraps-standard-error-tp3313322p3315463.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.
On Feb 20, 2011, at 2:20 PM, danielepippo wrote:
when I calculate the bootstrap standard errors
... of what?
the results are summarized in
a matrix with columns equal to the parameters of the model.
... but they are all the same within each row???
In my case I
have a matrix like this:
when I calculate the bootstrap standard errors the results are summarized in
a matrix with columns equal to the parameters of the model. In my case I
have a matrix like this:
[,1] [,2] [,3] [,4] [,5]
[,6][,7] [,8][,9]
[
On Feb 18, 2011, at 1:36 PM, danielepippo wrote:
Hi everyone,
if I have the bootsrap values of the parameters of my regression
model,
how can I calculate the standard error of them to refer to the
parameters?
The same way you would any other vector?
Thank you
--
View this message i
Hi everyone,
if I have the bootsrap values of the parameters of my regression model,
how can I calculate the standard error of them to refer to the parameters?
Thank you
--
View this message in context:
http://r.789695.n4.nabble.com/Bootstraps-standard-error-tp3313322p3313322.html
Sent from
7 matches
Mail list logo
