Re: [R] simple main effect.

Try going through this: http://www.personality-project.org/r/r.anova.html http://www.personality-project.org/r/r.anova.html Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com

Re: [R] read.table: mysterious line omissions

Jonathan wrote: Tried deleting several line above and below line 362; did that and changed the extension. Still finds the same number of rows (362). Strange! Jonathan The reason is quite apparent if you try page(temp) Ovine lentivirus, Peanut chlorotic streak virus, Pelargonium vein

[R] utf8 postscript cyrillic

Dear R users, I am running R version 2.10.0 (2009-10-26). I need to prepare an eps graphic with a legend with cyrillic words. I tried setting the encoding parameter of the postscript command, but in vain, nothing seems to work. I tried with CP1251, KOI8-R, UTF-8 and Cyrillic (UTF-8 turned out to

Re: [R] utf8 postscript cyrillic

'Nothing seems to work' is not at all helpful, and you seem under the misconception that PostScript supports UTF-8 (sic): life would be much easier for device writers if it did. The postscript() help page says More details of font families and encodings and especially handling text

[R] how to create three new variables? Thanks a lot!

Dear R experts: I have following data structure: student id, exam time and score. I want to create three new columes: 1 st before, 2nd before and 3rd before. For example, for student 1's 4th exam score is assumed to be influenced by his previous three cloest scores, which are 6, 9, and 10 in

Re: [R] simple main effect.

Thanks. But there is no simple main effect there. On Sun, Dec 20, 2009 at 10:00 AM, Tal Galili tal.gal...@gmail.com wrote: Try going through this: http://www.personality-project.org/r/r.anova.html http://www.personality-project.org/r/r.anova.html Contact

Re: [R] how to count the total number of (INCLUDING overlapping) occurrences of a substring within a string?

Use a zero lookaround expression. It will not consume its match. See ?regexp gregexpr(a(?=a), aaa, perl = TRUE) [[1]] [1] 1 2 attr(,match.length) [1] 1 1 On Sun, Dec 20, 2009 at 1:43 AM, Jonathan jonsle...@gmail.com wrote: Last one for you guys: The command: length(gregexpr('cus','hocus

Re: [R] how to count the total number of (INCLUDING overlapping) occurrences of a substring within a string?

On Sun, Dec 20, 2009 at 5:33 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: Use a zero lookaround expression.  It will not consume its match.  See ?regexp That should be lookahead, not lookaround. gregexpr(a(?=a), aaa, perl = TRUE) [[1]] [1] 1 2 attr(,match.length) [1] 1 1 On

[R] Res: how to create three new variables? Thanks a lot!

Hello! Here is a way: mydata - data.frame(third_before=c(6,10), second_before=c(9,10), first_before=c(10,8)) mydata   third_before second_before first_before 1    6 9   10 2   10    10    8 And if you already have a dataset and wants only

Re: [R] how to count the total number of (INCLUDING overl apping) occurrences of a substring within a string ?

Gabor Grothendieck ggrothendieck at gmail.com writes: Use a zero lookaround expression. It will not consume its match. See ?regexp gregexpr(a(?=a), aaa, perl = TRUE) [[1]] [1] 1 2 attr(,match.length) [1] 1 1 I wonder how you would count the number of occurrences of, for example,

Re: [R] how to count the total number of (INCLUDING overlapping) occurrences of a substring within a string?

Try this: findall(aba, ababacababab) [1] 1 3 7 9 gregexpr(a(?=ba), ababacababab, perl = TRUE) [[1]] [1] 1 3 7 9 attr(,match.length) [1] 1 1 1 1 findall(a.a, ababacababab) [1] 1 3 5 7 9 gregexpr(a(?=.a), ababacababab, perl = TRUE) [[1]] [1] 1 3 5 7 9 attr(,match.length) [1] 1 1 1 1 1 On

Re: [R] read.table: mysterious line omissions

Most likely an unbalanced quote. put the following option in the read.table: quote='', comment.char='' On Sat, Dec 19, 2009 at 11:42 PM, Jonathan jonsle...@gmail.com wrote: Hello again, I am simply trying to import a rectangular table of strings. The table's dimensions are 1990 x 2,

Re: [R] simple main effect.

Hi Or, Maybe I didn't understand you. Are you asking how can I read the output of a fitted (complex within and between) model for finding the simple main effect ? Tal Contact Details:--- Contact me: tal.gal...@gmail.com |

Re: [R] simple main effect.

I don't think so. I'm asking how can I see/analyse the simple main effect. I don't think it shows in the summary report of the aov() function. On Sun, Dec 20, 2009 at 3:35 PM, Tal Galili tal.gal...@gmail.com wrote: Hi Or, Maybe I didn't understand you. Are you asking how can I read the

Re: [R] how to create three new variables? Thanks a lot!

On Dec 20, 2009, at 4:34 AM, jie feng wrote: Dear R experts: I have following data structure: student id, exam time and score. I want to create three new columes: 1 st before, 2nd before and 3rd before. For example, for student 1's 4th exam score is assumed to be influenced by his

Re: [R] simple main effect.

Could you please write the aov formula you are using ? Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com/

Re: [R] simple main effect.

No problem. If I have a mixed model with A as within subject (A is exposure, has 2 levels) and drug(2 levels between subject) and strain(2 levels between subject). I reshape the data to fit the R aov() function. thus instead of looking like this: subjectdrug strain exposure_1 exposure2 1

Re: [R] how to count the total number of (INCLUDING overl apping) occurrences of a substring within a string ?

Gabor Grothendieck ggrothendieck at gmail.com writes: Try this: findall(aba, ababacababab) [1] 1 3 7 9 gregexpr(a(?=ba), ababacababab, perl = TRUE) [[1]] [1] 1 3 7 9 attr(,match.length) [1] 1 1 1 1 findall(a.a, ababacababab) [1] 1 3 5 7 9 gregexpr(a(?=.a), ababacababab, perl =

Re: [R] how to count the total number of (INCLUDING overlapping) occurrences of a substring within a string?

Try this: findall((.).\\1, ababacababab) [1] 1 2 3 5 7 8 9 10 gregexpr((.)(?=.\\1), ababacababab, perl = TRUE) [[1]] [1] 1 2 3 5 7 8 9 10 attr(,match.length) [1] 1 1 1 1 1 1 1 1 On Sun, Dec 20, 2009 at 9:33 AM, Hans W Borchers hwborch...@googlemail.com wrote: Gabor

Re: [R] Problem reading binaries created with fortran

Duncan Murdoch helped me sort it out last Friday. Thanks to all who responded. On Fri, Dec 18, 2009 at 11:20 AM, Don MacQueen m...@llnl.gov wrote: However, source code is available at http://water.usgs.gov/nrp/gwsoftware/modflow2005/modflow2005.html so it would seem that the details are

Re: [R] simple main effect.

Thanks Or, So I am failing to understand, When you put that aov expression into the summary() Why isn't what you are getting what you need ? Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me:

Re: [R] simple main effect.

Also, do you have any missing data ? Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebrew) | www.r-statistics.com/ (English)

Re: [R] simple main effect.

I don't have missing data. about what I need. Lets say the drug*strain interaction is significant - now I want to check for drug under the levels of strain - compare drug 1 and 2 only on strain 1 and then only on strain 2. Or I'd like to compare the strains under levels of exposure. This is the

[R] expression()

Hi All, I'm wondering if its possible to write degree in symbol. I would like y-label as Temperature (degreeF). where degree should be in symbols. Thanks in advance, #R Code library(lattice) data(barley) barchart(yield ~ variety | site, data = barley, groups = year, layout = c(1,6),

Re: [R] simple main effect.

On 12/20/2009 04:56 PM, Or Duek wrote: I don't have missing data. about what I need. Lets say the drug*strain interaction is significant - now I want to check for drug under the levels of strain - compare drug 1 and 2 only on strain 1 and then only on strain 2. Or I'd like to compare the

Re: [R] expression()

Hi, try this, ylab = expression(Temperature~(degree*F)) ?plotmath baptiste 2009/12/20 Kim Jung Hwa kimhwamaill...@gmail.com: Hi All, I'm wondering if its possible to write degree in symbol. I would like y-label as Temperature (degreeF). where degree should be in symbols. Thanks in

Re: [R] simple main effect.

For some reasion I wasn't able to use TukeyHSD - I think because I need to set the different levels under a second variable. Tukey only helps me when I have more than 2 levels of same variable. Thanl you. On Sun, Dec 20, 2009 at 6:32 PM, S Devriese sdmaill...@gmail.com wrote: On 12/20/2009 04:56

[R] basic proto question

Dear list, I made the following example of a proto object that contains some data and a spline interpolation. I don't understand why test$predict() fails with this error message: Error: evaluation nested too deeply: infinite recursion / options(expressions=)? Best regards, baptiste test -

Re: [R] basic proto question

The free variables in a proto method are looked up in the object that the method was defined in so by referencing predict within test$predict you are referring back to test$predict whereas you mean to refer to stats::predict. Change the line that calls predict to: stats::predict(.$spline(),

[R] as.Date question

All! This piece of code: zzz1 - as.POSIXct(1999-03-18, tz=CET) zzz2 - as.POSIXlt(1999-03-18, tz=CET) zzz1 == zzz2 as.Date(zzz1) as.Date(zzz2) yields TRUE for zzz1==zzz2, but the two dates returned by as.Date are different: as.Date(zzz1) [1] 1999-03-17 as.Date(zzz2) [1] 1999-03-18 I'm using

Re: [R] simple main effect.

Hi Or, I understand your question (and am not sure what the confusion is about actually). Just like you said, you want to know the effect of A at B=1, and the effect of A at B=2. In this case, you want to know if drug has a significant effect for those with strain one, and whether drug has a

[R] Problem with sequence in loop

Hi all, I've got a list of files from 1914 to 2000. For each file, I can call the read.table function as follows. file - read.table(1914.txt) But I want to make a script so that at the end of the loop, the script selects the file 1915.txt, 1916.txt and so on. I've used seq() to create

Re: [R] basic proto question

Thanks, it seems so obvious now! baptiste 2009/12/20 Gabor Grothendieck ggrothendi...@gmail.com: The free variables in a proto method are looked up in the object that the method was defined in so by referencing predict within test$predict you are referring back to test$predict whereas you

Re: [R] as.Date question

Hello, On 2009.12.20 18:06:17, MAL wrote: All! This piece of code: zzz1 - as.POSIXct(1999-03-18, tz=CET) zzz2 - as.POSIXlt(1999-03-18, tz=CET) zzz1 == zzz2 as.Date(zzz1) as.Date(zzz2) yields TRUE for zzz1==zzz2, but the two dates returned by as.Date are different:

[R] as.Date question

All! This piece of code: zzz1 - as.POSIXct(1999-03-18, tz=CET) zzz2 - as.POSIXlt(1999-03-18, tz=CET) zzz1 == zzz2 as.Date(zzz1) as.Date(zzz2) yields TRUE for zzz1==zzz2, but the two dates returned by as.Date are different: as.Date(zzz1) [1] 1999-03-17 as.Date(zzz2) [1] 1999-03-18 I'm

Re: [R] simple main effect.

Thank you very much Ista. Can you please be a bit more specific as to using the overall error. I don't know how to actually do it in R. thank you, Or. On Sun, Dec 20, 2009 at 7:09 PM, Ista Zahn istaz...@gmail.com wrote: Hi Or, I understand your question (and am not sure what the confusion is

[R] Recall: Problem with sequence in loop

I recall the above problem. It works. Just error on my part in the file output. Here is the script that works. - seq - paste(seq(1914, 1916, by=1), *.y, sep=.) for (i in 1:3){ list - list.files(~/ukcp09/txt/x.djf, seq[[i]]) file - lapply(list, read.table) mean - (Reduce(+, file))

Re: [R] simple main effect.

Well, the F value is just MS_effect/MS_error. So you run the original model, run the subset models, but disregard the F values. To calculate the correct F values, divide MS_effect from the subset model by the MS_error term from the original (full) model. It's not really a matter of learning how to

[R] How to rotate an axis?

Hi everybody, I'm trying to build a trilineare coordinate system with R. For that, I need to rotate an axis. Does anybody know, if it is possible to rotate an axis, created with the command axis(), about for instance 60 degrees? I'm looking foward to any ideas and hints and want to wish

[R] plotting polynomial regression line

Dear All, I am trying to plot polynomial regression line to a scatterplot. I did following so far: x=c(1:9335) y=read.table(gp.txt,header=T,sep=\t) length(y$PCC) # y$PCC has values between 1 to 0 in decreasing order [1] 9335 plot(x,y$PCC,col=red) #scatterplot between x and y$PCC

Re: [R] How to rotate an axis?

On 20.12.2009 19:06, Etienne Stockhausen wrote: Hi everybody, I'm trying to build a trilineare coordinate system with R. For that, I need to rotate an axis. Does anybody know, if it is possible to rotate an axis, created with the command axis(), about for instance 60 degrees? You cannot

Re: [R] How to rotate an axis?

On Dec 20, 2009, at 1:06 PM, Etienne Stockhausen wrote: Hi everybody, I'm trying to build a trilineare coordinate system with R. For that, I need to rotate an axis. Does anybody know, if it is possible to rotate an axis, created with the command axis(), about for instance 60 degrees?

Re: [R] plotting polynomial regression line

On 20.12.2009 19:35, Amit wrote: Dear All, I am trying to plot polynomial regression line to a scatterplot. I did following so far: x=c(1:9335) y=read.table(gp.txt,header=T,sep=\t) length(y$PCC) # y$PCC has values between 1 to 0 in decreasing order [1] 9335 plot(x,y$PCC,col=red)

Re: [R] plotting polynomial regression line

Hello Amit, On 2009.12.20 19:35:09, Amit wrote: Dear All, I am trying to plot polynomial regression line to a scatterplot. I did following so far: x=c(1:9335) y=read.table(gp.txt,header=T,sep=\t) length(y$PCC) # y$PCC has values between 1 to 0 in decreasing order [1] 9335

Re: [R] birats example with R2WinBugs

On 18.12.2009 17:29, M JH wrote: Hello, I am trying, and failing, to do the birats example from the WinBugs manual with R2Winbugs. I can manage the rats example OK. Also, I can manage birats in WinBugs. Here is the code I am running, the .bug program, and the error message. Any help would

Re: [R] plotting polynomial regression line

Amit, please provide gp.txt file. On Sun, Dec 20, 2009 at 1:47 PM, Jason Morgan jwm-r-h...@skepsi.net wrote: Hello Amit, On 2009.12.20 19:35:09, Amit wrote: Dear All, I am trying to plot polynomial regression line to a scatterplot. I did following so far: x=c(1:9335)

Re: [R] plotting polynomial regression line

On Dec 20, 2009, at 1:35 PM, Amit wrote: Dear All, I am trying to plot polynomial regression line to a scatterplot. I did following so far: x=c(1:9335) y=read.table(gp.txt,header=T,sep=\t) length(y$PCC) # y$PCC has values between 1 to 0 in decreasing order [1] 9335 plot(x,y$PCC,col=red)

Re: [R] read.table: mysterious line omissions

Thanks everyone! The problem was definitely with the apostrophe being interpreted as a single quote. Jonathan On Sun, Dec 20, 2009 at 7:49 AM, jim holtman jholt...@gmail.com wrote: Most likely an unbalanced quote. put the following option in the read.table: quote='', comment.char='' On

[R] How to put text outside an xyplot?

Dear R-users, I have a plot created with the code below. I tried to put some text to the right of the color key. I worked with grid.text and also viewport(), but couldn't achieve this. I know this should be simple, but I just couldn't figure out how to do it. How does this work? Cheers,

Re: [R] How to rotate an axis?

Hi, Do you want a ternary plot? http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=34 It's easy to rotate an axis with Grid graphics, library(grid) pushViewport(viewport(0.5,0.5, width=0.5, height=unit(3, lines))) grid.xaxis(at=seq(-0.5,0.5,by=0.1), vp=viewport(x=1, angle=-60)) HTH,

[R] Object of type 'closure' not subsettable

Hi all, How can I overcome the error object of type 'closure' not subsettable I ran the following script seq - paste(seq(1914, 1916, by=1), *.y, sep=.) # make sequence c - 3 # total number of files d2 - file # creates dummy file # Input sequence in loop for (i in 1:3){ list -

Re: [R] Object of type 'closure' not subsettable

On Sun, Dec 20, 2009 at 7:40 PM, Muhammad Rahiz muhammad.ra...@ouce.ox.ac.uk wrote: Hi all, How can I overcome the error object of type 'closure' not subsettable I ran the following script seq - paste(seq(1914, 1916, by=1), *.y, sep=.) # make sequence c - 3 # total number of files d2 -

Re: [R] Object of type 'closure' not subsettable

On Dec 20, 2009, at 2:40 PM, Muhammad Rahiz wrote: Hi all, How can I overcome the error object of type 'closure' not subsettable I ran the following script seq - paste(seq(1914, 1916, by=1), *.y, sep=.) # make sequence c - 3 # total number of files d2 - file # creates dummy file # Input

Re: [R] Run time error when executing sqlQuery using the 64-bit version of R with 64-bit RODBC package in a Solaris 10 Sparc machine.

On Dec 19, 2009, at 11:39 AM, Alex Chelminsky wrote: I have compiled and linked a 64 bit version of R (R 2.9.2) and the corresponding unix ODBC 64 bit package The red highlighted text below is the error I'm getting trying to when invoking a sqlQuery library(RODBC) channel -

[R] Problems in installing MCMCglmm package

Dear R-Helpers, I am having troubles with installing with MCMCglmm package and I get the following error with a package Matrix Warning in library(pkg, character.only = TRUE, logical.return = TRUE, lib.loc = lib.loc) : there is no package called 'Matrix' Error: package 'Matrix' could not be

[R] insert a dot to the numbers

Hi, Anybody can give me some hints on the following problem? s-c(110,112321) I want to insert a dot . after the third number and get the following results. 110. 112.321 Thanks a lot. -- - Jane Chang Queen's [[alternative HTML version deleted]]

Re: [R] insert a dot to the numbers

On Dec 20, 2009, at 3:50 PM, rusers.sh wrote: Hi, Anybody can give me some hints on the following problem? s-c(110,112321) I want to insert a dot . after the third number and get the following results. 110. 112.321 Thanks a lot. s - c(110, 112321) s / (10 ^ (nchar(s) -

Re: [R] simple main effect.

For simple effects in the presence of interaction there are several options included in the HH package. If you don't already have the HH package, you can get it with install.packages(HH) Graphically, you can plot them with the function interaction2wt(..., simple=TRUE) See the examples in

Re: [R] forecasting

Very good patrick, lol. does anyone have any useful ideas Patrick Burns wrote: DispersionMap wrote: What about validating forecasting results. I have 5 years of data and have been forecasting the following three years i.e. 2010 to 2012. How can i check my forecast. Why types of

Re: [R] insert a dot to the numbers

On 20-Dec-09 21:50:58, rusers.sh wrote: Hi, Anybody can give me some hints on the following problem? s-c(110,112321) I want to insert a dot . after the third number and get the following results. 110. 112.321 Thanks a lot. - Jane Chang Queen's

[R] Object is not a matrix Error

I'm trying to follow this guide here: http://www.ats.ucla.edu/stat/r/modules/dummy_vars.htm In which I'm creating categorical variables using the factor function. I am able to go through the example listed above and have everything work, however, when I try to input my own numbers, I get an

Re: [R] Object is not a matrix Error

Where is the object 'write'? SHouldn't you be using: lm(visits ~ (day.f)) On Sun, Dec 20, 2009 at 5:59 PM, John Paul Telthorst jpteltho...@gmail.comwrote: I'm trying to follow this guide here: http://www.ats.ucla.edu/stat/r/modules/dummy_vars.htm In which I'm creating categorical

Re: [R] Object is not a matrix Error

When you imported did you not import the headers? Joe King 206-913-2912 j...@joepking.com Never throughout history has a man who lived a life of ease left a name worth remembering. --Theodore Roosevelt -Original Message- From: r-help-boun...@r-project.org

[R] CRAN (and crantastic) updates this week

CRAN (and crantastic) updates this week New packages Updated packages New reviews --- * deal, by newfuntek http://crantastic.org/reviews/45 * bnlearn, by newfuntek http://crantastic.org/reviews/44 * gRain, by newfuntek

Re: [R] insert a dot to the numbers

Try this also: sapply(s, function(x)format(x, big.interval = nchar(x) - 3, big.mark = .)) On Sun, Dec 20, 2009 at 7:50 PM, rusers.sh rusers...@gmail.com wrote: Hi,  Anybody can give me some hints on the following problem? s-c(110,112321)  I want to insert a dot . after the third number

[R] Limit on number of times Realloc can be called?

Hello, I've been writing a program in C that will be called by R. I seem to have stumbled upon an odd error that seems to suggest there is a limit on the number of times Realloc (the R version as defined in the manual R-extenstions not the C version realloc) when I try to use the following

Re: [R] Limit on number of times Realloc can be called?

Adam Waldemar Kowalewski wrote: Hello, I've been writing a program in C that will be called by R. I seem to have stumbled upon an odd error that seems to suggest there is a limit on the number of times Realloc (the R version as defined in the manual R-extenstions not the C version realloc) when

[R] Can you delete me from the list I don't want to receive the emails from the forum anymore, thank you!

Date: Mon, 21 Dec 2009 12:33:26 +1100 From: a.kowalew...@ugrad.unimelb.edu.au To: r-help@r-project.org Subject: [R] Limit on number of times Realloc can be called? Hello, I've been writing a program in C that will be called by R. I seem to have stumbled upon an odd error that seems

Re: [R] Can you delete me from the list I don't want to receive the emails from the forum anymore, thank you!

A couple of week back someone posted this. It applies to you too: None of the people reading this at the moment can do this for you. Read the information on the page where you subscribed.: https://stat.ethz.ch/mailman/listinfo/r-help Q: How many psychiatrists does it take to unscrew a light

[R] Signif. codes

My question is about the Signif. codes and the p-value, specifically, the output when I run summary(nameofregression.lm) So you get this little key: Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 And on a regression I ran, next to the intercept data, I get '***' Coefficients:

[R] Reading multiple Input Files

Dear R helpers,   Suppose I am dealing with no of interest rates at a time and the no of interest rates I am selecting for my analysis is random i.e. it can be 2, can be 10 or even higher. The R-code I had written (with the guidance of R helpers obviously and I am really grateful to all of you)

Re: [R] Signif. codes

There will always be uncertainty in your estimates so you don't have 0 percent chance of being wrong, but remember that's your intercept, your regressors are not significant. Although you can say it is less than ..05, I mean if its significant at .001 (or something like that), that's

Re: [R] Signif. codes

No, so the probability means that's the probability of getting that data by chance, so a p-value of .9997 means there is a .9997 probability that the data could be acquired by chance. This is a very simplistic view and you should study the regression model better. Joe King 206-913-2912

Re: [R] Hello all, How can I get corss-validation MSE of SVM in e1071?

Hello, Max The caret package is so good, I am learning it, but one problem is that nearZeroVar function can be used to identify near zero–variance variables and it only identify, how can I remove those variables that were identify, because I have many zero- or near zero- ones, it is not

[R] what is criterion of removing independence?

Hello, all I have a lot of independents and one dependent, finally, I want to build one model using them, and predict the new samples value, that is regression. before it, I must remove some independents according to some criterion: 1. constant values independent. 2. variant near zero.