search for: sacha

Hello, Can you say me which is the key to add in the registry for Windows XP to Enable Plain Text Password? Best Regards

Dear Sacha, use glm() in this case. I'd rather code the covariable as TRUE / FALSE or as a factor. Best regards, ir. Thierry Onkelinx Statisticus / Statistician Vlaamse Overheid / Government of Flanders INSTITUUT VOOR NATUUR- EN BOSONDERZOEK / RESEARCH INSTITUTE FOR NATURE AND FOREST Team Biometrie &a...

Hello! I am doing an analysis on a questionnaire of hunters taken in 4 different districts of some mysterious foreign country. The aim of the study was to gather info on the factors that determine the hunting success of a peculiarly beautiful bird in that area. All variables are factors, i.e. they are variables such as "Use of Guns - yes / no", "Use of Dogs - yes / no" and

...tors (variables coded 0,1). I guess I cannot just smooth binary variables. By the way I code them as 0=no,1=yes, then mgcv will think those variables are numeric.? I have tried to change 0 and 1 in no and yes. It does not work. How to solve my problem. Here below my toy example. Many thanks. Best, Sacha ? ######################## y=c(0.02,0.1,0.04,0.07,0.09,0.1,0.054,0.076,0.0897,0.09,0.1,0.86,0.74,0.53,0.35,0.67,0.63,0.95,0.72,0.85) x1=c(1,0,0,1,1,0,0,1,0,1,1,1,0,1,0,0,0,1,0,0) x2=c(0,0,1,1,0,1,1,0,0,1,1,1,0,0,1,0,0,1,1,0) ? library(mgcv) reg = gam(y ~ s(x1) +s(x2), family=binomial(link="log...

Hi R-users I would like to extract the random effects ("1|SITE", "1+SPECIES|SITE" and "BA|SITE") from this model formula: Full_model <- formula (VAR ~ (1|SITE) + (1+SPECIES|SITE) + (BA|SITE) + HEIGHT + COND + NN_DIST) I tried: terms(Full_model) labels(terms(Full_model)) but I could not distinguish between random and fixed effects. thanks

...dexing. I need to create new mailboxes on the fly and perform fast searches across all the mailboxes. I intend to store all the messages in an SQL database. Any help getting started interfacing with dovecot (other than reading through _all_ the source code) would be greatly appreciated. Thanks, Sacha

...quot;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 below my...

hello list! I'm sorry, I just stumbled over this strange behaviour (at least I am not able to explain the behaviour, therefore I assume it to be a strange behaviour): attach(water) # I know, this is not recommended names(water[3:10]) [1] "temp" "pH" "DO" "BOD" "COD" "no3" "no2" "po4" for (i in

On Mon, 4 Apr 2016, varin sacha via R-help wrote: > 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 wh...

...t; , "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 have nega...

...n 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) > cor(...

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

...o 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.yahoo.co...

Hello, Sacha, le ven. 08 févr. 2019 18:00:22 +0100, a ecrit: > On  Debian GNU/Linux 9.7 (stretch) amd64, we have a bug on the last Xen > Hypervisor version: > > xen-hypervisor-4.8-amd64 4.8.5+shim4.10.2+xsa282 (Read: 4.8.5+shim4.10.2+xsa282-1+deb9u11) > The rollback on the previous package...

...riances 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) white.test(LinearModel) Hope this helps. Sacha ________________________________ De : Deepak Singh <sdeepakrhelp at gmail.com> ? : r-help at r-project.org Envoy? le : Lundi 4 avril 2016 10h40 Objet : [R] Test for Homoscedesticity in R Without BP Test Respected Sir, I am doing a project on multiple linear model fitting and in that p...

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

...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 results I...

...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: > co...

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

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 *