Dear R-Experts, I am trying to compute the BCa nonparametric bootstrap on regression coefficients. Here is the reproducible example : GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12) Quality.score <-c(12,11,13,14,15,16,12,10,9,9) Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7) fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score) confint(fit, level=.95) confint.default(fit, level=.95) confint(fit,level=.95,type="bca") I am not sure but I think that I can not get the nonparametric BCa bootstrap with the confint function. As you can see, I have tried the argument type="bca", I don?t get any error message, but the results don?t change, the results are exactly the same as confint(fit,level=.95). As I have understood, the default argument uses normal quantiles and the method for linear models uses T-quantiles instead. So, I have checked the boot package and the boot.ci function to calculate the BCa bootstrap on the regression coefficients, but I don?t really understand how to compute the code. So, any help from you would be highly appreciated. Best, S
BIG difficulties in Using boot.ci (bot package)
8 messages · varin sacha, Adams, Jean, Michael Dewey +2 more
2 days later
S, There is no mention of a type="bca" argument on the ?confint help file. You can look here for an example of using the boot.ci() function in the boot package: http://www.statmethods.net/advstats/bootstrapping.html ?Jean?
On Fri, Apr 10, 2015 at 11:01 AM, varin sacha <varinsacha at yahoo.fr> wrote:
Dear R-Experts, I am trying to compute the BCa nonparametric bootstrap on regression coefficients. Here is the reproducible example : GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12) Quality.score <-c(12,11,13,14,15,16,12,10,9,9) Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7) fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score) confint(fit, level=.95) confint.default(fit, level=.95) confint(fit,level=.95,type="bca") I am not sure but I think that I can not get the nonparametric BCa bootstrap with the confint function. As you can see, I have tried the argument type="bca", I don?t get any error message, but the results don?t change, the results are exactly the same as confint(fit,level=.95). As I have understood, the default argument uses normal quantiles and the method for linear models uses T-quantiles instead. So, I have checked the boot package and the boot.ci function to calculate the BCa bootstrap on the regression coefficients, but I don?t really understand how to compute the code. So, any help from you would be highly appreciated. Best, S
______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Hi Jean, Many thanks, I got it but there is still a problem. When trying to bootstrap the confidence intervals, I get these messages. boot.ci(results,type="bca",index=1) [1] "All values of t are equal to 5.75620151906917 \n Cannot calculate confidence intervals" NULL
boot.ci(results,type="bca",index=2)
[1] "All values of t are equal to 0.618293471234648 \n Cannot calculate confidence intervals" NULL
boot.ci(results,type="bca",index=3)
[1] "All values of t are equal to 0.148068842921784 \n Cannot calculate confidence intervals"
NULL
Here is the reproducible example :
GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12)
Quality.score <-c(12,11,13,14,15,16,12,10,9,9)
Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7)
fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score)
install.packages("boot")
library(boot)
bs=function(formula,data,indices){
d=data[indices,]
fit=lm(formula,data=Dataset)
return(coef(fit))
}
results=boot(data=Dataset,statistic=bs, R=2000,formula= GDP.LOG ~ Quality.score + Competitivness.score)
boot.ci(results,type="bca",index=1)
boot.ci(results,type="bca",index=2)
boot.ci(results,type="bca",index=3)
How to solve that problem ?
Best,
S.
________________________________
De : "Adams, Jean" <jvadams at usgs.gov>
Cc : "r-help at r-project.org" <r-help at r-project.org>
Envoy? le : Lundi 13 avril 2015 14h22
Objet : Re: [R] BIG difficulties in Using boot.ci (bot package)
S,
There is no mention of a type="bca" argument on the ?confint help file.
You can look here for an example of using the boot.ci() function in the boot package:
http://www.statmethods.net/advstats/bootstrapping.html
?Jean?
Dear R-Experts,
I am trying to compute the BCa nonparametric bootstrap on regression coefficients. Here is the reproducible example : GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12) Quality.score <-c(12,11,13,14,15,16,12,10,9,9) Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7) fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score) confint(fit, level=.95) confint.default(fit, level=.95) confint(fit,level=.95,type="bca") I am not sure but I think that I can not get the nonparametric BCa bootstrap with the confint function. As you can see, I have tried the argument type="bca", I don?t get any error message, but the results don?t change, the results are exactly the same as confint(fit,level=.95). As I have understood, the default argument uses normal quantiles and the method for linear models uses T-quantiles instead. So, I have checked the boot package and the boot.ci function to calculate the BCa bootstrap on the regression coefficients, but I don?t really understand how to compute the code. So, any help from you would be highly appreciated. Best, S
______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
See in line
On 13/04/2015 14:38, varin sacha wrote:
Hi Jean, Many thanks, I got it but there is still a problem. When trying to bootstrap the confidence intervals, I get these messages. boot.ci(results,type="bca",index=1) [1] "All values of t are equal to 5.75620151906917 \n Cannot calculate confidence intervals" NULL
boot.ci(results,type="bca",index=2)
[1] "All values of t are equal to 0.618293471234648 \n Cannot calculate confidence intervals" NULL
boot.ci(results,type="bca",index=3)
[1] "All values of t are equal to 0.148068842921784 \n Cannot calculate confidence intervals" NULL Here is the reproducible example : GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12) Quality.score <-c(12,11,13,14,15,16,12,10,9,9) Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7)
Did you mean at this point to put these into a data frame?
fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score)
install.packages("boot")
library(boot)
bs=function(formula,data,indices){
d=data[indices,]
fit=lm(formula,data=Dataset)
In the code you sent there does not seem to be an object called Dataset. Did you mean data = d?
return(coef(fit))
You do not need to say return( ) at the end, just coef(fit) would work.
} results=boot(data=Dataset,statistic=bs, R=2000,formula= GDP.LOG ~ Quality.score + Competitivness.score) boot.ci(results,type="bca",index=1) boot.ci(results,type="bca",index=2) boot.ci(results,type="bca",index=3) How to solve that problem ? Best, S.
________________________________ De : "Adams, Jean" <jvadams at usgs.gov> Cc : "r-help at r-project.org" <r-help at r-project.org> Envoy? le : Lundi 13 avril 2015 14h22 Objet : Re: [R] BIG difficulties in Using boot.ci (bot package) S, There is no mention of a type="bca" argument on the ?confint help file. You can look here for an example of using the boot.ci() function in the boot package: http://www.statmethods.net/advstats/bootstrapping.html ?Jean? Dear R-Experts, I am trying to compute the BCa nonparametric bootstrap on regression coefficients. Here is the reproducible example : GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12) Quality.score <-c(12,11,13,14,15,16,12,10,9,9) Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7) fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score) confint(fit, level=.95) confint.default(fit, level=.95) confint(fit,level=.95,type="bca") I am not sure but I think that I can not get the nonparametric BCa bootstrap with the confint function. As you can see, I have tried the argument type="bca", I don?t get any error message, but the results don?t change, the results are exactly the same as confint(fit,level=.95). As I have understood, the default argument uses normal quantiles and the method for linear models uses T-quantiles instead. So, I have checked the boot package and the boot.ci function to calculate the BCa bootstrap on the regression coefficients, but I don?t really understand how to compute the code. So, any help from you would be highly appreciated. Best, S ______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Hi Michael,
Thanks for your response. About the data frame not necessary. I correct the code according to your comments. I still get the following warnings :
[1] "All values of t are equal to 5.75620151906917 \n Cannot calculate confidence intervals"
NULL
I have found this on the Net :
"Note that boot.ci just gives a warning and returns NA values, if all values are equal. There is no error and if you can work with NAs, there is no need for the if condition.
The boot package assumes that the bootstrap
statistic has all ways the same dimension.
Whenever you have a statistic with less dimensions you get an NA
or 0 or whatever you want".
The reproducible code :
GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12)
Quality.score <-c(12,11,13,14,15,16,12,10,9,9)
Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7)
fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score)
install.packages("boot")
library(boot)
bs=function(formula,data,indices){
d=data[indices,]
fit=lm(formula,data=d)
(coef(fit))
}
results=boot(data=Dataset,statistic=bs, R=2000,formula= GDP.LOG ~ Quality.score + Competitivness.score)
boot.ci(results,type="bca",index=1)
boot.ci(results,type="bca",index=2)
boot.ci(results,type="bca",index=3)
How is it possible to avoid that warning ?
Best,
S
----- Mail original -----
De : Michael Dewey <lists at dewey.myzen.co.uk>
? : varin sacha <varinsacha at yahoo.fr>; "Adams, Jean" <jvadams at usgs.gov>
Cc : "r-help at r-project.org" <r-help at r-project.org>
Envoy? le : Lundi 13 avril 2015 16h00
Objet : Re: [R] BIG difficulties in Using boot.ci (bot package)
See in line
On 13/04/2015 14:38, varin sacha wrote:
Hi Jean, Many thanks, I got it but there is still a problem. When trying to bootstrap the confidence intervals, I get these messages. boot.ci(results,type="bca",index=1) [1] "All values of t are equal to 5.75620151906917 \n Cannot calculate confidence intervals" NULL
boot.ci(results,type="bca",index=2)
[1] "All values of t are equal to 0.618293471234648 \n Cannot calculate confidence intervals" NULL
boot.ci(results,type="bca",index=3)
[1] "All values of t are equal to 0.148068842921784 \n Cannot calculate confidence intervals" NULL Here is the reproducible example : GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12) Quality.score <-c(12,11,13,14,15,16,12,10,9,9) Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7)
Did you mean at this point to put these into a data frame?
fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score)
install.packages("boot")
library(boot)
bs=function(formula,data,indices){
d=data[indices,]
fit=lm(formula,data=Dataset)
In the code you sent there does not seem to be an object called Dataset. Did you mean data = d?
return(coef(fit))
You do not need to say return( ) at the end, just coef(fit) would work.
} results=boot(data=Dataset,statistic=bs, R=2000,formula= GDP.LOG ~ Quality.score + Competitivness.score) boot.ci(results,type="bca",index=1) boot.ci(results,type="bca",index=2) boot.ci(results,type="bca",index=3) How to solve that problem ? Best, S.
________________________________ De : "Adams, Jean" <jvadams at usgs.gov> Cc : "r-help at r-project.org" <r-help at r-project.org> Envoy? le : Lundi 13 avril 2015 14h22 Objet : Re: [R] BIG difficulties in Using boot.ci (bot package) S, There is no mention of a type="bca" argument on the ?confint help file. You can look here for an example of using the boot.ci() function in the boot package: http://www.statmethods.net/advstats/bootstrapping.html ?Jean? Dear R-Experts, I am trying to compute the BCa nonparametric bootstrap on regression coefficients. Here is the reproducible example : GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12) Quality.score <-c(12,11,13,14,15,16,12,10,9,9) Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7) fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score) confint(fit, level=.95) confint.default(fit, level=.95) confint(fit,level=.95,type="bca") I am not sure but I think that I can not get the nonparametric BCa bootstrap with the confint function. As you can see, I have tried the argument type="bca", I don?t get any error message, but the results don?t change, the results are exactly the same as confint(fit,level=.95). As I have understood, the default argument uses normal quantiles and the method for linear models uses T-quantiles instead. So, I have checked the boot package and the boot.ci function to calculate the BCa bootstrap on the regression coefficients, but I don?t really understand how to compute the code. So, any help from you would be highly appreciated. Best, S ______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
On 4/13/2015 9:06 AM, varin sacha wrote:
Hi Michael,
Thanks for your response. About the data frame not necessary. I correct the code according to your comments. I still get the following warnings :
[1] "All values of t are equal to 5.75620151906917 \n Cannot calculate confidence intervals"
NULL
I have found this on the Net :
"Note that boot.ci just gives a warning and returns NA values, if all values are equal. There is no error and if you can work with NAs, there is no need for the if condition.
The boot package assumes that the bootstrap
statistic has all ways the same dimension.
Whenever you have a statistic with less dimensions you get an NA
or 0 or whatever you want".
The reproducible code :
GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12)
Quality.score <-c(12,11,13,14,15,16,12,10,9,9)
Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7)
fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score)
install.packages("boot")
library(boot)
bs=function(formula,data,indices){
d=data[indices,]
fit=lm(formula,data=d)
(coef(fit))
}
results=boot(data=Dataset,statistic=bs, R=2000,formula= GDP.LOG ~ Quality.score + Competitivness.score)
boot.ci(results,type="bca",index=1)
boot.ci(results,type="bca",index=2)
boot.ci(results,type="bca",index=3)
How is it possible to avoid that warning ?
Best,
S
The first problem is that this is not reproducible code. This is what I get when I run your code on my computer:
GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12)
Quality.score <-c(12,11,13,14,15,16,12,10,9,9)
Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7)
fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score)
#####install.packages("boot")
library(boot)
bs=function(formula,data,indices){
+ d=data[indices,] + fit=lm(formula,data=d) + (coef(fit)) + }
results=boot(data=Dataset,statistic=bs, R=2000,formula= GDP.LOG ~ Quality.score + Competitivness.score)
Error in NROW(data) : object 'Dataset' not found
boot.ci(results,type="bca",index=1)
Error in boot.ci(results, type = "bca", index = 1) : object 'results' not found
boot.ci(results,type="bca",index=2)
Error in boot.ci(results, type = "bca", index = 2) : object 'results' not found
boot.ci(results,type="bca",index=3)
Error in boot.ci(results, type = "bca", index = 3) : object 'results' not found
A reproducible example means that when I run your code on my machine, I get the same results / warnings / errors that you get. I got something different. Dan
Daniel Nordlund Bothell, WA USA
Hi Daniel,
Sorry for that, once more.... ;=(
Here is the reproducible code and this time IT WORKS FINALLY !!!
GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12)
Quality.score <-c(12,11,13,14,15,16,12,10,9,9)
Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7)
df=data.frame(GDP.LOG,Quality.score,Competitivness.score)
fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score)
install.packages("boot")
library(boot)
bs=function(formula,data,indices){
d=data[indices,]
fit=lm(formula,data=d)
(coef(fit))
}
results=boot(data=df,statistic=bs, R=2000,formula= GDP.LOG ~
Quality.score + Competitivness.score)
boot.ci(results,type="bca",index=1)
boot.ci(results,type="bca",index=2)
boot.ci(results,type="bca",index=3)
Best,
S
----- Mail original -----
De : Daniel Nordlund <djnordlund at frontier.com>
? : r-help at r-project.org
Cc :
Envoy? le : Lundi 13 avril 2015 19h26
Objet : Re: [R] BIG difficulties in Using boot.ci (bot package)
On 4/13/2015 9:06 AM, varin sacha wrote:
Hi Michael,
Thanks for your response. About the data frame not necessary. I correct the code according to your comments. I still get the following warnings :
[1] "All values of t are equal to 5.75620151906917 \n Cannot calculate confidence intervals"
NULL
I have found this on the Net :
"Note that boot.ci just gives a warning and returns NA values, if all values are equal. There is no error and if you can work with NAs, there is no need for the if condition.
The boot package assumes that the bootstrap
statistic has all ways the same dimension.
Whenever you have a statistic with less dimensions you get an NA
or 0 or whatever you want".
The reproducible code :
GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12)
Quality.score <-c(12,11,13,14,15,16,12,10,9,9)
Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7)
fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score)
install.packages("boot")
library(boot)
bs=function(formula,data,indices){
d=data[indices,]
fit=lm(formula,data=d)
(coef(fit))
}
results=boot(data=Dataset,statistic=bs, R=2000,formula= GDP.LOG ~ Quality.score + Competitivness.score)
boot.ci(results,type="bca",index=1)
boot.ci(results,type="bca",index=2)
boot.ci(results,type="bca",index=3)
How is it possible to avoid that warning ?
Best,
S
The first problem is that this is not reproducible code. This is what I get when I run your code on my computer:
GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12)
Quality.score <-c(12,11,13,14,15,16,12,10,9,9)
Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7)
fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score)
#####install.packages("boot")
library(boot)
bs=function(formula,data,indices){
+ d=data[indices,] + fit=lm(formula,data=d) + (coef(fit)) + }
results=boot(data=Dataset,statistic=bs, R=2000,formula= GDP.LOG ~ Quality.score + Competitivness.score)
Error in NROW(data) : object 'Dataset' not found
boot.ci(results,type="bca",index=1)
Error in boot.ci(results, type = "bca", index = 1) : object 'results' not found
boot.ci(results,type="bca",index=2)
Error in boot.ci(results, type = "bca", index = 2) : object 'results' not found
boot.ci(results,type="bca",index=3)
Error in boot.ci(results, type = "bca", index = 3) : object 'results' not found
A reproducible example means that when I run your code on my machine, I get the same results / warnings / errors that you get. I got something different. Dan
Daniel Nordlund Bothell, WA USA ______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
That is probably the number one reason for requesting a reproducible example when writing to R-help. In the proce3ss of working that out, you often solve your own problem. Best of luck with your bootstrapping, Dan Daniel J. Nordlund, PhD Research and Data Analysis Division Services & Enterprise Support Administration Washington State Department of Social and Health Services
-----Original Message-----
From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of varin
sacha
Sent: Monday, April 13, 2015 11:55 AM
To: Daniel Nordlund; r-help at r-project.org
Subject: Re: [R] BIG difficulties in Using boot.ci (bot package)
Hi Daniel,
Sorry for that, once more.... ;=(
Here is the reproducible code and this time IT WORKS FINALLY !!!
GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12)
Quality.score <-c(12,11,13,14,15,16,12,10,9,9)
Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7)
df=data.frame(GDP.LOG,Quality.score,Competitivness.score)
fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score)
install.packages("boot")
library(boot)
bs=function(formula,data,indices){
d=data[indices,]
fit=lm(formula,data=d)
(coef(fit))
}
results=boot(data=df,statistic=bs, R=2000,formula= GDP.LOG ~
Quality.score + Competitivness.score)
boot.ci(results,type="bca",index=1)
boot.ci(results,type="bca",index=2)
boot.ci(results,type="bca",index=3)
Best,
S
----- Mail original -----
De : Daniel Nordlund <djnordlund at frontier.com>
? : r-help at r-project.org
Cc :
Envoy? le : Lundi 13 avril 2015 19h26
Objet : Re: [R] BIG difficulties in Using boot.ci (bot package)
On 4/13/2015 9:06 AM, varin sacha wrote:
Hi Michael, Thanks for your response. About the data frame not necessary. I
correct the code according to your comments. I still get the following warnings :
[1] "All values of t are equal to 5.75620151906917 \n Cannot
calculate confidence intervals"
NULL I have found this on the Net : "Note that boot.ci just gives a warning and returns NA values, if all
values are equal. There is no error and if you can work with NAs, there is no need for the if condition.
The boot package assumes that the bootstrap
statistic has all ways the same dimension.
Whenever you have a statistic with less dimensions you get an NA
or 0 or whatever you want".
The reproducible code :
GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12)
Quality.score <-c(12,11,13,14,15,16,12,10,9,9)
Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7)
fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score)
install.packages("boot")
library(boot)
bs=function(formula,data,indices){
d=data[indices,]
fit=lm(formula,data=d)
(coef(fit))
}
results=boot(data=Dataset,statistic=bs, R=2000,formula= GDP.LOG ~
Quality.score + Competitivness.score)
boot.ci(results,type="bca",index=1) boot.ci(results,type="bca",index=2) boot.ci(results,type="bca",index=3) How is it possible to avoid that warning ? Best, S
The first problem is that this is not reproducible code. This is what I get when I run your code on my computer:
GDP.LOG <-c(14,12,13,15.5,16,17,16.5,13.5,12.5,12)
Quality.score <-c(12,11,13,14,15,16,12,10,9,9)
Competitivness.score=c(8,6,7,5,6.5,7,8,4.5,6,7)
fit <- lm(formula = GDP.LOG ~ Quality.score + Competitivness.score)
#####install.packages("boot")
library(boot)
bs=function(formula,data,indices){
+ d=data[indices,] + fit=lm(formula,data=d) + (coef(fit)) + }
results=boot(data=Dataset,statistic=bs, R=2000,formula= GDP.LOG ~
Quality.score + Competitivness.score) Error in NROW(data) : object 'Dataset' not found
boot.ci(results,type="bca",index=1)
Error in boot.ci(results, type = "bca", index = 1) : object 'results' not found
boot.ci(results,type="bca",index=2)
Error in boot.ci(results, type = "bca", index = 2) : object 'results' not found
boot.ci(results,type="bca",index=3)
Error in boot.ci(results, type = "bca", index = 3) : object 'results' not found
A reproducible example means that when I run your code on my machine, I get the same results / warnings / errors that you get. I got something different. Dan -- Daniel Nordlund Bothell, WA USA
______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting- guide.html and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting- guide.html and provide commented, minimal, self-contained, reproducible code.