[R] R-3.0.2 on September 25

2013-08-21 Thread peter dalgaard
Just a quick note, mainly to warn off maintainers of the recommended packages, that we intend to release R-3.0.2 on Wednesday, September 25. We'll be following the usual schedule from http://developer.r-project.org/release-checklist.html Notice in particular that new versions of recommended pack

Re: [R] Windows "The system cannot find the path specified" error when submitting R in batch mode

2013-08-21 Thread Tal Galili
Hello James, Does this answer your question: http://stackoverflow.com/questions/14956887/problems-executing-script-from-command-line-in-r-error-message-cannot-find-pat ? Contact Details:--- Contact me: tal.gal...@gmail.com | R

Re: [R] italic(mu)

2013-08-21 Thread David Winsemius
On Aug 20, 2013, at 3:04 PM, David Winsemius wrote: > > On Aug 20, 2013, at 12:06 PM, Beaulieu, Jake wrote: > >> Hi, >> >> I would like to include the greek letter mu, in italics, in an axis title. >> The following gets close, but the mu isn't italicized. >> >> x <- 1:10 >> y <- 1:10 >> plo

Re: [R] download.file error - corrupt: Can't read SAT; charset=binary'

2013-08-21 Thread Tal Galili
Dear Uwe, My apologies - my original code had mode="wb", which also *does not* work and produce the errors I've mentioned in the previous e-mail. The code I pasted had mode ="w", since that is one of the versions I've played with (which also does not work). With regards, Tal On Wed, Aug 21,

[R] cross-sectional analysis of a financial time series

2013-08-21 Thread Brijesh Gulati
Hi: I have a financial series data. For instance, one can take YHOO from the quantmod package. >library("quantmod") >getSymbols("YHOO") As you can see this series has date along with close prices. I want to do a cross-sectional analysis of the time series and want to see if there are any sea

Re: [R] R will not find certain entries in an array created by an interval

2013-08-21 Thread arun
Hi, Please check: ?FAQ 7.31 http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f  intervalloweredge[3]-3.845 #[1] -4.440892e-16 which(as.character(intervalloweredge)=="3.845") #[1] 3 which(as.character(intervalloweredge)=="3.86") #[1] 4 #or which(rou

Re: [R] Replacing data values from ref table

2013-08-21 Thread arun
Hi, You may try ?merge() or ?join() from library(plyr) dat1<- read.table(text=" Col1 Col2 red 1.23 blu 4.56 ",sep="",header=TRUE,stringsAsFactors=FALSE) dat2<- read.table(text=" Col1 Col2 red SQU blu CIR ",sep="",header=TRUE,stringsAsFactors=FALSE) library(plyr) res1<-join(dat2,dat1,by="Col1")[

[R] Lag selection in Unit Root Test

2013-08-21 Thread jpm miao
Hi, I attempt to conduct the unit root test with the two functions unitrootTest(x, lags = 1, type = c("nc", "c", "ct"), title = NULL, description = NULL) adfTest(x, lags = 1, type = c("nc", "c", "ct"), title = NULL, description = NULL) According to the document, "lagsis the maxim

[R] Interpreting the result of 'cutree' from hclust/heatmap.2

2013-08-21 Thread Gundala Viswanath
I have the following code that perform hiearchical clustering and plot them in heatmap. __ library(gplots) set.seed(538) # generate data y <- matrix(rnorm(50), 10, 5, dimnames=list(paste("g", 1:10, sep=""), paste("t", 1:5, sep=""))) # the actual data is much larger that the above # perform hiear

Re: [R] how to code y~x/(x+a) in lm() function

2013-08-21 Thread Ben Bolker
On 13-08-21 05:17 PM, Rolf Turner wrote: > > > Thott about this a bit more and have now decided that I don't understand > after all. > > Doesn't > > glm(1/y~x,family=gaussian(link="inverse")) > > fit the model > > 1/y ~ N(1/(a+bx), sigma^2) > > whereas what the OP wanted was > >

Re: [R] X axis label as months

2013-08-21 Thread Jim Lemon
On 08/22/2013 09:04 AM, Prof J C Nash (U30A) wrote: Thanks. The staxlab works with the plot generated by plot(tt, data,) but not with plot(tsdata, ...) i.e., it seems to be the tsplot that somehow changes something and the axis commands have no effect. I probably should have pointed out

Re: [R] Narrowing values collected from .txt file

2013-08-21 Thread Morway, Eric
Hello John, I tried the dput function a couple of different ways, with some different arguments set, and it dumped the following (or less): "structure(3L, class = c("textConnection", "connection"), conn_id = )" Concerned that I miscopied some of the reproducible example to my original post, I ju

Re: [R] R.oo in R-3.0.1

2013-08-21 Thread Anamika Chaudhuri
Hi All: I am getting some warnings trying to run the following code and I am not sure what they mean. > set.seed(111) > k<-63 > n<-60 > x<-NULL > p<-rbeta(k,3,3)# so that the mean nausea rate is alpha/(alpha+beta) > #min<-10 > #max<-60 > #n<-as.integer(runif(k,min,max)) > for(i in 1:k) + x<-cbind(

Re: [R] X axis label as months

2013-08-21 Thread Prof J C Nash (U30A)
Thanks. The staxlab works with the plot generated by plot(tt, data,) but not with plot(tsdata, ...) i.e., it seems to be the tsplot that somehow changes something and the axis commands have no effect. I probably should have pointed out that my main concern is that I'm not getting any e

Re: [R] X axis label as months

2013-08-21 Thread Jim Lemon
On 08/22/2013 07:56 AM, Prof J C Nash (U30A) wrote: There are several items on the web about putting month names as tick labels on x axis of time plots, but I found when I tried them they did not work for me. After an hour or so of driving myself silly looking for a bug in my code, I've prepared

Re: [R] X axis label as months

2013-08-21 Thread Gabor Grothendieck
On Wed, Aug 21, 2013 at 5:56 PM, Prof J C Nash (U30A) wrote: > There are several items on the web about putting month names as tick labels > on x axis of time plots, but I found when I tried them they did not work for > me. After an hour or so of driving myself silly looking for a bug in my > code

[R] X axis label as months

2013-08-21 Thread Prof J C Nash (U30A)
There are several items on the web about putting month names as tick labels on x axis of time plots, but I found when I tried them they did not work for me. After an hour or so of driving myself silly looking for a bug in my code, I've prepared a reproducible example below, where I did find a w

[R] R.oo in R-3.0.1

2013-08-21 Thread Grach, Alexander
Fellow R users, I'm trying to upgrade from 2.11 to 3.0.1. I have a legacy S3 class defining method "DAILY" within one of our "in-house" packages using R.oo. This method does not seem to get 'resolved' under 3.0.1, specifically the following does not work: CODE WITHIN PACKAGE (e.g. MyPackage

Re: [R] ave function

2013-08-21 Thread arun
Hi Robert, source("shareB101") ##Clean is the dataset  res1<-with(Clean,aggregate(GRADE,list(TERM,INST_NUM),FUN=function(x) cbind(shapiro.test(x)$p.value,shapiro.test(x)$statistic)) )  head(res1) #  Group.1 Group.2  x.1  x.2 #1  201001  689809 1.720329e-07 9.307362e-01 #2  201201

Re: [R] download.file error - corrupt: Can't read SAT; charset=binary'

2013-08-21 Thread Uwe Ligges
On 21.08.2013 22:35, Tal Galili wrote: Dear R-help mailing list memebers, I encountered the following error, and I'd be happy for ideas on how to fix it: # using R 3.0.1 on Windows 7: exe_URL = 'http://pandoc.googlecode.com/files/pandoc-1.11.1.msi' exe_filename <- file.path(tempdir(), file.n

[R] download.file error - corrupt: Can't read SAT; charset=binary'

2013-08-21 Thread Tal Galili
Dear R-help mailing list memebers, I encountered the following error, and I'd be happy for ideas on how to fix it: # using R 3.0.1 on Windows 7: exe_URL = 'http://pandoc.googlecode.com/files/pandoc-1.11.1.msi' exe_filename <- file.path(tempdir(), file.name.from.url(exe_URL)) download.file(exe_UR

Re: [R] Coordinate scales for pairs plot

2013-08-21 Thread Ted Harding
On 21-Aug-2013 19:08:29 David Winsemius wrote: > > On Aug 21, 2013, at 10:30 AM, (Ted Harding) wrote: > >> Greetings all. >> >> I suspect this question has already been asked. Apologies >> for not having taced it ... >> >> In the default pairs plot produces by the function pairs(), >> the coord

Re: [R] lattice: bwplot: getting two plots into one.

2013-08-21 Thread Anna Zakrisson Braeunlich
Dear Mr. Heiberger, thank you for your reply. I think that you may have misunderstood my question. What I want is to superpose the two panels into one. Lets take this simplified code: bwplot(var1 ~ factor1|factor2, data=mydata, panel=panel.bwplot.intermediate.hh) What I want is this: bwplot(var1

Re: [R] rcode problem

2013-08-21 Thread John Kane
Your data arrived but no code. The R help list is fussy about what type of files you can attache. For an alternative approach have a look at http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example John Kane Kingston ON Canada > -Original Message- > From

Re: [R] Narrowing values collected from .txt file

2013-08-21 Thread John Kane
I am not sure what happened but it may be that you accidentally sent your message in html but the text connection data seems unuseable. It is much better to use ?dput to supply sample data. Have a look at http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example Joh

Re: [R] Coordinate scales for pairs plot

2013-08-21 Thread David Winsemius
On Aug 21, 2013, at 10:30 AM, (Ted Harding) wrote: > Greetings all. > > I suspect this question has already been asked. Apologies > for not having taced it ... > > In the default pairs plot produces by the function pairs(), > the coordinate scales alternate between top and bottom and > right an

Re: [R] lattice: bwplot: getting two plots into one.

2013-08-21 Thread Richard M. Heiberger
That is also easy. See the ## install.packages("HH") ## if you don't yet have it. library(HH) demo("bwplot.examples", package="HH") These examples show how to use the tools based on panel=panel.bwplot.intermediate.hh and you will easily be able to modify the code for your specific case. Rich

[R] Coordinate scales for pairs plot

2013-08-21 Thread Ted Harding
Greetings all. I suspect this question has already been asked. Apologies for not having taced it ... In the default pairs plot produces by the function pairs(), the coordinate scales alternate between top and bottom and right and left sides. For example, in a 5x5 plot for variables X1, X2, X3, X

Re: [R] data import: strange experience

2013-08-21 Thread David Carlson
You should be able to figure it out if you just print out the four factor levels that read.table() missed. The main differences are that read.table() includes ' in the quote= argument and it recognizes # as a comment (and therefore discards it and everything after it): setdiff(levels(dfcsv$Var), l

Re: [R] Problems installing FAME package on Windows 7.

2013-08-21 Thread Uwe Ligges
Probably something with your environment and/or your toolset. Does the CRAN binary work for you? i.e. just use install.packages("fame") and try it out.. Best, Uwe Ligges On 20.08.2013 11:07, Andreas Dibiasi wrote: I have problems to install the FAME package. I have installed the fame cli

Re: [R] how to code y~x/(x+a) in lm() function

2013-08-21 Thread Ben Bolker
Rolf Turner xtra.co.nz> writes: > > On 21/08/13 11:23, Ye Lin wrote: > > T > > hanks for your insights Rolf! The model I want to fit is y=x/a+x with > > no intercept, so I transformed it to 1/y=1+a/x as they are the same. > > For crying out loud, they are ***NOT*** the same. The equations y =

Re: [R] odfWeave post processing error

2013-08-21 Thread Tom Li
> Error: XML content does not seem to be XML: 'content_1.xml'In addition: > Warning message:In file.remove("content.xml") : cannot remove file > 'content.xml', reason 'No such file or directory' Check out the XML file. When I run into errors with odfWeave, it's usually because a non-unicode

Re: [R] lattice: bwplot: getting two plots into one.

2013-08-21 Thread Richard M. Heiberger
This is the fundamentals bwplot(var1 ~ factor1 | factor2, data=mydata, panel=function(...) { panel.abline(h=0, lty=2, col="gray90") panel.bwplot(...) }) You can add the font control if you need to. The major change I made was to put the panel.abline first. It is part of the background

Re: [R] multiple core

2013-08-21 Thread Uwe Ligges
On 20.08.2013 21:05, Andrea Lamont wrote: Hello: I have somewhat of a generic question. How do you use multiple CPU/core processors in R? For those familiar with Mplus software, I'm looking for something analogous to PROCESSORS command in Mplus to speed up processing. I searched previous h

Re: [R] data import: strange experience

2013-08-21 Thread SH
Thanks Peter. It works with read.delim. David: Thanks for your comments. To answer your questions. I don't have 'NA' and all balanced. The number of mssing levels were 4 and it happened only to those four levels. Yes, there is commas embedded and some characters (e.g., '-', space, some wired

Re: [R] data import: strange experience

2013-08-21 Thread SH
Hi Sarah, Thanks for a prompt feedback. I knew it will be very vague without example. However, I only used two commands to import data and had no 'apparent' errors. The original data have about 19000 obs and I was able to reduce about 3200. I wonder if I can attach the data file (size: 109K) w

Re: [R] data import: strange experience

2013-08-21 Thread David Carlson
This is not really enough information to diagnose the problem. What are the missing factor levels? Were the missing levels combined with another level or do you have missing values (NA) for those observations? Do the extra factor levels include embedded commas? There are differences between read.ta

Re: [R] Bringing tk_choose.files to front

2013-08-21 Thread mmiikkii
Dear Jimbo, Unfortunately I got no answer to this problem so I decided to put the instructions in user manual for the program I am writing...not the best solution, but I cannot think of anything better. If you find out how to solve this in a better way, please share the knowledge. Best

Re: [R] data import: strange experience

2013-08-21 Thread peter dalgaard
On Aug 21, 2013, at 16:46 , Sarah Goslee wrote: > Hi, > > We don't know anything about your data or your file, so it's utterly > impossible to offer useful suggestions. > > The very best thing you can do is condense your problem into a > reproducible example, with fake data if necessary. Otherw

Re: [R] Filter Dataframe for Alarm for particular column(s).

2013-08-21 Thread R_Antony
Ok. I got it. From: Akkara, Antony (GE Energy, Non-GE) Sent: Wednesday, August 21, 2013 3:52 PM To: 'arun kirshna [via R]' Subject: RE: Filter Dataframe for Alarm for particular column(s). What will be the code, if we required all the row number ? Thanks, Antony. From: arun kirsh

[R] lattice: bwplot: getting two plots into one.

2013-08-21 Thread Anna Zakrisson Braeunlich
Dear all, I have problems gettting the following as one graph instead of the two created by this script. I know that with this dummy data, the resulting graph will look like a catastrophy. I have other data I wish to apply this too. Just wanted to supply some dummy data to work with. There may

Re: [R] Filter Dataframe for Alarm for particular column(s).

2013-08-21 Thread R_Antony
What will be the code, if we required all the row number ? Thanks, Antony. From: arun kirshna [via R] [mailto:ml-node+s789695n467097...@n4.nabble.com] Sent: Saturday, July 06, 2013 12:54 AM To: Akkara, Antony (GE Power & Water, Non-GE) Subject: Re: Filter Dataframe for Alarm for particula

Re: [R] data import: strange experience

2013-08-21 Thread Sarah Goslee
Hi, We don't know anything about your data or your file, so it's utterly impossible to offer useful suggestions. The very best thing you can do is condense your problem into a reproducible example, with fake data if necessary. Otherwise you're limited by the ability of the list to guess what you'

[R] data import: strange experience

2013-08-21 Thread SH
Dear List: I had some strange experience in importing data. I wonder if anyone of you had the same problem before and would greatly appreciate your suggestion in advance. The original data set in excel format. Here is a brief summary of the procedure I did: 1. I saved the original excel data as

[R] Narrowing values collected from .txt file

2013-08-21 Thread Morway, Eric
useRs- The output generated from a groundwater model post-processor contains millions of lines of text. Using the custom R function shown below, I can quickly gather values from this file. As you can see in the textConnection provided below (which is only a small snippet from the file), the outp

Re: [R] Rmpi unable to load

2013-08-21 Thread Prof Brian Ripley
I think you need to ask whoever installed Rmpi on your machine or provided a binary package: CRAN does not distribute a Windows binary, and one reason is that Rmpi has to be compiled against a specific MPI installation. Or ask the package maintainer, as the posting guide suggested. In any ca

[R] Rmpi unable to load

2013-08-21 Thread Stephen Clark
Hello, I am trying to run the optimise_ppso_robust parallel swarm optimisation function in R3.0.0 on a 6 core AMD Windows 8 64bit PC. To do this I need to load the Rmpi package. When I do this I get the following message in a dialogue box: The program can't start because mpich2mpi.dll is missi

Re: [R] Visualising SOM as hexagonal heatmaps

2013-08-21 Thread Roger Bivand
Adam Szałucha gmail.com> writes: > > I have a kohonen object generated with grid parameters grid = > somgrid(20, 20, "hexagonal"). Standard plot shows the codebook vectors > for all variables. I want to have a hexagonal heatmap for each > variable. > > I found in Applied Spatial Data Analysis w

[R] bootstrap from binary confidence intervals (Hmisc binconf)

2013-08-21 Thread Folkes, Michael
Hello all, I've generated a series of CI's of proportions from numerous samples using Hmisc binconf(). Now I wish to bootstrap from those values for some additional products. (ie sample to get a proportion, multiply that proportion with another number). We know it's not normally distributed. How

Re: [R] Confidence Intervals for survreg (survival)

2013-08-21 Thread Luis Antunes
Very simple, indeed. Thank you very much for your help. Luis > Subject: Re: [R] Confidence Intervals for survreg (survival) > From: dwinsem...@comcast.net > Date: Tue, 20 Aug 2013 11:50:52 -0700 > CC: r-help@r-project.org > To: ljantu...@hotmail.com > > > On Aug 20, 2013, at 7:04 AM, Luis Antun