Re: [R] Can not find lme

2012-08-08 Thread Pascal Oettli
Hello, Function "lme" is provided by the "nlme" package. HTH Pascal Le 08/08/2012 15:22, Santini Silvana a écrit : Dear all, Can anyone help me, my R software can not run a nested linear regression by using the lme funcion. The message that appears is Error: could not find function "lme" I a

Re: [R] process evaluation packages (slightly off topic)

2012-08-08 Thread Petr PIKAL
Thanks Bert However this package is quite simple and actually it is just only plotting routine for already prepared list. This can be done by reading appropriate file but need a step before - prepare timeline from oriented graph which I did not find yet and the only way I know is pencil/paper a

[R] RQuantLib: SET_VECTOR_ELT() can only be applied to a 'list', not a 'symbol'

2012-08-08 Thread Cren
# Hi all, # trying to run the following example code # from 'RQuantLib' package... HullWhite <- list(term = 0.055, alpha = 0.03, sigma = 0.01, gridIntervals = 40) Price <- rep(as.double(100),24) Type <- rep(as.character("C"), 24) Date <- seq(as.Date("2006-09-15"), by = '3 months

Re: [R] Decimal number

2012-08-08 Thread Petr PIKAL
Probably a problem in your setting or envirenment or print function Here is what I get > x1<-64.90 > x2<-17.7025 > c(x1,x2) [1] 64.9000 17.7025 > x<-c(x1,x2) > x [1] 64.9000 17.7025 > Regards Petr > > HI > > >i have a little problem please help me to solve it > > >this is the code in R: > >

Re: [R] RQuantLib: SET_VECTOR_ELT() can only be applied to a 'list', not a 'symbol'

2012-08-08 Thread Cren
Cren wrote > > # trying to run the following example code > # from 'RQuantLib' package... > # Obviously, run require(RQuantLib) # before executing the example :) -- View this message in context: http://r.789695.n4.nabble.com/RQuantLib-SET-VECTOR-ELT-can-only-be-applied-to-a-list-not-a-symb

Re: [R] Can not find lme

2012-08-08 Thread Jeff Newmiller
Perhaps you should read the help for that package, where you would discover that it does not contain an "lme" function. library(help="lme4") You could type help.search("lme") and discover that there is an "lme" function in the "nlme" package. At that point you could perhaps confirm what you

Re: [R] Problem with global variable building a package

2012-08-08 Thread Eva Prieto Castro
Hello, Greg: That is one of the solutions I have tried; the problem is that I can't obtain a package in this way, because an glo-internal.R file is generated (glo is the name of the package) and when I try to create the package an error message appears (a parse error). The content of the glo-i

[R] which R version to connect with .net?

2012-08-08 Thread agent dunham
Dear community, I want to call R from .net. I've been reading this article: http://www.codeproject.com/Articles/25819/The-R-Statistical-Language-and-C-NET-Foundations. It says that first I need to download the COM server. When I do this, I see that the COM server is from 2008; http://cran.es.r

Re: [R] Problem with global variable building a package

2012-08-08 Thread Eva Prieto Castro
Hi Michael, I did not know those alternatives you mention (setCompilerOptions() and setRmetricsOptions()), so I will study them (apart to find out why I can not create a package when I use a new environment, as I was telling Greg). Thanks again for your help. Regards. Eva --- El mar, 7/8/1

Re: [R] label_wrap_gen question

2012-08-08 Thread vd3000
Hi, Brian So, I could only rename the label by myself such as Light and heavy\ngood vehicles (diesel) -\nGVX in order to get Light and heavy good vehicles (diesel) - GVX Right? -- View this message in context: http://r.789695.n4.nabble.com/label-wrap-gen-question-tp4639364p4639539.htm

[R] GEV distribution fitted by L-moment graph

2012-08-08 Thread Al Ehan
Hi, I have been having difficulties in finding packages/ codes that simplify plotting of a GEV fitted to dataset (by L-moments) that would print out graph comprising dataset versus gumbel reduced variate n return period at the same. Anyone can help me on this? Thanks. [[alternative HTML v

Re: [R] r-forge down?

2012-08-08 Thread Barry Rowlingson
On Wed, Aug 8, 2012 at 12:40 AM, Duncan Murdoch wrote: > This is presumably a temporary condition. It has been up recently, but > it is summer, and some people take vacations, so I won't give an > expected uptime. There's a note on that web page now explaining there's maintenance going on. ht

Re: [R] test if elements of a character vector contain letters

2012-08-08 Thread Liviu Andronic
On Tue, Aug 7, 2012 at 10:26 PM, Marc Schwartz wrote: > since there are alpha-numerics present, whereas the first option will: > >> grepl("[^[:alnum:]]", "ab%") > [1] TRUE > > > So, use the first option. > And I should start reading more carefully. The above works fine for me. I ended up defining

Re: [R] RQuantLib: SET_VECTOR_ELT() can only be applied to a 'list', not a 'symbol'

2012-08-08 Thread Enrico Schumann
Am 08.08.2012 09:54, schrieb Cren: # Hi all, # trying to run the following example code # from 'RQuantLib' package... HullWhite <- list(term = 0.055, alpha = 0.03, sigma = 0.01, gridIntervals = 40) Price <- rep(as.double(100),24) Type <- rep(as.character("C"), 24) Date <- se

Re: [R] Problem with global variable building a package

2012-08-08 Thread Uwe Ligges
On 08.08.2012 08:53, Eva Prieto Castro wrote: Hello, Greg: That is one of the solutions I have tried; the problem is that I can't obtain a package in this way, because an glo-internal.R file is generated (glo is the name of the package) and when I try to create the package an error message

Re: [R] process evaluation packages (slightly off topic)

2012-08-08 Thread Petr PIKAL
Thanks Rui igraph seems to be more related to my problem, so I will try to play with it. Regards Petr > > Hello, > > There's also, apparently for simple needs, plotrix::gantt.chart > See an example in > http://addictedtor.free.fr/graphiques/sources/source_74.R > > To transfer graphs to time

Re: [R] RQuantLib: SET_VECTOR_ELT() can only be applied to a 'list', not a 'symbol'

2012-08-08 Thread Cren
Enrico Schumann-2 wrote > > I cannot reproduce this error. I get... > > > sessionInfo() > *R version 2.15.1* (2012-06-22) > # Thank you for testing, Enrico (Italian? ), # it seems an updating issue. # I am trying to update everything possible to the latest # version because of compatibility.

Re: [R] help to program my function

2012-08-08 Thread Petr PIKAL
Hi Maybe it is time for you to read some basic stuff like R intro. It seems to me that you expect R to behave like some other language you know but probably your expectation is wrong. See inline > > HI > > >i have a problem please help me to solve it: > http://r.789695.n4.nabble.com/file/n4

Re: [R] Overlapping a Plot with Dataframe

2012-08-08 Thread Jim Lemon
On 08/07/2012 10:57 PM, mhimanshu wrote: Hello Everyone, I am trying to overlap a plot with a data set in the form of a data frame. Its very easy to overlap the data using "points" function. But the only problem I am facing is "Standard deviation bar" on the plot. data<- data.frame( x = c(3

Re: [R] Overlapping a Plot with Dataframe

2012-08-08 Thread Himanshu
Thanks Jim, Yes its working now. Thanks, Himanshu On Wed, Aug 8, 2012 at 12:10 PM, Jim Lemon wrote: > On 08/07/2012 10:57 PM, mhimanshu wrote: > >> Hello Everyone, >> >> I am trying to overlap a plot with a data set in the form of a data frame. >> >> Its very easy to overlap the data using "po

Re: [R] Interpreting predictions of svm

2012-08-08 Thread Jessica Streicher
You should give us the data is what you should do :) Aside from that: you can only make probability predictions if you activated it when making the model. On 07.08.2012, at 17:23, Camomille wrote: > Hi, I have some difficulties in interpreting the prediction of a svm model > using the package

Re: [R] R versus SAS

2012-08-08 Thread Duncan Murdoch
On 12-08-07 10:46 PM, Kjetil Halvorsen wrote: I found this on CrossValidated: "A medical statistician once told me, that they use SAS because if they make mistakes due to software bugs and it comes to lawsuits, SAS will recompensate them. R comes without warranty." I suspect that statistician

Re: [R] Error with convUL (PBSmapping)

2012-08-08 Thread niandra
Thank you. Il giorno 07/ago/2012, alle ore 20.36, R. Michael Weylandt ha scritto: > This looks like a fairly low level problem that you might need the > maintainer to look at. To get contact information, type > > maintainer("PBSmapping") > > at the R prompt. > > Best, > Michael > > On Tue, Au

[R] Pass Conditional & loop argument into a function

2012-08-08 Thread greatest.possible.newbie
Dear R community, I have been spending hours to solve my problem myself but I'm just unable to do it. So excuse me if I demand your time. My problem is, that inside my function, I want to pass conditional arguments (that are the result of a loop) into a function. I am sure that my approach is kind

Re: [R] Outer product from matrix by row and a vector

2012-08-08 Thread David Winsemius
On Aug 7, 2012, at 11:26 PM, Petr Savicky wrote: On Wed, Aug 08, 2012 at 06:03:19AM +, Ingo Reinhold wrote: Hi Jeff, David, what I'm trying to do is speed the currently used nest for loop a<-matrix(c(1:6),ncol=3) b<-c(1,2,3,4) result<-matrix(rep(0, times=dim(a)[1]*length(b)),nrow=dim(a)

[R] Cannot create Java virtual machine (-1)

2012-08-08 Thread Bos, Roger
All, I am trying to get the XLConnect package to work on one of my servers and I am getting an error that XLConnectJars cannot create the java virtual machine. The server is a 32-bit Windows Server 2003 machine with R 2.15.1. The version of Java installed on the server is J2SE Environment 5.

Re: [R] Repeated Aggregation with data.table

2012-08-08 Thread David Winsemius
On Aug 7, 2012, at 9:28 PM, arun wrote: HI, Try this: fun1<-function(x,.expr){ .expr<-expression(list(mean.z=mean(z),sd.z=sd(z))) z1<-eval(.expr) } #or fun1<-function(x,.expr){ .expr<-expression(list(mean.z=mean(z),sd.z=sd(z))) z1<-.expr } dat[,eval(z1),list(x)] dat[,eval(z1),list

[R] Hearsay -- was R vs SAS

2012-08-08 Thread John C Nash
I would ask contributors to this list to at least note that the sort of FUD (Fear Uncertainty and Doubt) of the type mentioned by the poster is hearsay. Do we have any contractual promises from SAS Inc. or other companies that they will guarantee their software? Have there been documented compen

Re: [R] Pass Conditional & loop argument into a function

2012-08-08 Thread S Ellison
First thing; are you trying to fit a model specified as y ~ X + X^2 + X^3 ? ... because if you are you're unlikely to get anything useful. That uses formula syntax in which ^ does not have the arithmetic power meaning; see section 11.1 'Defining statistical models; formulae' in 'an introduction

Re: [R] lm with a single X and step with several Xi-s, beta coef. quite different:

2012-08-08 Thread Ista Zahn
Hi, Sounds like suppression -- see e.g., http://www.jstor.org/stable/2988294?seq=1 for a discussion. Since this is not an R question but a statistical one, it may be more appropriate to post this question to a statistics forum such as http://stats.stackexchange.com/ Best, Ista On Tue, Aug 7, 20

[R] Coloring specific cells

2012-08-08 Thread namit
Hi , I have a data frame XYZ. i want to color the cells where ever we find "G" it should be green,"Y" yellow,"R" redetc and i want to send it from R to ppt.can any one help me on this. My data frame looks as below. XYZ: GY RTG Y RT AA 0

[R] dimnames in array

2012-08-08 Thread aleksandr russell
Hello, I'm working with an array; I'm trying to make it so that an array of dim(42,2,2) has names whose length corresponds to that of the array, and am hoping someone with experience with this can see what I'm not doing correctly: data11 = array(0,c(41,2,2)) y = lsoda(x0,times,fhn$fn.ode,pars

Re: [R] Can any one help me on this Issue

2012-08-08 Thread namit
Thanks Arunthis problem is resolved -- View this message in context: http://r.789695.n4.nabble.com/Can-any-one-help-me-on-this-Issue-tp4638664p4639579.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing l

[R] map axis on projected shapefiles

2012-08-08 Thread penguins
Hi, I have overlayed 2 projected shapefiles using the "plot" function. When I plot them the numbers next to the x-axis ticks are printed in scientific format (e.g. -4e+05, -3e+05,...,), and the problem is these dont match up to the projected point locations I also wish to overlay from my dataset

Re: [R] ERROR : cannot allocate vector of size (in MB & GB)

2012-08-08 Thread Rantony
Hi Arun, But I am using Windows(XP). From: arun kirshna [via R] [mailto:ml-node+s789695n4639435...@n4.nabble.com] Sent: Tuesday, August 07, 2012 10:49 PM To: Akkara, Antony (GE Energy, Non-GE) Subject: Re: ERROR : cannot allocate vector of size (in MB & GB) HI, If you are using linux

[R] sweave problem with special danish characters

2012-08-08 Thread Jonas Hal
Dear all, I use Sweave to create my reports. my pdf output contains only NA's whenever my text contains special Danish characters like ø, ø and å. I have upgraded R from 2.13 to 2.15 and the scripts worked fine under R 2.13. Sweave("Test.Rnw", encoding = "utf8") texi2dvi("Test.tex", pdf = TRUE)

Re: [R] Overlapping a Plot with Dataframe

2012-08-08 Thread mhimanshu
Hello John, in simple term, I have a Plot as an Output. Now I want to overlap the plot with a Dataframe having error bar. -- View this message in context: http://r.789695.n4.nabble.com/Overlapping-a-Plot-with-Dataframe-tp4639396p4639555.html Sent from the R help mailing list archive at Nabble

[R] how to call R commands from .net

2012-08-08 Thread agent dunham
Dear community, I want to call R from .net. Precisely I´d really like to call rpart from .net, plot decision trees and predict with them. I have installed the COM server (http://cran.es.r-project.org/contrib/extra/dcom/). Just to note, because it is from 2008, the only way to work with this ser

[R] party with the mob - fitted model lines in the terminal nodes

2012-08-08 Thread tudor
Hi useRs, I built a logistic-regression-based tree using mob (from package party) and am wondering if anyone could pinpoint to me where I can find the documentation for how the fitted model lines in the terminal nodes are computed. I am partitioning a data set based on a logistic model with inter

[R] How can we compare two vectors?

2012-08-08 Thread Sri krishna Devarayalu Balanagu
x=c(5, 8, 28, 29, 30) y=c(5, 8, 28, 29, 31) How can we compare these two vectors, whether each element is exactly matched with the elements in the other vector ? How can we get the non matched elements from both the vectors? Can anyone help? Notice: The information

Re: [R] Save & reload list objects

2012-08-08 Thread robgriffin247
So it was simple, or so it seemed. Now to throw a spanner in the works! I got this error when testing it on my more powerful desktop (using one of my lists put in to a new object) : #put list in to new object called "temp" temp<-rg.lmer #save save(temp,file="C:\\Users\\...\\R & Data\\rg.lmer.txt

[R] time series, uneven length

2012-08-08 Thread joel.green
I have 4 univariate time series that I believe have correlation between them, I want to create a VAR model between them all. However I have an issue as 3 of them are the same length, however the 4th is smaller. meaning that i cannot use the 4th variable , is there anyway R can get round this is

Re: [R] GEV distribution fitted by L-moment graph

2012-08-08 Thread J. R. M. Hosking
On 2012-08-08 04:22, Al Ehan wrote: Hi, I have been having difficulties in finding packages/ codes that simplify plotting of a GEV fitted to dataset (by L-moments) that would print out graph comprising dataset versus gumbel reduced variate n return period at the same. Anyone can help me on this?

Re: [R] Repeated Aggregation with data.table

2012-08-08 Thread arun
HI David, Thanks for testing. It's a bit strange.  Yesterday, the function was working perfectly.  Today, it is not working in my system.  Not sure what happened. A.K.  - Original Message - From: David Winsemius To: arun Cc: Elliot Joel Bernstein ; R help Sent: Wednesday, August

[R] help, please! matrix operations inside 3 nested loops

2012-08-08 Thread Fridolin
hello, this is my script: #1) read in data: daten<-read.table('K:/Analysen/STRUCTURE/input_STRUCTURE_tab_excl_5_282_559.txt', header=TRUE, sep="\t") daten<-as.matrix(daten) #2) create empty matrix: indxind<-matrix(nrow=617, ncol=617) indxind[1:20,1:19] #3) compare cells to each other, score: fo

Re: [R] Coloring specific cells

2012-08-08 Thread Jessica Streicher
If its just a few tables, can't you just copy the text and color it in PowerPoint? If you'd like Latex output i could help maybe, but i'm not sure what you expect. On 08.08.2012, at 14:19, namit wrote: > Hi , > > > I have a data frame XYZ. i want to color the cells where ever we find "G" it

Re: [R] How can we compare two vectors?

2012-08-08 Thread Michael Weylandt
Either == or %in% if you care about order or not. Michael On Aug 8, 2012, at 7:00 AM, Sri krishna Devarayalu Balanagu wrote: > x=c(5, 8, 28, 29, 30) > y=c(5, 8, 28, 29, 31) > How can we compare these two vectors, whether each element is exactly matched > with the elements in the other vector

Re: [R] How can we compare two vectors?

2012-08-08 Thread Jessica Streicher
> x=c(5, 8, 28, 29, 30) > y=c(5, 8, 28, 29, 31) > x==y [1] TRUE TRUE TRUE TRUE FALSE > which(x==y) [1] 1 2 3 4 > which(x!=y) [1] 5 > x[x!=y] [1] 30 > y[x!=y] [1] 31 RTFM! - just kidding ;) also, this might not work for floating point numbers, theres a section in the documentation/fa

Re: [R] sweave problem with special danish characters

2012-08-08 Thread Duncan Murdoch
On 08/08/2012 6:15 AM, Jonas Hal wrote: Dear all, I use Sweave to create my reports. my pdf output contains only NA's whenever my text contains special Danish characters like ø, ø and å. I have upgraded R from 2.13 to 2.15 and the scripts worked fine under R 2.13. Sweave("Test.Rnw", encoding

Re: [R] How can we compare two vectors?

2012-08-08 Thread ramakanth reddy
Hi krishna, 1.How can we compare these two vectors,whether each element is exactly matched with the elements in the other vector? x=c(5, 8, 28, 29, 30) y=c(5, 8, 28, 29, 31) z<-(x %in% y) > z [1] TRUE TRUE TRUE TRUE FALSE (OR) z<-identical(x,y) > z [1] FALSE (OR) to get values that are ma

Re: [R] sapply and matrix command

2012-08-08 Thread alijk1989 [via R]
Hi, I have made some progress speeding up my code. This is what I have at the moment: M1=sum(sapply(1:m, function(k){sum(sapply(1:m,function(j){w[k]*w[j]*LGD^2 (pmnorm(c(qnorm(Q[k]),qnorm(Q[j])),c(0,0),equicorr(2,rho[k,j]))-Q[k]*Q[j])}))})) I tried setting up a function as so: f1 <- function

Re: [R] map axis on projected shapefiles

2012-08-08 Thread penguins
Since posting this I've figured out the problem was due to the shapefiles I was reading in. As they were already projected they were in meters, as oppose to decimal degrees. By opening the shapefiles in ArcGIS, reprojecting them to WGS1984 and exporting, I could then reload them in R in decimal deg

Re: [R] Overlapping a Plot with Dataframe

2012-08-08 Thread John Kane
John Kane Kingston ON Canada > -Original Message- > From: bioinfo.himan...@gmail.com > Sent: Wed, 8 Aug 2012 02:12:00 -0700 (PDT) > To: r-help@r-project.org > Subject: Re: [R] Overlapping a Plot with Dataframe > > Hello John, > > in simple term, I have a Plot as an Output. > Now I wan

Re: [R] sweave problem with special danish characters

2012-08-08 Thread Ista Zahn
Hi Johas, I've run into a similar problem due to Emacs on OSX setting my locale to "C". If you're on Mac or Linux you can try setting Sys.setlocale(category = "LC_ALL", locale = "UTF-8") before calling Sweave() Best, Ista This doesn't æ ø å and pdf document only contains NA On Wed, Aug 8, 20

Re: [R] dimnames in array

2012-08-08 Thread David Winsemius
On Aug 8, 2012, at 3:19 AM, aleksandr russell wrote: Hello, I'm working with an array; I'm trying to make it so that an array of dim(42,2,2) has names whose length corresponds to that of the array, and am hoping someone with experience with this can see what I'm not doing correctly: data11

Re: [R] sweave problem with special danish characters

2012-08-08 Thread Jonas Hal
I am using win xp, R 2.15.0 and texi2dvi from the tools package. Regards Jonas -Oprindelig meddelelse- Fra: Duncan Murdoch [mailto:murdoch.dun...@gmail.com] Sendt: 8. august 2012 16:43 Til: Jonas Hal Cc: r-help@r-project.org Emne: Re: [R] sweave problem with special danish characters On

Re: [R] Save & reload list objects

2012-08-08 Thread Ista Zahn
Hi Rob, On Wed, Aug 8, 2012 at 8:50 AM, robgriffin247 wrote: > So it was simple, or so it seemed. > Now to throw a spanner in the works! > > I got this error when testing it on my more powerful desktop (using one of > my lists put in to a new object) : > > #put list in to new object called "temp"

Re: [R] sweave problem with special danish characters

2012-08-08 Thread Duncan Murdoch
On 08/08/2012 10:58 AM, Jonas Hal wrote: I am using win xp, R 2.15.0 and texi2dvi from the tools package. I'd suggest checking the encoding of your file (if I use latin1, I get an NA like you did). If that's not it, then you should upgrade to the latest release, or to the current patched ver

Re: [R] sapply and matrix command

2012-08-08 Thread R. Michael Weylandt
On Wed, Aug 8, 2012 at 9:17 AM, alijk1989 [via R] wrote: > > > Hi, > > I have made some progress speeding up my code. This is what I have at the > moment: > > > M1=sum(sapply(1:m, function(k){sum(sapply(1:m,function(j){w[k]*w[j]*LGD^2 > (pmnorm(c(qnorm(Q[k]),qnorm(Q[j])),c(0,0),equicorr(2,rho[k,j]

Re: [R] time series, uneven length

2012-08-08 Thread R. Michael Weylandt
On Wed, Aug 8, 2012 at 7:53 AM, joel.green wrote: > I have 4 univariate time series that I believe have correlation between them, > I want to create a VAR model between them all. > > However I have an issue as 3 of them are the same length, however the 4th is > smaller. meaning that i cannot use

Re: [R] help, please! matrix operations inside 3 nested loops

2012-08-08 Thread R. Michael Weylandt
On Wed, Aug 8, 2012 at 9:06 AM, Fridolin wrote: > hello, this is my script: > > #1) read in data: > daten<-read.table('K:/Analysen/STRUCTURE/input_STRUCTURE_tab_excl_5_282_559.txt', > header=TRUE, sep="\t") > daten<-as.matrix(daten) > > #2) create empty matrix: > indxind<-matrix(nrow=617, ncol=617

Re: [R] help to program my function

2012-08-08 Thread hafida
HI peter there is the function that i want to programm (joint in pdf folder). my data is > dataexpv Ti1 265.792 26 1579.523 26 2323.704 28 68.855 28 426.076 28 110.297 28 108.298 28 1067.609 30 17.0510 30 22.6611 30 21.0212 30 175.8813 30 139.0714 30 144.1215

Re: [R] help, please! matrix operations inside 3 nested loops

2012-08-08 Thread Berend Hasselman
Fridolin wrote > > hello, this is my script: > > #1) read in data: > daten<-read.table('K:/Analysen/STRUCTURE/input_STRUCTURE_tab_excl_5_282_559.txt', > header=TRUE, sep="\t") > daten<-as.matrix(daten) > > #2) create empty matrix: > indxind<-matrix(nrow=617, ncol=617) > indxind[1:20,1:19] >

[R] Odp: help, please! matrix operations inside 3 nested loops

2012-08-08 Thread Petr PIKAL
Hi > > hello, this is my script: > > #1) read in data: > daten<-read.table('K:/Analysen/STRUCTURE/input_STRUCTURE_tab_excl_5_282_559.txt', > header=TRUE, sep="\t") > daten<-as.matrix(daten) If there is any column with nonnumeric values it will transfer all numeric values from daten data.frame

Re: [R] help to program my function

2012-08-08 Thread Petr PIKAL
Hi please cc your post also to Rhelp list, others may give you better/quicker answer. > > > HI peter > there is the function that i want to programm (joint in pdf folder). No pdf allowed. > my data is > dataexpv Ti1 265.792 26 1579.523 26 2323.704 > 28 68.855 28 426.07

Re: [R] dimnames in array

2012-08-08 Thread aleksandr russell
Hi again, Thanks for your suggestion. I think I can clarify the thing by looking at the bare bones: data11 = array(0,c(41,2,2)) m8<-cbind(.29,1:41) m9<-as.array(m8,dim=c(41,2)) data11[,1,]=m9 data11[,2,]=m9 varnames=c("V","R") colnames(dat

[R] R help

2012-08-08 Thread chong hui qi
Hi I'm trying to increase the memory limit but Im facing this problem "Error: unexpected symbol in ""C:\\Program Files\\R\\R-2.15.0\\bin\\i386\\Rgui.exe" --max-mem-size=500M" " Hope you can help me out Thank You [[alternative HTML version deleted]] __

Re: [R] sapply and matrix command

2012-08-08 Thread alijk1989 [via R]
Hi Michael, Thanks for your response. Here is a simple example of what I am trying to do: w=rep(0.02,10) Q=rep(0.02,10) rho=matrix(0.5,nrow=10,ncol=10) m=10 LGD=0.45 M1=sum(sapply(1:m, function(k){sum(sapply(1:m,function(j){w[k]*w[j]*LGD^2*(pmnorm(c(qnorm(Q[k]),qnorm(Q[j])),c(0,0),equicorr(2,r

[R] ggplot2 with separate average lines

2012-08-08 Thread Alan Miller
Hi I'm just starting off with R but is this correct? I have this data set. > data Type ID Size 1 Reqd 244808 1024 2 Reqd 244045 512 3 Reqd 245427 800 4 Reqd 245423 1024 5 Reqd 244983 1024 6 Used 244808 615 7 Used 244045 33 8 Used 245427 261 9 Used 245423 461 10 Used 244983 1

Re: [R] Confidence bands around LOESS

2012-08-08 Thread John Kane
I may have missed something entirely but I think ggplot2 will do this for you pretty well automatically library(ggplot2) mydata <- read.csv("http://doylesdartden.com/smoothing.csv";, sep=",") p <- ggplot(mydata , aes( X_Axis_Parameter, Y_Axis_Parameter )) + geom_point() + geom_

Re: [R] Setting Number of Displayed Digits

2012-08-08 Thread Rich Shepard
On Tue, 7 Aug 2012, arun wrote: formatC() may work for you. Thanks for the lesson. Rich __ 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

Re: [R] dimnames in array

2012-08-08 Thread David Winsemius
On Aug 8, 2012, at 8:43 AM, aleksandr russell wrote: Hi again, Thanks for your suggestion. I think I can clarify the thing by looking at the bare bones: data11 = array(0,c(41,2,2)) m8<-cbind(.29,1:41) m9<-as.array(m8,dim=c(41,2)) data11[,1,]=m9 data11[

Re: [R] Confidence bands around LOESS

2012-08-08 Thread Bert Gunter
Please Note: loess() is a fitting algorithm, so no parameters are estimated and no _confidence_ intervals nor bands can be produced. What predict.loess() _does_ give is standard errors for the fits/predictions, and you can add and subtract as many of these as you like to produce standard error b

Re: [R] ggplot2 with separate average lines

2012-08-08 Thread John Kane
Looks like it works :) I'd suggest not using attach(). It can cause problems in some cases. gglot2 is set up to read the relevant data frame so p <- ggplot( data , aes( Whatever. )) + and so on However "data" is a defined function in R so something like mydata or df1 is a better choice.

Re: [R] R help

2012-08-08 Thread Jeff Newmiller
There seem to be too many quotes in your error message... the path and name of the program should be quoted, the argument should not (or should be separately quoted). 500M isn't especially large these days... be sure to check your existing (default) value using the memory.limit() function befor

Re: [R] Pass Conditional & loop argument into a function

2012-08-08 Thread greatest.possible.newbie
Thank you for your improvement suggestions. I forgot to write the I() and it's clear that I have to specify the data (which I was indicating with the '...' without explicitly writing it). In this special case where I need a formula you really helped me with the function as.formula() to transform t

Re: [R] Pass Conditional & loop argument into a function

2012-08-08 Thread Rui Barradas
Hello, The other problem with the boxplot question is that it's ill defined. Those data[,1], ..., data[, n] variables, what are they? All of them numeric? If so, then where is the factor? Some time ago there was a question on multiple boxplots with cut() defining the factor. I wrote a benera

Re: [R] Pass Conditional & loop argument into a function

2012-08-08 Thread Jeff Newmiller
The quotes shown on a string when you print it are artifacts of making the string usable in source code... they are not present in the string as it is stored in memory. (Same goes for escape sequences, such that "\n" occupies one byte of UTF8 storage.) You can use the cat function to send the s

Re: [R] Interpreting predictions of svm

2012-08-08 Thread Camille Hénaff
Jessica- Thanks for responding quickly to me; I have fund a solution to my problem this morning. Kind regards, Camille [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEA

Re: [R] label_wrap_gen question

2012-08-08 Thread Brian Diggs
On 8/8/2012 12:30 AM, vd3000 wrote: Hi, Brian So, I could only rename the label by myself such as Light and heavy\ngood vehicles (diesel) -\nGVX in order to get Light and heavy good vehicles (diesel) - GVX Right? First an etiquette comment: it is customary to quote context when replying t

[R] getting started with R and mysql

2012-08-08 Thread Michael Friendly
I have a mysql database I installed on my Ubuntu server, and I'm trying to figure out how to access and analyze it from a Win XP machine running R 2.15.1. I thought that RMySQL was the way to go. (Is there an easier way?) I tried to follow the instructions at http://www.r-bloggers.com/installi

Re: [R] getting started with R and mysql

2012-08-08 Thread Curt Seeliger
> I have a mysql database I installed on my Ubuntu server, and I'm trying > to figure out how to access and analyze it > from a Win XP machine running R 2.15.1. > > I thought that RMySQL was the way to go. (Is there an easier way?) We use RODBC. Is there a reason that won't work? cur -- Curt

[R] Advice: How to best ensure column values match in different vectors?

2012-08-08 Thread DG Christensen
Hello all, I would like some advice on how to order elements in a vector. Background: my company is running a k-means clustering model on our historical data warehouse of products, which will produce a matrix of cluster centers. Then, on our production web servers, we will take newly created prod

[R] Installing HDF5 package

2012-08-08 Thread sihemsaad
Dear all I am a new user of R and I have been trying to install the HDF5 package but I don't seem to find any easily. I tried looking in the archives but I didn't succeed in installing the packages there. I guess they must be for Unix If any body succeeded in installing this library, please give

Re: [R] help to program my function

2012-08-08 Thread hafida
hi peter the pdf folder is in http://r.789695.n4.nabble.com/file/n4639434/aj.pdf -- View this message in context: http://r.789695.n4.nabble.com/help-to-program-my-function-tp4639434p4639629.html Sent from the R help mailing list archive at Nabble.com. _

Re: [R] Cannot create Java virtual machine (-1)

2012-08-08 Thread Martin Studer
Hi Roger, note that since version 0.1-7 of XLConnect Java 1.6 is required - it won't work with Java 1.5. However, the problem you describe is another one. It's most likely related to your Java setup/configuration and therefore independent of XLConnect. You can verify this by running the following

[R] Wilcoxon test

2012-08-08 Thread BELLAY Juliette
Dear list, I am facing a problem in my statistical analyses on R. My experiments are about plants, I record there growth after each cutting (every 3 weeks). 'BC' is for the plant, and '1' to '5' is the time of cutting and recording. The data and R script are : "" BourdCoup <- c(21, 7.2, 9.2, 0,

Re: [R] summing and combining rows

2012-08-08 Thread arun
HI, From the ?aggregate(), formula: a formula, such as ‘y ~ x’ or ‘cbind(y1, y2) ~ x1 + x2’,   where the ‘y’ variables are numeric data to be split into   groups according to the grouping ‘x’ variables (usually   factors). So, I converted your data to factors for the grou

Re: [R] summing and combining rows

2012-08-08 Thread arun
Hello, I tried with ddply  ddply(dat1,.(Data,Plot,Lat,LatCat,Elevation,Type,Area,Density),summarize,sum(Stems))   Data  Plot  Lat LatCat Elevation    Type Area  Density ..1 1  VTM 39C16 39.76282   Lat6  1500 Conifer 751.5347  0.0   0 2  VTM 39C16 39.76282   Lat6  1500 Conifer

[R] mgcv and gamm4: REML, GCV, and AIC

2012-08-08 Thread john benson
Hi, I've been using gamm4 to build GAMMs for exploring environmental influences on genetic ancestry. Things have gone well and I have 2 very straightforward questions: 1. I've used method=REML. Am I correct that this is an alternative method for estimating the smooth functions in GAMMs rather

[R] Showing multiple Probability Plots on the same graph using the plot(cenros) command from the NADA library

2012-08-08 Thread plebel
I am currently conducting a performance analysis for an urban stormwater management facility which employs Low Impact Development practices. I have the Event Mean Concentrations (EMCs) of several water quality constituents (ex: Total Phosphorus (TP), Total Oxidized Nitrogen (OxN), etc) at inflow an

[R] decimal points midline

2012-08-08 Thread array chip
Hi, does anyone know how to make decimal points midline when plotting? The journal to which we are going to submit a manuscript require this particular formatting, and it gives direction how to do this on PC: hold down ALT key and type 0183 on the number pad Thanks John [[alternative H

Re: [R] Installing HDF5 package

2012-08-08 Thread John Kane
You have not told us what you have done or what your OS is. In any case R is case sensitive so try: install.packages("hdf5") John Kane Kingston ON Canada > -Original Message- > From: saadsihe...@gmail.com > Sent: Wed, 8 Aug 2012 09:25:06 -0700 (PDT) > To: r-help@r-project.org > Subject

Re: [R] Installing HDF5 package

2012-08-08 Thread Uwe Ligges
On 08.08.2012 20:54, John Kane wrote: You have not told us what you have done or what your OS is. In any case R is case sensitive so try: install.packages("hdf5") Which won't work on Windows since binaries are not available under that platform. Another package brandnew on CRAN is h5r. Bi

Re: [R] Installing HDF5 package

2012-08-08 Thread Martin Morgan
On 08/08/2012 12:10 PM, Uwe Ligges wrote: On 08.08.2012 20:54, John Kane wrote: You have not told us what you have done or what your OS is. In any case R is case sensitive so try: install.packages("hdf5") Which won't work on Windows since binaries are not available under that platform. An

Re: [R] decimal points midline

2012-08-08 Thread John Kane
Can you point to an example? It sounds like the journal is still using typewriters. John Kane Kingston ON Canada > -Original Message- > From: arrayprof...@yahoo.com > Sent: Wed, 8 Aug 2012 11:52:27 -0700 (PDT) > To: r-help@r-project.org > Subject: [R] decimal points midline > > Hi, do

Re: [R] decimal points midline

2012-08-08 Thread Uwe Ligges
On 08.08.2012 21:22, John Kane wrote: Can you point to an example? It sounds like the journal is still using typewriters. Anyway, setting options(OutDec="ALT0183") and replace ALT0183 by the keystrokes your described below or by options(OutDec="\xB7") Best, Uwe Ligges John Kane Ki

Re: [R] Alternating between "for loops"

2012-08-08 Thread Claudia Penaloza
Answers inserted in BLUE below On Thu, Aug 2, 2012 at 5:34 PM, Claudia Penaloza wrote: > Thank you very much for all your suggestions. I am very sorry my code is > so crude (it gives me a headache too!), I'm very new to R programing. I > will have to look at your suggestions/questions very carefu

Re: [R] NADA Package: Referencing Data Frame Columns

2012-08-08 Thread MacQueen, Don
Hi Rich, I may not have the complete picture here, but I do see what looks to me like a problem with your chem.cast. Specifically, since it has only a single detection indicator column (ceneq1), it implies that within any single sample either all the analytes were detected, or all were not. Not w

Re: [R] decimal points midline

2012-08-08 Thread array chip
Thanks Uwe, options(OutDec="\xB7") worked! From: Uwe Ligges To: John Kane roject.org> Sent: Wednesday, August 8, 2012 12:43 PM Subject: Re: [R] decimal points midline On 08.08.2012 21:22, John Kane wrote: > Can you point to an example?  It sounds like the

  1   2   >