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 <djnordl...@frontier.com> À : r-help@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@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@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.