Re: [R] RODBC package not found (in English)
Dear Xavier,
Perhaps, you should read the faq:
http://cran.r-project.org/doc/FAQ/R-FAQ.html
And also "un livre en français sur R* : R pour les débutants *par Emmanuel
Paradis ou *Introduction à R par *Julien Barnier"
> install.packages("RODBC")
Un dernier petit mot : Méfiez vous! La patience des membres de la mailing
list R est limitée surtout pour les questions qui ont été posées mille fois
et dont les réponses sont dans le FAQ.
A bon entendeur, Bonne journée
Guillaume
2015-04-20 13:01 GMT+02:00 CHIRIBOGA Xavier :
> Dear members,
>
>
>
> What can I do if I get this message: ?
>
>
>
> library(RODBC)
> Error in library(RODBC) : any package called ‘RODBC’ was found
>
>
>
>
>
> Thanks in advcance,
>
>
>
> Xavier
>
> __
> 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.
>
[[alternative HTML version deleted]]
__
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.
Re: [R] simulate an gradually increase of the number of subjects based on two variables
I omit to precise that I already try to generate data based on the mean and
sd of two variables.
x=rnorm(20,1,5)+1:20
y=rnorm(20,1,7)+41:60
simu<-function(x,y,n) {
simu=vector("list",length=n)
for(i in 1:n) {
x=c(x,rnorm(1,mean(x),sd(x)))
y=c(y,rnorm(1,mean(y),sd(y)))
simu[[i]]$x<-x
simu[[i]]$y<-y
}
return(simu)
}
test=simu(x,y,60)
lapply(test, function(x) cor.test(x$x,x$y))
As you could see, the correlation is disappearing with increasing N.
Perhaps, a bootstrap with lm or cor.test could solve my problem.
2012/3/13 guillaume chaumet
> Dear R list,
> I have a population with two groups. I want to simulate an gradually
> increase of the number of subjects for group 1 based on mean and sd of two
> variables (correlated).
> Bootstrap ?
> Sample ?
> Simulation ? (
>
> I just search some clues.
> Thank you
>
> Guillaume
>
>
>
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
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] simulate an gradually increase of the number of subjects based on two variables
Dear R list, I have a population with two groups. I want to simulate an gradually increase of the number of subjects for group 1 based on mean and sd of two variables (correlated). Bootstrap ? Sample ? Simulation ? ( I just search some clues. Thank you Guillaume [[alternative HTML version deleted]] __ R-help@r-project.org mailing list 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.
Re: [R] as.POSIXct month problem
I am on Ubuntu Karmic.
2009/10/26 Peter Dalgaard
> guillaume chaumet wrote:
>
>> Thank you for your quick response
>>
>>
>> sessionInfo()
>>>
>> R version 2.9.2 (2009-08-24)
>> x86_64-pc-linux-gnu
>>
>> locale:
>>
>> LC_CTYPE=fr_FR.UTF-8;LC_NUMERIC=C;LC_TIME=fr_FR.UTF-8;LC_COLLATE=fr_FR.UTF-8;LC_MONETARY=fr_FR.UTF-8;LC_MESSAGES=fr_FR.UTF-8;LC_PAPER=fr_FR.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=fr_FR.UTF-8;LC_IDENTIFICATION=C
>>
>> attached base packages:
>> [1] stats graphics grDevices utils datasets methods base
>>
>> I need to change fr by en?
>>
>>
> Probably, since your input is in English.
>
> However, it is still a bit puzzling. I have (on Fedora)
>
> > x <- c("1jan1960", "2jan1960", "31mar1960", "30jui1960")
>
> > strptime(x, "%d%b%Y")
> [1] "1960-01-01" "1960-01-02" "1960-03-31" NA
> > Sys.setlocale("LC_TIME","fr_FR.UTF-8")
> [1] "fr_FR.UTF-8"
>
> > strptime(x, "%d%b%Y")
> [1] "1960-01-01" "1960-01-02" "1960-03-31" "1960-07-30"
> > x <- c("1jan1960", "2jan1960", "31mar1960", "30jul1960")
> > strptime(x, "%d%b%Y")
> [1] "1960-01-01" "1960-01-02" "1960-03-31" NA
>
> So there are understandable issues with July (and February, April, May, and
> August), but I don't see any reason to get four NAs. However, you are at the
> mercy of your OS's locale database (you're not saying which OS, just that it
> is a linux variant).
>
> --
> O__ Peter Dalgaard Øster Farimagsgade 5, Entr.B
> c/ /'_ --- Dept. of Biostatistics PO Box 2099, 1014 Cph. K
> (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918
> ~~ - (p.dalga...@biostat.ku.dk) FAX: (+45) 35327907
>
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
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.
Re: [R] as.POSIXct month problem
That's work with
Sys.setlocale("LC_TIME","en_US.UTF-8")
Thank you
2009/10/26 Jorge Ivan Velez
> Hi Guillaume,
>
> What is your sessionInfo() ? It works me:
>
> > x <- c("1jan1960", "2jan1960", "31mar1960", "30jul1960")
> > z <- strptime(x, "%d%b%Y")
> > z
> [1] "1960-01-01" "1960-01-02" "1960-03-31" "1960-07-30"
>
> > sessionInfo()
> R version 2.9.2 RC (2009-08-23 r49375)
> i386-pc-mingw32
>
> locale:
> LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
> States.1252;LC_MONETARY=English_United
> States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252
>
> attached base packages:
> [1] stats graphics grDevices utils datasets methods base
>
> HTH,
> Jorge
>
>
> On Mon, Oct 26, 2009 at 4:31 PM, guillaume chaumet <> wrote:
>
>> Hi everybody
>> When I try example of strptime
>> x <- c("1jan1960", "2jan1960", "31mar1960", "30jul1960")
>> z <- strptime(x, "%d%b%Y")
>>
>> The result is;
>>
>> > z
>> [1] NA NA NA NA
>>
>> I have got the same result with complete form of month but not with
>> numeric
>> form.
>>
>> Any idea?
>>
>>[[alternative HTML version deleted]]
>>
>> __
>> R-help@r-project.org mailing list
>> 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.
>>
>
>
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
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.
Re: [R] as.POSIXct month problem
Thank you for your quick response
>sessionInfo()
R version 2.9.2 (2009-08-24)
x86_64-pc-linux-gnu
locale:
LC_CTYPE=fr_FR.UTF-8;LC_NUMERIC=C;LC_TIME=fr_FR.UTF-8;LC_COLLATE=fr_FR.UTF-8;LC_MONETARY=fr_FR.UTF-8;LC_MESSAGES=fr_FR.UTF-8;LC_PAPER=fr_FR.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=fr_FR.UTF-8;LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
I need to change fr by en?
2009/10/26 Jorge Ivan Velez
> Hi Guillaume,
>
> What is your sessionInfo() ? It works me:
>
> > x <- c("1jan1960", "2jan1960", "31mar1960", "30jul1960")
> > z <- strptime(x, "%d%b%Y")
> > z
> [1] "1960-01-01" "1960-01-02" "1960-03-31" "1960-07-30"
>
> > sessionInfo()
> R version 2.9.2 RC (2009-08-23 r49375)
> i386-pc-mingw32
>
> locale:
> LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
> States.1252;LC_MONETARY=English_United
> States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252
>
> attached base packages:
> [1] stats graphics grDevices utils datasets methods base
>
> HTH,
> Jorge
>
>
> On Mon, Oct 26, 2009 at 4:31 PM, guillaume chaumet <> wrote:
>
>> Hi everybody
>> When I try example of strptime
>> x <- c("1jan1960", "2jan1960", "31mar1960", "30jul1960")
>> z <- strptime(x, "%d%b%Y")
>>
>> The result is;
>>
>> > z
>> [1] NA NA NA NA
>>
>> I have got the same result with complete form of month but not with
>> numeric
>> form.
>>
>> Any idea?
>>
>>[[alternative HTML version deleted]]
>>
>> __
>> R-help@r-project.org mailing list
>> 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.
>>
>
>
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
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] as.POSIXct month problem
Hi everybody
When I try example of strptime
x <- c("1jan1960", "2jan1960", "31mar1960", "30jul1960")
z <- strptime(x, "%d%b%Y")
The result is;
> z
[1] NA NA NA NA
I have got the same result with complete form of month but not with numeric
form.
Any idea?
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
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] building experimental paradigm with R as "Brainard/Pelli Psych Toolbox"
Hi R users, I would known if any attempt of building psychological experimental paradigm with R as "Brainard/Pelli Psych Toolbox" on Matlab or e-prime was done. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list 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] lme equivalent of formula aov
Hi R people, In lme package, I'm searching to find equivalent formula of: >aov(frt~consistency*length*context+Error(subject/(consistency*length*context),data=agE1B4) I try this : >lme(fixed=frt~consistency*length*context, random=~1|subject, data=agE1B4) but I did not obtain same F value with anova(). Thanks in advance Guillaume [[alternative HTML version deleted]] __ R-help@r-project.org mailing list 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] aov within subject and contrasts
Dear R list, again a newbie question here, so I wish I do not exasperated reader. This is example of my data.frame subject length consistency context acc frt 1 1200high scene_c 0.8181818 396.6642 2 2200high scene_c 1.000 595.7778 3 3200high scene_c 0.9090909 510.7315 4 4200high scene_c 0.900 503. 5 5200high scene_c 0.400 523. 6 6200high scene_c 1.000 811.5556 7 7200high scene_c 1.000 661.6402 8 8200high scene_c 1.000 395. 9 9200high scene_c 1.000 514.0909 10 10200high scene_c 0.900 654.6012 11 11200high scene_c 0.900 400. 12 12200high scene_c 0.889 631.1250 ... 491800high scene_c 1.000 376.6667 502800high scene_c 1.000 606.2727 513800high scene_c 0.900 541.7778 524800high scene_c 1.000 483.7273 535800high scene_c 0.556 472.4087 ... My experimental plan contain only "within variables": length (200 or 800 ms) consistency (high or low) and context(scene_c, scene_f, or scene_n) and so I used aov(frt~length*consistency*context + Error(subject/(length*consistency*context))) with one row per observation. How I could obtain simply contrasts for each variables and for interaction? Do I need using gmodels library with lme? Thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list 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.
Re: [R] main title x title and y title with ggplot2
Thierry and Ingo,
Thanks for these smart responses. It works fine.
Guillaume
2008/3/5, Ingo Michaelis <[EMAIL PROTECTED]>:
>
>
>
> ONKELINX, Thierry inbo.be> writes:
>
> >
> > Guillaume,
> >
> > You'll have to add the appropriate scales.
> >
> > ggplot(data, aes(x,y,fill)) + geom_point() + scale_x_continuous("your
> > xlabel") + scale_y_continuous("your ylabel")
> >
> > I suppose you can add a main title in a similar way, but I haven't found
> > that yet. But I shure that Hadley will answer this.
> >
> > Thierry
>
>
> Guillaume,
>
> I had the same problem and found a solution in some forums. Try this:
>
> p<-ggplot(data, aes(x,y,fill)) + geom_point() + scale_x_continuous("your
>
> xlabel") + scale_y_continuous("your ylabel")
>
>
> p$title<-"your title"
>
> print(p)
>
>
> Ingo.
>
>
> __
> R-help@r-project.org mailing list
> 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.
>
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
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.
Re: [R] main title x title and y title with ggplot2
Thierry First thank you for the celerity of your response. Second I use ggplot2 like this : >ggplot(data, aes(x,y,fill)) + geom_point() + etc. Where did you your xlab and ylab when using ggplot2 like that? Guillaume 2008/3/5, ONKELINX, Thierry <[EMAIL PROTECTED]>: > > Guillaume, > > Have a look at the ggplot book on p. 29 > (http://had.co.nz/ggplot2/book.pdf). > > HTH, > > Thierry > > > > ir. Thierry Onkelinx > Instituut voor natuur- en bosonderzoek / Research Institute for Nature > and Forest > Cel biometrie, methodologie en kwaliteitszorg / Section biometrics, > methodology and quality assurance > Gaverstraat 4 > 9500 Geraardsbergen > Belgium > tel. + 32 54/436 185 > [EMAIL PROTECTED] > www.inbo.be > > Do not put your faith in what statistics say until you have carefully > considered what they do not say. ~William W. Watt > A statistical analysis, properly conducted, is a delicate dissection of > uncertainties, a surgery of suppositions. ~M.J.Moroney > > -----Oorspronkelijk bericht- > Van: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] > Namens guillaume chaumet > Verzonden: woensdag 5 maart 2008 8:51 > Aan: r-help@r-project.org > Onderwerp: [R] main title x title and y title with ggplot2 > > > Hi R people, > I'm a R newbie and I'm trying to put main title, x title and y title in > my > graph with no success. > Any idea? > I'm sorry for this newbie question.. > > Guillaume > > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list > 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. > [[alternative HTML version deleted]] __ R-help@r-project.org mailing list 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] main title x title and y title with ggplot2
Hi R people, I'm a R newbie and I'm trying to put main title, x title and y title in my graph with no success. Any idea? I'm sorry for this newbie question.. Guillaume [[alternative HTML version deleted]] __ R-help@r-project.org mailing list 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] anova repeated measure
Dear R people, I'm new user of R and I'm trying to do a Anova with 3 different tests (intra subject), 6 times of measures (intra subject) and 3 groups (inter subject). I want to obtain an Huynh-Feldt epsilon on accuracy with: >anova(lm(acc.t1~gr.dem),test="Spherical") for test 1 # gr.dem is my group variable >anova(lm(acc.t2~gr.dem),test="Spherical") for test 2 >anova(lm(acc.t3~gr.dem),test="Spherical") for test 3 acc.t1, acc.t2 and acc.t3 was data frames with 6 columns each the results was something like this: Analysis of Variance Table Greenhouse-Geisser epsilon: 0.2746 Huynh-Feldt epsilon:0.2915 Df F num Df den Df Pr(>F) G-G Pr H-F Pr (Intercept) 1 7352.7017 6150 1.640e-182 6.894e-52 6.269e-55 gr.dem10.6222 61500.712300.511700.52064 Residuals25 And so what could I do with my epsilon? Previous aov() revealed me that there no effects of my groups but is my 6 times measures factor and my test factor were always significant? Guillaume [[alternative HTML version deleted]] __ R-help@r-project.org mailing list 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] anova (repeated measure)
Hi R people, I'm trying to do a Anova with 3 different tests (intra subject), 6 times of measures (intra subject) and 3 groups (inter subject). I want to obtain an Huynh-Feldt epsilon on accuracy with >anova(lm(acc.t1~gr.dem),test="Spherical") for test 1 # gr.dem is my group variable >anova(lm(acc.t2~gr.dem),test="Spherical") for test 2 >anova(lm(acc.t3~gr.dem),test="Spherical") for test 3 acc.t1, acc.t2 and acc.t3 was data frames with 6 columns each the results was something like this: Analysis of Variance Table Greenhouse-Geisser epsilon: 0.2746 Huynh-Feldt epsilon:0.2915 Df F num Df den Df Pr(>F) G-G Pr H-F Pr (Intercept) 1 7352.7017 6150 1.640e-182 6.894e-52 6.269e-55 gr.dem10.6222 61500.712300.511700.52064 Residuals25 And so what could I do with my epsilon? Previous aov() revealed me that there no effects of my groups but is my 6 times measures factor and my test factor were always significant? Guillaume [[alternative HTML version deleted]] __ R-help@r-project.org mailing list 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.
