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 <li...@dewey.myzen.co.uk>
À : varin sacha <varinsa...@yahoo.fr>; "Adams, Jean" <jvad...@usgs.gov>
Cc : "r-help@r-project.org" <r-help@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" <jvad...@usgs.gov>
>
> Cc : "r-help@r-project.org" <r-help@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@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.
>
--
Michael
http://www.dewey.myzen.co.uk/home.html
______________________________________________
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.
