Re: [R] Logical vectors

On Wed, Nov 3, 2010 at 10:50 PM, Stephen Liu wrote: > Hi folks, > > Pls help me to understand follow; > > An Introduction to R > > 2.4 Logical vectors > http://cran.r-project.org/doc/manuals/R-intro.html#R-and-statistics > > 1) >> x > [1] 1 2 3 4 5 a vector, x, is defined with 5 elements, {1, 2,

Re: [R] Logical vectors

On Wed, 3 Nov 2010, Stephen Liu wrote: [snip] 2) x [1] 1 2 3 4 5 temp <- x > 1 temp [1] FALSE TRUE TRUE TRUE TRUE Why NOT temp [1] TRUE FALSE FALSE FALSE FALSE ? Maybe because of the definition of ">" (greater (!) than)? Or do you expect 1 to be greater than 1 and not greate

[R] Logical vectors

Hi folks, Pls help me to understand follow; An Introduction to R 2.4 Logical vectors http://cran.r-project.org/doc/manuals/R-intro.html#R-and-statistics 1) > x [1] 1 2 3 4 5 > temp <- x != 1 > temp [1] FALSE TRUE TRUE TRUE TRUE > 2) > x [1] 1 2 3 4 5 > temp <- x > 1 > temp [1] FALSE TRU

Re: [R] avoiding too many loops - reshaping data

Beware of facile comparisons of this sort -- they may be apples and nematodes. I cannot speak to the others, but (1) tapply does not yield a data frame and (2) tapply actually **is** a (efficient, disguised) loop (at the interpreter level, essentially). I suspect what makes it so much faster is th

[R] Best Fit line trouble with rsruby

Hello, I am using R, through rsruby, to create a graph and best fit line for a set of data points, regarding data collected in a Chemistry class. The problem is that although the graph functions perfectly properly, the best fit line will not work. I initially used code I pretty much copied from a

[R] cross-validation for choosing regression trees

Dear All, We came across a problem when using the "tree" package to analyze our data set. First, in the "tree" function, if we use the default value "mindev=0.01", the resulting regression tree has a single node. So, we set "mindev=0", and obtain a tree with 931 terminal nodes. However, when we

Re: [R] Auto-killing processes spawned by foreach::doMC

Hi, I am also interest in ways to in R send signals to other R sessions/processes, ideally in (what appears to be) an OS-independent way. For what it is worth, related question have been asked before, cf. R-devel thread 'Sending signals to current R process from R running under MS Windows (c.f. E

[R] Problems with mgui

Hello I am using the mgui function in the following way: mgui ( graf_cuenta_margen_interfaz,title=c("Gráficas","Histogramas valoración (No lineal) Cuenta de Margen"),exec="Graficar",argText=list(fecha_adelante="Fecha adelante"),closeOnExec=TRUE,output=NULL,,helps=list(fecha_adelante=paste("La valo

Re: [R] multivariate Poisson distribution

On 03-Nov-10 21:06:53, Jourdan Gold wrote: > Hello, from a search of the archives and functions, I am looking for > information on creating random correlated counts from a multivariate > Poisson distribution. I can not seem to find a function that does this. > Perhaps, it has not yet been created.

Re: [R] Rd installation (not markup language) primer?

On 03/11/2010 6:28 PM, ivo welch wrote: I have a set of functions that I always load on startup. for example, there is my now infamous "is.defined()" function. I would like to add some documentation for these functions, so that I can do a ?is.defined inside R. The documentation tells me

Re: [R] multivariate Poisson distribution

Jourdan Gold uoguelph.ca> writes: > > > Hello, from a search of the archives and functions, > I am looking for information on creating random > correlated counts from a multivariate Poisson distribution.  > I can not seem to find a function that > does this. Perhaps, it has not yet  been c

[R] Rd installation (not markup language) primer?

I have a set of functions that I always load on startup. for example, there is my now infamous "is.defined()" function. I would like to add some documentation for these functions, so that I can do a ?is.defined inside R. The documentation tells me how to mark up Rd files is very good, but I

Re: [R] avoiding too many loops - reshaping data

Here is the summary of methods. tapply is the fastest! library(reshape) system.time(for(i in 1:1000)cast(melt(mydf, measure.vars = "value"), city ~ brand,fun.aggregate = sum)) user system elapsed 18.400.00 18.44 library(reshape2) system.time(for(i in 1:1000)dcast(mydf,city ~ brand, su

[R] how to handle 'g...@gtdata' ?

I have a few questions about GenABEL, gwaa data. 1) is there a universal way that most GenABEL people use to add more individuals into a 'gwaa' data? For example, I have a 'gwaa' data, but I need to add some dummy parents, for 'g...@phdata', it's easy to add these rows, but for 'g...@gtdata', I t

[R] Auto-killing processes spawned by foreach::doMC

Hi all, Sometimes I'll find myself "ctrl-c"-ing like a madman to kill some code that's parallelized via foreach/doMC when I realized that I just set my cpu off to do something boneheaded, and it will keep doing that thing for a while. In these situations, since I interrupted its normal execution,

Re: [R] biding rows while merging at the same time

On Nov 3, 2010, at 5:38 PM, Dimitri Liakhovitski wrote: Never mind - I found it in "reshape" package: rbind.fill I wonder if it's still in reshape2. Look in plyr. -- David. Dimitri __ R-help@r-project.org mailing list https://stat.ethz.ch/mailma

Re: [R] biding rows while merging at the same time

Just merge(df1, df2, all = TRUE) does it, yes? Dimitri Liakhovitski wrote: Hello! I have 2 data frames like this (well, actually, I have 200 of them): df1<-data.frame(location=c("loc 1","loc 2","loc 3"),date=c("1/1/2010","1/1/2010","1/1/2010"), a=1:3,b=11:13,c=111:113) df2<-data.frame(locati

Re: [R] Loop

On Nov 3, 2010, at 5:03 PM, Matevž Pavlič wrote: Hi, Thanks for the help and the manuals. Will come very handy i am sure. But regarding the code i don't hink this is what i wantbasically i would like to repeat bellow code : w1<-table(lit$W1) w1<-as.data.frame(w1) It appears you are

Re: [R] biding rows while merging at the same time

Never mind - I found it in "reshape" package: rbind.fill I wonder if it's still in reshape2. Dimitri On Wed, Nov 3, 2010 at 5:34 PM, Dimitri Liakhovitski wrote: > Hello! > > I have 2 data frames like this (well, actually, I have 200 of them): > > df1<-data.frame(location=c("loc 1","loc 2","loc >

[R] biding rows while merging at the same time

Hello! I have 2 data frames like this (well, actually, I have 200 of them): df1<-data.frame(location=c("loc 1","loc 2","loc 3"),date=c("1/1/2010","1/1/2010","1/1/2010"), a=1:3,b=11:13,c=111:113) df2<-data.frame(location=c("loc 1","loc 2","loc 3"),date=c("2/1/2010","2/1/2010","2/1/2010"), a=4:6,c=

[R] multivariate Poisson distribution

Hello, from a search of the archives and functions, I am looking for information on creating random correlated counts from a multivariate Poisson distribution.  I can not seem to find a function that does this. Perhaps, it has not yet  been created. Has anyone created an R package that do

Re: [R] Loop

Hi, Thanks for the help and the manuals. Will come very handy i am sure. But regarding the code i don't hink this is what i wantbasically i would like to repeat bellow code : w1<-table(lit$W1) w1<-as.data.frame(w1) write.table(w1,file="w1.csv",sep=";",row.names=T, dec=".") w1<- w1[order(w

Re: [R] avoiding too many loops - reshaping data

Want to thank everyone once more for pointing in "reshape" direction. Saved me about 16 hours of looping! Dimitri On Wed, Nov 3, 2010 at 4:38 PM, Dimitri Liakhovitski wrote: > In reshape2 this does the job: > > dcast(mydf,city~brand,sum) > > > On Wed, Nov 3, 2010 at 4:37 PM, Dimitri Liakhovitski

[R] multi-level cox ph with time-dependent covariates

Dear all, I would like to know if it is possible to fit in R a Cox ph model with time-dependent covariates and to account for hierarchical effects at the same time. Additionally, I'd like also to know if it would be possible to perform any feature selection on this model fit. I have a data set th

Re: [R] avoiding too many loops - reshaping data

In reshape2 this does the job: dcast(mydf,city~brand,sum) On Wed, Nov 3, 2010 at 4:37 PM, Dimitri Liakhovitski wrote: > Thanks a lot! > > Yes - I just found the reshape package too - and guess what, my math was > wrong! > reshape2 seems like the more up-to-date version of reshape. > > Wonder w

Re: [R] avoiding too many loops - reshaping data

Thanks a lot! Yes - I just found the reshape package too - and guess what, my math was wrong! reshape2 seems like the more up-to-date version of reshape. Wonder what's faster - xtabs or dcast... Dimitri On Wed, Nov 3, 2010 at 4:32 PM, Henrique Dallazuanna wrote: > Try this: > >  xtabs(value ~ c

Re: [R] avoiding too many loops - reshaping data

Try this: xtabs(value ~ city + brand, mydf) On Wed, Nov 3, 2010 at 6:23 PM, Dimitri Liakhovitski < dimitri.liakhovit...@gmail.com> wrote: > Hello! > > I have a data frame like this one: > > > mydf<-data.frame(city=c("a","a","a","a","a","a","a","a","b","b","b","b","b","b","b","b"), > brand=c("x

Re: [R] avoiding too many loops - reshaping data

Hadley's reshape package (google for it) can do this. There's a nice intro on the site. > library(reshape) > cast(melt(mydf, measure.vars = "value"), city ~ brand, fun.aggregate = sum) city x y z 1a 3 23 450 2b 12 42 231 Although the numbers differ slightly? I've heard of

[R] avoiding too many loops - reshaping data

Hello! I have a data frame like this one: mydf<-data.frame(city=c("a","a","a","a","a","a","a","a","b","b","b","b","b","b","b","b"), brand=c("x","x","y","y","z","z","z","z","x","x","x","y","y","y","z","z"), value=c(1,2,11,12,111,112,113,114,3,4,5,13,14,15,115,116)) (mydf) What I need to get i

[R] Loop

Hi all, I managed to do what i want (with the great help of thi mailing list) manually . Now i would like to automate it. I would probably need a for loop for to help me with this...but of course I have no idea how to do that in R. Bellow is the code that i would like to be replicated fo

Re: [R] NFFT on a Zoo?

> From: ggrothendi...@gmail.com > Date: Wed, 3 Nov 2010 15:27:13 -0400 > To: flym...@gmail.com > CC: r-help@r-project.org; rpy-l...@lists.sourceforge.net > Subject: Re: [R] NFFT on a Zoo? > > On Wed, Nov 3, 2010 at 2:59 PM, Bob Cunningham wrote: > > I

Re: [R] programming questions

On Wed, Nov 3, 2010 at 1:04 PM, ivo welch wrote: > thanks, eric---I need a little more clarification.  *yes, I write > functions and then forget them.  so I want them to be self-sufficient. >  I want to write functions that check all their arguments for > validity.)  For example, > >  my.fn <- fun

Re: [R] FW: optim works on command-line but not inside a function

Well, the function should not be able to be infinite as IRT.llZetaLambdaCorrNan is a sum of products of either one or zero and log(1+exp(x)) or log(1+exp(-x)) (these logs are always bigger or equal to log(1)=0) Further more, I bounded x to be finite to fix my problem (as I expected that it might tr

Re: [R] NFFT on a Zoo?

On 11/03/2010 12:27 PM, Gabor Grothendieck wrote: On Wed, Nov 3, 2010 at 2:59 PM, Bob Cunningham wrote: I have an irregular time series in a Zoo object, and I've been unable to find any way to do an FFT on it. More precisely, I'd like to do an NFFT (non-equispaced / non-uniform time FFT) o

Re: [R] programming questions

thanks, eric---I need a little more clarification. *yes, I write functions and then forget them. so I want them to be self-sufficient. I want to write functions that check all their arguments for validity.) For example, my.fn <- function( mylist ) { stop.if.not( is.defined(mylist) ) #

Re: [R] optim works on command-line but not inside a function

Damokun wrote: > > Dear all, > > I am trying to optimize a logistic function using optim, inside the > following functions: > #Estimating a and b from thetas and outcomes by ML > > IRT.estimate.abFromThetaX <- function(t, X, inits, lw=c(-Inf,-Inf), > up=rep(Inf,2)){ > optRes <- optim(inits

Re: [R] programming questions

On Nov 3, 2010, at 3:32 PM, ivo welch wrote: thanks, barry and eric. I didn't do a good job---I did an awful job. is.defined(never.before.seen$anything) ## if a list does not exist, anything in it does not exist either Except the $ function return NULL rather than an error and you alrea

Re: [R] deleteing all but some observations in data.frame

Note that these methods don't 'delete' observations. They all create brand new objects that are subsets of the test.df object. You can effectively 'delete' the observations by replacing the original data.frame with the returned object... so test.df <- head(test.df, 20) Erik Iverson wrote: It

Re: [R] deleteing all but some observations in data.frame

It depends on which 20 you want. If you have a data.frame called 'test.df', you can do: #first 20 test.df[20, ] -or- head(test.df, 20) #random 20 test.df[sample(nrow(test.df), 20), ] None of this was tested, but it should be a start. --Erik Matevž Pavlič wrote: Hi, I am sure that can

[R] deleteing all but some observations in data.frame

Hi, I am sure that can be done in R How would i delete all but let say 20 observations in data.frame? Thank you, M [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listin

Re: [R] programming questions

alas, should R not come with an is.defined() function? ?exists a variable may never have been created, and this is different from a variable existing but holding a NULL. this can be the case in the global environment or in a data frame. > is.null(never.before.seen) Error: objected 'ne

Re: [R] smooth: differences between R and S-PLUS

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Nicola > Sturaro Sommacal (Quantide srl) > Sent: Wednesday, November 03, 2010 10:41 AM > To: r-help@r-project.org > Subject: [R] smooth: differences between R and S-PLUS > > Hi!

Re: [R] programming questions

thanks, barry and eric. I didn't do a good job---I did an awful job. alas, should R not come with an is.defined() function? a variable may never have been created, and this is different from a variable existing but holding a NULL. this can be the case in the global environment or in a data fram

Re: [R] NFFT on a Zoo?

On Wed, Nov 3, 2010 at 2:59 PM, Bob Cunningham wrote: > I have an irregular time series in a Zoo object, and I've been unable to > find any way to do an FFT on it.  More precisely, I'd like to do an NFFT > (non-equispaced / non-uniform time FFT) on the data. > > The data is timestamped samples fro

Re: [R] Lattice plots for images

On Wed, Nov 3, 2010 at 8:13 AM, Neba Funwi-Gabga wrote: > Hello UseRs, > I need help on how to plot several raster images (such as those obtained > from a kernel-smoothed intensity function) in a layout > such as that obtained from the lattice package. I would like to obtain > something such as ob

Re: [R] getting p-values from fitted ARIMA

On Wed, 3 Nov 2010, h0453...@wu.ac.at wrote: Hi I fitted an ARIMA model using the function arima(). The output consists of the fitted coefficients with their standard errors. However i need information about the significance of the coefficients, like p-values. I hope you can help me on that

[R] NFFT on a Zoo?

I have an irregular time series in a Zoo object, and I've been unable to find any way to do an FFT on it. More precisely, I'd like to do an NFFT (non-equispaced / non-uniform time FFT) on the data. The data is timestamped samples from a cheap self-logging accelerometer. The data is weakly re

Re: [R] getting p-values from fitted ARIMA

Hi Stefan, Take a look at https://stat.ethz.ch/pipermail/r-help/2009-June/202173.html HTH, Jorge On Wed, Nov 3, 2010 at 2:50 PM, <> wrote: > Hi > > I fitted an ARIMA model using the function arima(). The output consists of > the fitted coefficients with their standard errors. > > However i nee

Re: [R] '=' vs '<-'

@all: Does it seem reasonable to add a discussion of '=' vs. '<-' to the FAQ? It seems a regular question and something of a "hot" topic to debate. @KM Here are links I've accumulated to prior discussions on this topic. I am pretty certain they are all unique. http://blog.revolutionanalytics

[R] getting p-values from fitted ARIMA

Hi I fitted an ARIMA model using the function arima(). The output consists of the fitted coefficients with their standard errors. However i need information about the significance of the coefficients, like p-values. I hope you can help me on that issue... ciao Stefan

Re: [R] programming questions

On Wed, Nov 3, 2010 at 6:17 PM, ivo welch wrote: > yikes.  this is all my fault.  it was the first thing that I ever > defined when I started using R. > >   is.defined <- function(name) exists(as.character(substitute(name))) > > I presume there is something much better... You didn't do a good jo

Re: [R] '=' vs '<-'

On Wed, Nov 3, 2010 at 6:04 PM, km wrote: > Hi all, > > can we use '=' instead of '<-' operator for assignment in R programs? Yes, mostly, you can also use 'help' to ask such questions: > help("=") The operators ‘<-’ and ‘=’ assign into the environment in which they are evaluated. T

Re: [R] '=' vs '<-'

Yes, but <- is preferred. Note, there are also some differences. You can do the following: > a <- 10 > b = 10 > identical(a,b) [1] TRUE And you can also do > myFun <- function(x, y = 100){ + result <- x*y + result} > myFun(x = 20) [1] 2000 But, you cannot use '<-' to define the arguments of a

Re: [R] spliting first 10 words in a string

Hi all, Thanks for all the help. I realize i have a lot to learn in R but i love it. m From: steven mosher [mailto:mosherste...@gmail.com] Sent: Tuesday, November 02, 2010 11:45 PM To: Matevž Pavlič Cc: David Winsemius; Gaj Vidmar; r-h...@stat.math.ethz.ch Subject: Re: [R] spliting

Re: [R] programming questions

For data frames you can also use with() in your example: with(d, exists("z")) -- Jonathan P. Daily Technician - USGS Leetown Science Center 11649 Leetown Road Kearneysville WV, 25430 (304) 724-4480 "Is the room still a room when its empty? Does the room, the t

Re: [R] programming questions

yikes. this is all my fault. it was the first thing that I ever defined when I started using R. is.defined <- function(name) exists(as.character(substitute(name))) I presume there is something much better... /iaw On Wed, Nov 3, 2010 at 2:12 PM, Erik Iverson wrote: > > > ivo welch wrote:

Re: [R] programming questions

On 03/11/2010 2:05 PM, ivo welch wrote: quick programming questions. I want to "turn on" more errors. there are two traps I occasionally fall into. * I wonder why R thinks that a variable is always defined in a data frame. > is.defined(d) [1] FALSE > d= data.frame( x=1:5,

Re: [R] programming questions

ivo welch wrote: quick programming questions. I want to "turn on" more errors. there are two traps I occasionally fall into. * I wonder why R thinks that a variable is always defined in a data frame. > is.defined(d) [1] FALSE > d= data.frame( x=1:5, y=1:5 ) > is.defined(

[R] programming questions

quick programming questions. I want to "turn on" more errors. there are two traps I occasionally fall into. * I wonder why R thinks that a variable is always defined in a data frame. > is.defined(d) [1] FALSE > d= data.frame( x=1:5, y=1:5 ) > is.defined(d$z) [1] TRUE

[R] '=' vs '<-'

Hi all, can we use '=' instead of '<-' operator for assignment in R programs? regards, KM __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and

Re: [R] Orthogonalization with different inner products

Suppose one wanted to consider random variables X_1,...X_n and from each subtract off the piece which is correlated with the previous variables in the list. i.e. make new variables Z_i so that Z_1=X_1 and Z_i=X_i-cov(X_i,Z_1)Z_1/var(Z_1)-...- cov(X_i,Z__{i-1})Z__{i-1}/var(Z_{i-1}) I have code

[R] smooth: differences between R and S-PLUS

Hi! I am studying differences between R and S-PLUS smooth() functions. I know from the help that they worked differently, so I ask: - exist a package that permit to have the same results? - alternatively, someone know how can I obtain the same results in R, using a self made script? I know that

Re: [R] How to unquote string in R

lord12 wrote: s= "Hey" a = "Hello" table = rbind(s,a) write.table(table,paste("blah",".PROPERTIES",sep = ""),row.names = FALSE,col.names = FALSE) In my table, how do I output only the words and not the words with the quotations? You read the help page for the function you're using :). Fro

Re: [R] density() function: differences with S-PLUS

Dear William, thank you very much for your reply. I see it only after my reply to Joshua. Unfortunately I cannot try until tomorrow, because I don't have S-PLUS on this machine. Thanks again. Nicola 2010/11/3 William Dunlap > Did you get my reply (1:31pm PST Tuesday) > to your request? It

Re: [R] How to unquote string in R

Put the quote = FALSE argument in write.table On Wed, Nov 3, 2010 at 2:13 PM, lord12 wrote: > > > s= "Hey" > a = "Hello" > table = rbind(s,a) > write.table(table,paste("blah",".PROPERTIES",sep = ""),row.names = > FALSE,col.names = FALSE) > > In my table, how do I output only the words and not th

[R] How to unquote string in R

s= "Hey" a = "Hello" table = rbind(s,a) write.table(table,paste("blah",".PROPERTIES",sep = ""),row.names = FALSE,col.names = FALSE) In my table, how do I output only the words and not the words with the quotations? -- View this message in context: http://r.789695.n4.nabble.com/How-to-unquote-

Re: [R] Line numbers in Sweave

On 03/11/2010 11:52 AM, Yihui Xie wrote: Well, I know this is true for the default Sweave(), and my problem actually comes from the pgfSweave package: I used some tricks to "cheat" R to parse and deparse the code chunks so that the output can be automatically formatted (and preserve the comments,

Re: [R] Tukey's table

Hi: Try this: Trat <- c(2:30) # number of treatments gl <- c(2:30, 40, 60, 120) # Write a one-line 2D function to get the Tukey distribution quantile: f <- function(x,y) qtukey(0.95, x, y) outer(Trat, gl, f) It's slow (takes a few seconds) but it seems to work. HTH, Dennis

Re: [R] "non-numeric argument to binary operator" error while reading ncdf file

Charles Novaes de Santana wrote: > Thank you everybody for the help! The solution of my problem is here: > > http://climateaudit.org/2009/10/10/unthreaded-23/ > > "The mv variable is the designated NA for the variable and it appears that > somebody screwed that up in the file. This workaround worke

Re: [R] Line numbers in Sweave

Well, I know this is true for the default Sweave(), and my problem actually comes from the pgfSweave package: I used some tricks to "cheat" R to parse and deparse the code chunks so that the output can be automatically formatted (and preserve the comments, [1]). The price to pay for not being "hone

Re: [R] Lattice plots for images

On Nov 3, 2010, at 11:13 AM, Neba Funwi-Gabga wrote: Hello UseRs, I need help on how to plot several raster images (such as those obtained from a kernel-smoothed intensity function) in a layout such as that obtained from the lattice package. I would like to obtain something such as obtained

Re: [R] Lattice plots for images

Have you tried using the 'mai' argument to par()? Something like: par(mfrow=c(3,3), mai=c(0,0,0,0)) I've used this in conjunction with image() to plot raster data in a tight grid. -Matt On Wed, 2010-11-03 at 11:13 -0400, Neba Funwi-Gabga wrote: > Hello UseRs

Re: [R] dll problem with C++ function

Just a shot in the dark... Do you properly close the input/output files at the end of your function? If not and the file remains open, it may throw an error upon new attempt to read it. It is possible that dyn.unload, among other things, closes all open connections and hence upon re-load everything

[R] [klaR package] [NaiveBayes] warning message numerical 0 probability

Hi, I run R 2.10.1 under ubuntu 10.04 LTS (Lucid Lynx) and klaR version 0.6-4. I compute a model over a 2 classes dataset (composed of 700 examples). To that aim, I use the function NaiveBayes provided in the package klaR. When I then use the prediction function : predict(my_model, new_data). I g

Re: [R] longer object length is not a multiple of shorter object length

On Nov 3, 2010, at 11:17 AM, Stephen Liu wrote: - Original Message From: David Winsemius To: Stephen Liu Cc: r-help@r-project.org Sent: Wed, November 3, 2010 11:03:18 PM Subject: Re: [R] longer object length is not a multiple of shorter object length - snip - v <- 2*x + y + 1

Re: [R] longer object length is not a multiple of shorter object length

- Original Message From: David Winsemius To: Stephen Liu Cc: r-help@r-project.org Sent: Wed, November 3, 2010 11:03:18 PM Subject: Re: [R] longer object length is not a multiple of shorter object length - snip - >> v <- 2*x + y + 1 >> Warning message: >> In 2 * x + y : >> longer obje

[R] Lattice plots for images

Hello UseRs, I need help on how to plot several raster images (such as those obtained from a kernel-smoothed intensity function) in a layout such as that obtained from the lattice package. I would like to obtain something such as obtained from using the "levelplot" or "xyplot" in lattice. I current

Re: [R] longer object length is not a multiple of shorter object length

On Nov 3, 2010, at 11:00 AM, Stephen Liu wrote: Hi folks, I'm following An Introduction to R http://cran.r-project.org/doc/manuals/R-intro.html#R-and-statistics to learn R. Coming to; 2.2 Vector arithmetic v <- 2*x + y + 1 Warning message: In 2 * x + y : longer object length is not a mul

[R] longer object length is not a multiple of shorter object length

Hi folks, I'm following An Introduction to R http://cran.r-project.org/doc/manuals/R-intro.html#R-and-statistics to learn R. Coming to; 2.2 Vector arithmetic > v <- 2*x + y + 1 Warning message: In 2 * x + y : longer object length is not a multiple of shorter object length What does it mean?

[R] dll problem with C++ function

Dear fellow R-users, I have the problem of being unable to repeatedly use a C++-function within R unless by dyn.unloading/dyn.loading it after each .C call. The C++-code (too large to attach) compiles without problem using R CMD SHLIB. It loads (using dyn.load("myfun.so")) and executes (via .

[R] Granger causality with panel data (econometrics question)

Hi folks, I am trying to perform a Granger causality analysis with panel data. There are some packages around for panel data analysis and Granger causality. However, I have found neither a package for both panel data and Granger causality nor any R procedures (homogenous/heterogenous causality

Re: [R] optim works on command-line but not inside a function

As the error message says, the values of your function must be finite in order to run the algorithm. Some part of your loop is passing arguments (inits maybe... you only tried (0,0) in the CLI example) that cause IRT.llZetaLambdaCorrNan to be infinite. -- Jon

[R] optim works on command-line but not inside a function

Dear all, I am trying to optimize a logistic function using optim, inside the following functions: #Estimating a and b from thetas and outcomes by ML IRT.estimate.abFromThetaX <- function(t, X, inits, lw=c(-Inf,-Inf), up=rep(Inf,2)){ optRes <- optim(inits, method="L-BFGS-B", fn=IRT.llZetaLa

Re: [R] memory allocation problem

The optim function is very resource hungry. I have had similar problems in the past when dealing with extremely large datasets. What is perhaps happening is that each 'step' of the optimization algorithm stores some info so that it can compare to the next 'step', and while the original vector m

Re: [R] install vegan

On 03.11.2010 10:39, Carolin wrote: Dear all, I am trying to install Vegan, but I allways get the following error message: Warning in install.packages(choose.files("", filters = Filters[c("zip", : 'lib = "C:/Programme/R/R-2.12.0/library"' is not writable Error in install.packages(choose.f

Re: [R] Colour filling in panel.bwplot from lattice

Am 03.11.2010 12:52 (UTC+1) schrieb Deepayan Sarkar: On Wed, Nov 3, 2010 at 4:25 PM, Rainer Hurling wrote: Am 03.11.2010 10:23 (UTC+1) schrieb Deepayan Sarkar: On Wed, Nov 3, 2010 at 4:11 AM, Dennis Murphywrote: Hi: I don't know why, but it seems that in bwplot(voice.part ~ height, da

Re: [R] Drawing circles on a chart

Thanks Barry ... actually the intention was to have areas of the circle depicting the value (radius imputed) -Original Message- From: b.rowling...@googlemail.com [mailto:b.rowling...@googlemail.com] On Behalf Of Barry Rowlingson Sent: 03 November 2010 15:02 To: Santosh Srinivas Cc: r-help@

Re: [R] Colour filling in panel.bwplot from lattice

On Wed, Nov 3, 2010 at 4:25 PM, Rainer Hurling wrote: > Am 03.11.2010 10:23 (UTC+1) schrieb Deepayan Sarkar: >> >> On Wed, Nov 3, 2010 at 4:11 AM, Dennis Murphy  wrote: >>> >>> Hi: >>> >>> I don't know why, but it seems that in >>> >>> bwplot(voice.part ~ height, data = singer, >>> main = "NOT THE

Re: [R] Recoding -- test whether number begins with a certain number

On Wed, Nov 3, 2010 at 10:01 AM, Marcel Gerds wrote: >  Dear R community, > > I have a question concerning recoding of a variable. I have a data set > in which there is a variable devoted to the ISCO code describing the > occupation of this certain individual > (http://www.ilo.org/public/english/b

Re: [R] Multiple imputation for nominal data

The aregImpute function in the Hmisc package can do this through predictive mean matching and canonical variates (Fisher's optimum scoring algorithm). Frank - Frank Harrell Department of Biostatistics, Vanderbilt University -- View this message in context: http://r.789695.n4.nabble.com/Mu

[R] bad optimization with nnet?

Hy, I try to give an example of overfitting with multi-layer perceptron. I have done following small example : library(nnet) set.seed(1) x <- matrix(rnorm(20),10,2) z <- matrix(rnorm(10),10,1) rx <- max(x)-min(x) rz <- max(z)-min(z) x <- x/rx z <- z/rz erreur <- 10^9 for(i in 1:100){ temp.mo

Re: [R] density() function: differences with S-PLUS

Dear Joshua, first of all, thank you very much for reply. I hoped that someone who's familiar with both S+ and R can reply to me, because I spent some hours to looking for a solution. If someone else would try, this is the SPLUS code and output, while below there is the R code. I obtain the same

Re: [R] Colour filling in panel.bwplot from lattice

Am 03.11.2010 10:23 (UTC+1) schrieb Deepayan Sarkar: On Wed, Nov 3, 2010 at 4:11 AM, Dennis Murphy wrote: Hi: I don't know why, but it seems that in bwplot(voice.part ~ height, data = singer, main = "NOT THE RIGHT ORDER OF COLOURS\n'yellow' 'blue' 'green' 'red' 'pink' 'violet' 'brown' 'gold'"

[R] Tukey's table

Hi, I'm building Tukey's table using qtukey function. It happens that I can't get the values of Tukey's one degree of freedom and also wanted to eliminate the first column. The program is: Trat <- c(1:30) # number of treatments gl <- c(1:30, 40, 60, 120) # degree free

Re: [R] save() with 64 bit and 32 bit R

On Wed, 3 Nov 2010, Andrew Collier wrote: hi, i have been using a 64 bit desktop machine to process a whole lot of data which i have then subsequently used save() to store. i am now wanting to use this data on my laptop machine, which is a 32 bit install. i suppose that i should not be surprise

Re: [R] save() with 64 bit and 32 bit R

Andrew Collier wrote: hi, i have been using a 64 bit desktop machine to process a whole lot of data which i have then subsequently used save() to store. i am now wanting to use this data on my laptop machine, which is a 32 bit install. i suppose that i should not be surprised that the 64 bit dat

[R] Recoding -- test whether number begins with a certain number

Dear R community, I have a question concerning recoding of a variable. I have a data set in which there is a variable devoted to the ISCO code describing the occupation of this certain individual (http://www.ilo.org/public/english/bureau/stat/isco/isco88/major.htm). Every type of occupation begi

[R] install vegan

Dear all, I am trying to install Vegan, but I allways get the following error message: Warning in install.packages(choose.files("", filters = Filters[c("zip", : 'lib = "C:/Programme/R/R-2.12.0/library"' is not writable Error in install.packages(choose.files("", filters = Filters[c("zip", :

[R] boxplot of timeseries with different lengths

Hello List, I have a time serie of observations representing the activity of some users in different time periods, like: > table(obs1) user1 user2 user3 user31 user33 user4 user5 user6 user7 user8 user82 user83 user85 user9 1 1 3 1 1 1

Re: [R] Drawing circles on a chart

On Wed, Nov 3, 2010 at 2:07 AM, Santosh Srinivas wrote: > Dear Group, > Inside each cell there should be a circle (sphere preferable) with radius of > mod(data value). The color should be either red or green depending on -ve or > +ve and the intensity should be based on the value of the datapoint

  1   2   >