search for: varin

...ot;,?"Huber M-estimator", "Tukey", "MM", "fast tau"),inset=0.02,lwd=2,col=c("black","blue","red","purple","green", "yellow"),cex=.9)} ############### Le samedi 31 mars 2018 ? 21:52:55 UTC+2, varin sacha via R-help <r-help at r-project.org> a ?crit : Many thanks Duncun, Best, Le samedi 31 mars 2018 ? 18:05:53 UTC+2, Duncan Murdoch <murdoch.duncan at gmail.com> a ?crit : On 31/03/2018 11:57 AM, varin sacha via R-help wrote: > Dear R-experts, > > Here bel...

...1469 [75,] 5 1 0.742553847 attr(,"Ordered") [1] FALSE attr(,"call") silhouette.default(x = cutree(ar, k = 5), dist = daisy(ruspini)) attr(,"class") [1] "silhouette" Thanks for your attention, Cristiano --------------------------------- Cristiano Varin sammy@unive.it http://www.dst.unive.it/~sammy/ [[alternative HTML version deleted]]

Dear R-experts, Here below my R code. The warning message is not a problem to me but there is an error message more problematic. I understand the error message but I don't know if it is possible to correct the error and if yes, how to correct it. Many thanks. n <- 60 b <- runif(n, 0, 5) a <- runif(n, 0, 5) z <- rnorm(n*0.95,2,3) + rnorm(n*0.05,2,9) y_model <- 0.1 * b - 0.5 *

...MM" , "fast tau" ) , inset = 0.02 , lwd = 2 , col = c( "black","blue","red","purple","green", "yellow") , cex = .9 ) ###################### On Fri, 6 Apr 2018, varin sacha via R-help wrote: > R-experts, > > I have fitted many different lines. The fast-tau estimator (yellow line) seems strange to me?because this yellow line is not at all in agreement with the other lines (reverse slope, I mean the yellow line has a positive slope and the other ones hav...

Varin In the teaching you to fish department, you can find those dependencies Duncan appeared to pull out of thin air by looking at the CRAN contributed packages web pages, which can easily be found with Google. Start at the matie page and follow the dependency links and look at dependencies. In particu...

Hi, a is of length 60. b is of length 60. z is of length 57. What do you expect to have happen when you create y_model ? What happens to those other 3 observations? Sarah On Tue, Oct 27, 2020 at 3:07 PM varin sacha via R-help <r-help at r-project.org> wrote: > > Dear R-experts, > > Here below my R code. The warning message is not a problem to me but there is an error message more problematic. I understand the error message but I don't know if it is possible to correct the error and...

...unction myCor and all went well, with a warning relative to BCa intervals. myCor <- function(data, index){ mine(data[index, ])$MIC } results=boot(data = cbind(C,D), statistic = myCor, R = 2000) boot.ci(results,type="all") Hope this helps, Rui Barradas On 12/10/2017 3:19 PM, varin sacha via R-help wrote: > 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) &gt...

On 31/03/2018 11:57 AM, varin sacha via R-help wrote: > 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. &g...

...EASE do read the 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.ya...

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 dif...

...ences between: mine(C, D) and mine(cbind(C, D)) The first returns a value, the second returns a symmetric matrix. Just like cor() David L. Carlson Department of Anthropology Texas A&M University -----Original Message----- From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of varin sacha via R-help Sent: Sunday, December 10, 2017 11:07 AM To: Rui Barradas <ruipbarradas at sapo.pt>; R-help Mailing List <r-help at r-project.org> Subject: Re: [R] Confidence intervals around the MIC (Maximal information coefficient) Hi Rui, Many thanks. The R code works BUT the resu...

...mber, or 0 to exit 1: boot(e, func, R = 1000) <enter 1> Selection: 1 Called from: top level Browse[1]> print(r) [1] 2 Browse[1]> t.star[r,] [1] NA NA NA NA NA NA NA NA NA i[2,] [1] 14 15 22 22 21 14 8 2 12 22 10 15 9 7 9 13 12 23 1 20 15 7 5 10 On 2024-01-13 5:22 p.m., varin sacha via R-help wrote: > Dear Duncan, > Dear Ivan, > > I really thank you a lot for your response. > So, if I correctly understand your answers the problem is coming from this line: > > coef(lm(Score~ Time + factor(Country)),data=data[idx,]) > > This line should be: &...

...=data[idx,])) } B= boot(e, func, R=1000) ? boot.ci(B, index=2, type="perc") ############################################################# Le samedi 13 janvier 2024 ? 21:56:58 UTC+1, Ivan Krylov <ikrylov at disroot.org> a ?crit : ? Sat, 13 Jan 2024 20:33:47 +0000 (UTC) varin sacha via R-help <r-help at r-project.org> ?????: > coef(lm(Score~ Time + factor(Country)),data=data[idx,]) Wrong place for the data=... argument. You meant to give it to lm(...), but in the end it went to coef(...). Without the data=... argument, the formula passed to lm() picks up the...

On 28/03/2018 6:27 PM, varin sacha via R-help wrote: > 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) &...

Dear R-experts, Here below the reproducible example. I am trying to get the average of the 100 results coming from the "lst" function. I have tried lst$mean and mean(lst). It does not work. Any help would be highly appreciated. #################### ?## R script for getting MedAe and MedAeSQ from HBR model on Testing data install.packages("robustbase") install.packages(

Dear R-experts, Here below, my R code working BUT I get a strange result I was not expecting! Indeed, the?95% percentile bootstrap CIs is (-54.81, -54.81 ). Is anything going wrong? Best, ########################################## Score=c(345,564,467,675,432,346,476,512,567,543,234,435,654,411,356,658,432,345,432,345, 345,456,543,501) ? Country=c("Italy", "Italy",

? Sat, 13 Jan 2024 20:33:47 +0000 (UTC) varin sacha via R-help <r-help at r-project.org> ?????: > coef(lm(Score~ Time + factor(Country)),data=data[idx,]) Wrong place for the data=... argument. You meant to give it to lm(...), but in the end it went to coef(...). Without the data=... argument, the formula passed to lm() picks up the...

...sing R version 1.1.1 on Windows NT. (I'm sorry if this is a stupid question but I didn't find the answer in the help file or the FAQs.) Thanks in advance for any help. Martina Erdbr?gge Here is my code with results: > lm44 <- lm(betamat2[, 2] ~ 1, data = fndaten, weights = fndaten$varin) > step(lm44, ~ (A + B + C + D)^4) Start: AIC= -20.5 betamat2[, 2] ~ 1 Df Sum of Sq RSS AIC + C 1 1.3500 2.5706 -25.2555 + B 1 1.2323 2.6882 -24.5393 <none> 3.9206 -20.5016 + A 1 0.0132 3.9074 -18.5554 + D 1 0.0072 3....

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)

Hello, After the plot just do abline(v = median(A)) As for how to plot points, see, well, ?points(). Hope this helps, Rui Barradas Em 22-10-2017 16:33, varin sacha via R-help escreveu: > 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. > > > &gt...