Re: [R] Average curve in R

2012-06-02 Thread Rui Barradas
Hello, Just to add a note, since R defaults the x axis to the integer sequence 1:length(what.to.plot), this would do it: plot(rowMeans(datQ), type="l") If the op wants to plot this mean values line together with the other 4, one of the most forgotten plot instructions is the matrix plot inst

Re: [R] how to add a 'label' column

2012-06-02 Thread Jim Lemon
On 06/02/2012 05:47 AM, Jason Love wrote: Hello R users, I'd like to ask a question about how to add a new column. So, below is my situation. In order to perform the repeated ANOVA, I first imported the following table. score=read.csv("patients_tests.csv"); subject test1 test2

Re: [R] Add rank column to data frame as in SQL...

2012-06-02 Thread Gabor Grothendieck
On Fri, Jun 1, 2012 at 11:23 AM, Brigid Mooney wrote: > Hopefully this is an easy problem... > > I'm trying to add a partitioned rank column to a data frame where the > rank is calculated separately across a partition by categories, the > way you could easily do in SQL.  I found this solution in t

Re: [R] how to add a 'label' column

2012-06-02 Thread Jason Love
Thanks all for the helpful tips. Jason On Sat, Jun 2, 2012 at 7:14 AM, Jim Lemon wrote: > On 06/02/2012 05:47 AM, Jason Love wrote: > >> Hello R users, >> I'd like to ask a question about how to add a new column. So, below is my >> situation. >> >> In order to perform the repeated ANOVA, I firs

[R] How can I export a paired t-test output table to an excel file?

2012-06-02 Thread Jason Love
Hi R users, Could anyone let me know how to export a paired t-test output table (see below) to an excel file? Jason, with(score2,pairwise.t.test(values,ind, + p.adjust.method="holm", paired=T)) Pairwise comparisons using paired t tests data: values and ind test1

Re: [R] How can I export a paired t-test output table to an excel file?

2012-06-02 Thread Nicolas Iderhoff
Hello Jason, if you put your t-test into an r-object, like here (taken from the documentary example): t.test <- pairwise.t.test(Ozone, Month) you can get the table with t.test$p.value also you can get a lot of other stuff. you can see what's available with "str(t.test)". hope that's what you'

Re: [R] Silencing the output of install.packages()

2012-06-02 Thread Prof Brian Ripley
On 01/06/2012 08:46, Uwe Ligges wrote: On 01.06.2012 07:17, Tejas Kale wrote: Dear Uwe Many thanks for your reply. I agree with you but I need the silencing of output for a particular reason. I am working on a statistical package called VOStat which uses a Java based GUI to get the data and

Re: [R] Converting sp to maps ?/Map projection in sp-package?

2012-06-02 Thread Roger Bivand
Jari Haukka helsinki.fi> writes: > > Dear All, > > I hava maps in sp fomat and I like to use different map projection. > > SHouls I somehow convert sp to maps and make projection and ploting there? Please install the rgdal package, and use the spTransform() methods there to change the project

Re: [R] How can I export a paired t-test output table to an excel file?

2012-06-02 Thread David Winsemius
On Jun 2, 2012, at 8:58 AM, Nicolas Iderhoff wrote: Hello Jason, if you put your t-test into an r-object, like here (taken from the documentary example): Would need to attach(airqualty) for that to succeed t.test <- pairwise.t.test(Ozone, Month) Better to use an object name that is not t

Re: [R] raw input prompt

2012-06-02 Thread David Winsemius
On Jun 1, 2012, at 2:01 PM, meatloafthefrog wrote: Ah. I'm very new to this and I definitely wasn't clear enough about what I'm trying to do... sorry. I am extracting information (in this case, DNA sequences to be exported to a fasta file) from just two columns of a table that has a lot of

[R] mgcv (bam) very large standard error difference between versions 1.7-11 and 1.7-17, bug?

2012-06-02 Thread Martijn Wieling
Dear useRs, I reran an analysis with bam (mgcv, version 1.7-17) originally conducted using an older version of bam (mgcv, version 1.7-11) and this resulted in the same estimates, but much lower standard errors (in some cases 20 times as low) and lower p-values. This obviously results in a larger s

Re: [R] raw input prompt

2012-06-02 Thread J Toll
On Fri, Jun 1, 2012 at 1:01 PM, meatloafthefrog wrote: > I want a prompt that asks the user for the file name of the table, > the response to which will be the name of the table in R. > > So I did something like: > > file_name = function()                                         #This > function,

Re: [R] raw input prompt

2012-06-02 Thread Bert Gunter
... or On Sat, Jun 2, 2012 at 11:29 AM, J Toll wrote: > On Fri, Jun 1, 2012 at 1:01 PM, meatloafthefrog > wrote: >> I want a prompt that asks the user for the file name of the table, >> the response to which will be the name of the table in R. >> >> So I did something like: >> >> file_name = fun

[R] Double-buffering problem, this time with an example.

2012-06-02 Thread Daniel Carr
Most of my animations that used to work on windows() version 2.11.1 and earlier now flash as if the double buffering is turned off or buffer swapping is triggered by other events than in the past. The simplified example below using symbols should illustrate the problem in a windows environment.

Re: [R] plot background - excel gradient style background ?

2012-06-02 Thread Greg Snow
The short answer is "yes". The question as answer is "what is wrong with your data that you feel the need to hide/distort information and distract from the story of the data?". The hopefully thought provoking answer is "fortune(197)". The answer to the question not asked that should have been is

Re: [R] Determining frequency and period of a wave

2012-06-02 Thread Greg Snow
You could use the 'nls' function to fit a sine (or cosine) function to the data. On Fri, Jun 1, 2012 at 7:54 PM, Aaron Patterson wrote: > Hello!  I'm collecting data on a refrigerator that I'm using to cure > meat.  Specifically I am collection humidity and temperature readings. > The temperature

Re: [R] Double-buffering problem, this time with an example.

2012-06-02 Thread Bert Gunter
For the record, I confirm Dan's description on Windows: > sessionInfo() R version 2.15.0 (2012-03-30) Platform: x86_64-pc-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=English_United States.1252 [2] LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 [4] LC_NUMERIC=C [5]

Re: [R] mgcv (bam) very large standard error difference between versions 1.7-11 and 1.7-17, bug?

2012-06-02 Thread Simon Wood
the fREML results must be a bug. I'd modified the bam covariance and edf code in 1.7-17 and it looks like I must have messed something up. Any chance you could send me the data off line? Also can you try summary(...,freq=FALSE) under 1.7-17 and see if you get the same results as before? On 0

Re: [R] Double-buffering problem, this time with an example.

2012-06-02 Thread Duncan Murdoch
On 12-06-02 4:02 PM, Daniel Carr wrote: Most of my animations that used to work on windows() version 2.11.1 and earlier now flash as if the double buffering is turned off or buffer swapping is triggered by other events than in the past. The simplified example below using symbols should illustrat

[R] mtext bug

2012-06-02 Thread Matthew Johnson
Sir, I have hit a the limits of my understanding of text / par / opar etc... I have a few related xts data frames which have multiple columns, and have written a for-loop to make a set of charts - however i do not get any margin text when i run the loop. when i test the code outside of the loop

Re: [R] mtext bug

2012-06-02 Thread Duncan Murdoch
On 12-06-02 6:21 PM, Matthew Johnson wrote: Sir, I have hit a the limits of my understanding of text / par / opar etc... I have a few related xts data frames which have multiple columns, and have written a for-loop to make a set of charts - however i do not get any margin text when i run the lo

Re: [R] mtext bug

2012-06-02 Thread Matthew Johnson
thanks very much - fixed. On 3 June 2012 08:26, Duncan Murdoch wrote: > On 12-06-02 6:21 PM, Matthew Johnson wrote: > >> Sir, >> >> I have hit a the limits of my understanding of text / par / opar etc... >> >> I have a few related xts data frames which have multiple columns, and have >> written

Re: [R] Missing times chron

2012-06-02 Thread Gabor Grothendieck
On Fri, Jun 1, 2012 at 11:02 AM, Rick Admiraal wrote: > Dear all, > > As a novice user of R I ran into a problem that's quite hard for me to > resolve. I have a database containing data of a clinical trial in which > patients are included that survived or died: > > x <- matrix(data=c(1:5,0, "1/1/2

Re: [R] mtext bug

2012-06-02 Thread Matthew Johnson
Sir, I am now having a similar problem with the barplot: i cannot get the margin text to 'stick'. png(paste(drive, "MoMUnit.png", sep="")) par(mar=c(10, 3, 2, 2), opar=(1,0,0,0)) mainnameUnit = expression(paste("MoM %", Delta, " Apptmnt Prices", sep="")) ymin <- min(percentDiff[,c(6, 9, 33, 36, s

[R] multiple variance structure in lmer giving zero variances

2012-06-02 Thread Sara Krause
Hi all, I’m hoping someone might be able to help me out. Forgive me if my mistake is something simple. I am new to mixed models, new to R, and new to lme4 and am struggling to figure everything out. I have two questions that I am hoping someone can answer. 1) Am I using the correct random struc

[R] R 2.15.0 and 2.14.2 crash on the same code that runs on 2.14.1

2012-06-02 Thread Ebrahim Jahanshiri
Dear list, My code runs on 2.14.1 smoothly (with no error) but makes R crash on 2.14.2 and 2.15.0 with the usual windows message about "R for windows GUI front-end has stopped working..." and as such I dont have any access to R to get any error or warning message (if any). It crashes on random par

Re: [R] Adding a column into the file

2012-06-02 Thread John
On Fri, 1 Jun 2012 09:27:58 -0700 (PDT) pigpigmeow wrote: > Dear all, > I have a lot of problems on R-programming. > for example > my csv. file is .. > Date wrfRH wrfsolar wrfwindspeed wrfrain wrftd wrfta > 21/10/2010 92.97 22.11 53.27 0 1546.337861 61.00852664 > 22/10/2010 87.35

[R] Add a dim to an array

2012-06-02 Thread Patrick Hausmann
Dear list, I'm trying to add a new dim to a multidimensional array. My array looks like this a1 <- array(1:8, c(2, 2, 2)) dimnames(a1) <- list(A = c("A1", "A2"), B = c("B1", "B2"), D = c("D1", "D2")) I would like to add a new dim 'group' with the valu

Re: [R] Average curve in R

2012-06-02 Thread QAMAR MUHAMMAD UZAIR
Arun and Rui Barradas: thankyou very much for sorting out my issue. I got what I wanted. Rainer Schuermann and Jeff Newmiller: thankyou very much for your suggestion. I will take care of it in future. regards - On Fri, 1 Jun 2012 19:47:05 -0700 (PDT) arun wr

[R] Log-normal probability plot

2012-06-02 Thread eliza botto
Dear R users, You can literally safe my life my telling me the solution of my problem. I have created matrix of a data frame with 3 columns, with each column representing data of different year. 2 1 5 3 4 4 4 3

[R] Memory-random forest

2012-06-02 Thread MYRIAM TABASSO
Hi all, I am working with random forest on the large data set( 39). I have a warning message : "In matrix(integer(nrnodes * nt), ncol = nt) : Reached total allocation of 3766Mb: see help(memory.size)". Can you have any suggestions? Thanks, Myriam [[alternative HTML version deleted]] _

[R] Compare data between two groups/countries on 5-point Likert scale questionnare?

2012-06-02 Thread Frankfoot
Hi everyone, I'm stuck on my dissertation which is due next week. I'm a business major student and my project is a comparative research on corporate social responsibility (CSR) between Chinese and German firms. According to my literature review, it's quite obvious that German firms have much better

Re: [R] points() function will not plot

2012-06-02 Thread Sarah Goslee
Looks like you've swapped x and y. Sarah On Friday, June 1, 2012, David Winsemius wrote: > > On Jun 1, 2012, at 4:09 PM, jcrosbie wrote: > > I can not get the point function to work with in my code. >> >> The code: >> >> AggOfTempMatrix$**CumSumPercentSize<-c( >> 0.05265450, 0.05738490, 0.0586

[R] igraph and igraph0

2012-06-02 Thread jalantho...@verizon.net
Could someone tell me the difference between igraph and igraph0? I searched the CRAN web site, but cannot find an explantion of the differences. Thanks, Alan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://st

Re: [R] R and igraph problems

2012-06-02 Thread jalantho...@verizon.net
Why is igraph 0.5.5 not on the igraph sourceforgwe home page? [[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/posti

Re: [R] Help with this web scrape function

2012-06-02 Thread Sven D
Would like to see this at the top again, as I think it did not get pick up when I posted it initially. Apologies, last try. -- View this message in context: http://r.789695.n4.nabble.com/Help-with-this-web-scrape-function-tp4632137p4632169.html Sent from the R help mailing list archive at Nab

Re: [R] About "Plot.new"

2012-06-02 Thread czar
@mr.parashar4 - If this not resolved yet, here is what can help plot.new is trying to create a new object out of the standard point plot. I believe you have not created a standard point plot, i.e. in your case, plot(xm,ym). This will show the graphical plot window. Keep that window open and then ru

Re: [R] Add a dim to an array

2012-06-02 Thread Rui Barradas
Hello, Try a3 <- array(dim=c(2, 2, 2, 2)) dn <- dimnames(a1) dn$group <- c("low", "high") dimnames(a3) <- dn a3[] <- a1 # to use [] keeps the dimensions a3 # it fills a3 recycling a1 (two copies) dim(a2) dim(a3) As you can see, 'a2' is not of the right dimensions, and 'a3' has now

[R] R 15.0 Warning message.

2012-06-02 Thread Jinyan Huang
Dear list, I have install R15.0 on Ubuntu system using the source code. It is install to my home directory. I have no root privilege. When I start R, I got this warning message: During startup - Warning message: Setting LC_CTYPE failed, using "C" How can I fix this? Thanks. __

Re: [R] Add a dim to an array

2012-06-02 Thread Petr Savicky
On Sat, Jun 02, 2012 at 05:17:36PM +0200, Patrick Hausmann wrote: > Dear list, > > I'm trying to add a new dim to a multidimensional array. My array looks > like this > > a1 <- array(1:8, c(2, 2, 2)) > dimnames(a1) <- list(A = c("A1", "A2"), > B = c("B1", "B2"), >

[R] Modify pdf file with Illustrator.

2012-06-02 Thread Jinyan Huang
Dear list, Some times I want to modify the pdf file which produced by R with illustrator. But when I use Illustrator open the pdf file, it often makes the pdf some little changed. Anyone have some suggestions? Is it better to use other file type, not pdf? Or when I produced the pdf, I should set s