search for: varinsacha

...posting guide > http://www.r-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] ------------------------------ Message: 3 Date: Tue, 12 Dec 2023 21:19:12 +0000 (UTC) From: varin sacha <varinsacha at yahoo.fr> To: "r-help at r-project.org" <r-help at r-project.org>, Ben Bolker <bbolker at gmail.com> Subject: Re: [R] ggplot2: Get the regression line with 95% confidence bands Message-ID: <68588390.888662.1702415952477 at mail.yahoo.com> Content-T...

Sorry, didn't cc this to the list. -------- Forwarded Message -------- Subject: Re: [R] Strange results : bootrstrp CIs Date: Sat, 13 Jan 2024 17:37:19 -0500 From: Duncan Murdoch <murdoch.duncan at gmail.com> To: varin sacha <varinsacha at yahoo.fr> You can debug things like this by setting options(error = recover). That will drop into the debugger when the error occurs. Examine t.star, r, and res[[r]] and you will likely see what the problem was. My guess is that one of the bootstrap samples had a different selection of...

Dear R-Experts, Thanks to Prof. Bonnett, I have got an R script working to calculate confidence intervals around the semipartial correlation coefficients. Now, I would like to calculate BCa bootstrap CIs using the boot library and the boot.ci(results, type="all") function. How could I modify my R script (here below reproducible example) to get the BCa bootstrap CIs ? CIsemipartcorr

Dear R-experts, Here below is my code, I would like to add a vertical line on my plot, showing the median and I would like to place some values on this graph as well, i.e. 4.3 and -8.4. How can I do ? Many thanks for your reply. A=c(1,2.3,4,3.5,4.3,2.5,6.3,-0.1,-1.5,3.7,-2.3,-3.5,5.4,3.2, -10.5,-8.4,-9.4) d <- density(A) plot(d) median(A)

Hi Deepak, In econometrics there is another test very often used : the white test. The white test is based on the comparison of the estimated variances of residuals when the model is estimated by OLS under the assumption of homoscedasticity and when the model is estimated by OLS under the assumption of heteroscedastic. The White test with R install.packages("bstats") library(bstats)

Dear R-experts, Here below my reproducible R code. I want to add many straight lines to a plot using "abline" The last fit (fast Tau-estimator, color yellow) will not appear on the plot. What is going wrong ? Many thanks for your reply. ########## Y=c(2,4,5,4,3,4,2,3,56,5,4,3,4,5,6,5,4,5,34,21,12,13,12,8,9,7,43,12,19,21)

Dear R-Experts, Here below is my R code (reproducible example) to calculate the confidence intervals around the spearman coefficient. ########## C=c(2,4,5,6,3,4,5,7,8,7,6,5,6,7,7,8,5,4,3,2) D=c(3,5,4,6,7,2,3,1,2,4,5,4,6,4,5,4,3,2,8,9) cor(C,D,method= "spearman") library(boot) myCor=function(data,index){ cor(data[index, ])[1,2] } results=boot(data=cbind(C,D),statistic=myCor, R=2000)

Dear R-experts, I can not install the package "matie". If somebody can tell me where the problem is. > install.packages("matie") Installing package into ?/Users/Caro/Library/R/3.3/library? (as ?lib? is unspecified) essai de l'URL 'https://cran.rstudio.com/bin/macosx/mavericks/contrib/3.3/matie_1.2.tgz' Content type 'application/x-gzip' length 80151

Hi Rui, Many thanks. The R code works BUT the results I get are quite weird I guess ! MIC = 0.2650 Normal 95% CI = (0.9614, 1.0398) The MIC is not inside the confidence intervals ! Is there something wrong in the R code ? Here is the reproducible example : ########## C=c(2,4,5,6,3,4,5,7,8,7,6,5,6,7,7,8,5,4,3,2) D=c(3,5,4,6,7,2,3,1,2,4,5,4,6,4,5,4,3,2,8,9) library(minerva) mine(C,D)$MIC