[R] snowfall

Hello, Just wondering why I am unable to run this in parallel. A dput of my dataset is attached at the end. Please use to create my data object. I want to run this function in parallel (not sure if this is an efficient implementation): #Function to calculate the time to maturity for the option r

[R] From polynomial to function

Dear all. I would like to use legendre polynomials which is something pretty easy in R. x<-legendre.polynomials(2)[[3]] > x -0.5 + 1.5*x^2 > str(x) Class 'polynomial' num [1:3] -0.5 0 1.5 As you can see from the code above str(x) returns that x is of class polynomial. I want to use that polyn

Re: [R] ThinkCell type waterfall charts in R?

On 01/12/2011 05:54 AM, ang wrote: Hi Jim, I looked through the plotrix documentation, and the waterfall plot comes from the stackpoly function right? I'm not sure if I can modify the stackpoly to create the plot I want, since stackpoly is a line plot and fills the area under with color. I ha

Re: [R] From polynomial to function

Alaios wrote: > x<-legendre.polynomials(2)[[3]] >> x > -0.5 + 1.5*x^2 > >> str(x) > Class 'polynomial' num [1:3] -0.5 0 1.5 > > As you can see from the code above str(x) returns that x is of class > polynomial. I want to use that polynomial as a function. The reason for > that is that I would b

Re: [R] plot: skip a range of axis

On 01/12/2011 03:46 PM, Yuan Jian wrote: Hi, I am using plot to show scatter points in 2_D. in my data, there is no data between -1 and +1 in x-axis. I want to skip this region, i.e. x axis becomes [-Inf:-1, 1:Inf]. can any one tell me how to do? Hi Yu, Try the gap.plot function in the plotrix

[R] vector or list of matrices corresponding to an observation

Dear all, I observe for each observation several joint distributions of two multinomial random variables (5x5 matrices). Right now, data are arranged so that I have 20 columns for each joint distribution by observation, which is not practical. I would like to work with matrices that would be inde

Re: [R] list concatenation

Bert Gunter writes: > Lists are (isomorphic to) trees with (possibly) labelled nodes. A > completely general solution in which two trees have possibly different > topologies and different labels would therefore involve identifying > the paths to leaves on each tree, e.g. via depth first search us

Re: [R] Problems creating a PNG file for a dendrogram: "Error in plot.window(...) : need finite 'xlim' values"

That was a simple solution. Turns out you were correct, "plot(p)" was the problem. Simply removing it and everything worked perfectly. Thanks Bill, Peter and David for your help. On Jan 11, 2011, at 8:29 PM, David Winsemius wrote: On Jan 11, 2011, at 9:27 PM, David Winsemius wrote: >

Re: [R] how to change strip text of effect plot

On Wed, Jan 12, 2011 at 10:36 AM, Joshua Wiley wrote: > Hmm, I felt like it was the clearest, most direct way.  Then again, > things of this nature (overide defaults using arguments rather than > changing what you feed the functions) do seem to be a common request > both for lattice and ggplot2.  

Re: [R] Degrees of freedom

Hi: Look at the links in the following blog entry: http://blog.lib.umn.edu/moor0554/canoemoore/2010/09/lmer_p-values_lrt.html and this discussion, found on the R wiki: http://rwiki.sciviews.org/doku.php?id=guides:lmer-tests Also see Ben Bolker's GLMM wiki page, which discusses many of the unreso

[R] R syntax for 95% prediction interval on a left-truncated normal variable

Hi all, I am searching for a R procedure to calculate a 95% prediction interval (i.e. the interval in which following observations of the variable will occur, not a confidence interval) on a variable that is the natural logarithm of a ratio that is always equal ore superior to 1, so the natural-lo

Re: [R] how to change strip text of effect plot

Dear Deepayan, Josh, and Ronggui, I've recently changed plot.eff() so that it returns an object, normally "printed" by print.plot.eff(). You can therefore manipulate the lattice object, as Deepayan suggests. This change is currently in the development version of the effects package on R-Forge,

[R] graphics: 3D regression plane

Hello Masters, wishing you all a great 2011 I was also going to ask if anyone knows a quick and efficient way to plot a regression plane (z~x*y). I have tried the regr2.plot{HH} function but it is only an educational tool and has poor graphical properties. I also tried to run the following script

[R] Issue loading and executing own function.R with JRI, any ideas?

Hey, how did you manage to load and call your own function with JRI? I managed to execute the build-in R functions with JRI, but when I call r.eval("load(path-to-file)") or r.eval("source(path-to-file)"), but my java program terminates :( Thanks for your help! -- View this message in context:

Re: [R] Multilevel pseudo maximum likelihood

Caterina, Did you get an answer to this question? I'm trying to do something similar. Jason -- View this message in context: http://r.789695.n4.nabble.com/Multilevel-pseudo-maximum-likelihood-tp878413p3213583.html Sent from the R help mailing list archive at Nabble.com. __

[R] Sum by column

Dear List, I have a question of convenience, I am looking to sum the values of one column based on another column - a example may help explain better! ED ECOCODE 21.809467 AA0101 36.229566 PA1201 51.861284 PA1201 11.36232PA1201 27.264634 PA12

Re: [R] plot: skip a range of axis

thanks Jim, I found gap.plot seperates x axis or y axis into two boxes. do you know any plot tool that can skip a range in x-axis or y-axis without lines?   regards YU --- On Wed, 12/1/11, Jim Lemon wrote: From: Jim Lemon Subject: Re: [R] plot: skip a range of axis To: "Yuan Jian" Cc: r-help

[R] Weighted Likelihood Estimation of NIG Dist.

Hi All, I've put together some script which gives me the parameters of a modified NIG distribution. Now I'd like to include a weighted vector within the maximization function. The code below gives me parameter estimates which are implicitly equally weighted. nig.par.fit <- try(optim(vega,

[R] non-parametric discriminant analysis

Hi, I used linear discriminant analysis (lda) to classify and to cross-validate samples of two plant species based on morphometric data and to identify the variables that best discriminate between the two species. Because some of the variables are not, and can not be transformed to be, normal

Re: [R] Sum by column

There are two functions you need to become familiar with: ?tapply ?ave If you wanted these summed values to be placed in another column of the same dataframe, you would use ave. If you wanted a new structure (somewhat shorter) you would use tapply with sum as the function. E. g: tapply(eco

Re: [R] Sum by column

Hi Peter, R has some fairly flexible ways of passing values of some variable (X) by another (the INDEX) to different FUNctions. Here is an example using your data: ## your email data, in convenient form dat <- structure(list(ED = c(21.809467, 36.229566, 51.861284, 11.36232, 27.264634, 12.261986,

[R] flexmix: predictions on new data from flexmix object

Dear R Users, R Core Team, I currently wonder how to predict the probability of an event with new data resulting from a finite mixture. I read the documentation of the flexmix package and the examples of applications provided on CRAN but I could not find how to predict (except "manually" but

[R] Integrate and subdivisions limit

Dear all, I have some issues with integrate in R thus I would like to request your help. I am trying to calculate the integral of f(x)*g(x). The f(x) is a step function while g(x) is a polynomial. If f(x) (step function) changes its value only few times (5 or 6 'steps') everything is calulated ok

[R] Fwd: vector or list of matrices corresponding to an observation

You're right Jim. I have now started to work with list of matrices. I'm not sure that it is a clean and nice code but it works. I was wondering if there was a synthetic but advanced tutorial on list() and associated functions? Stephane -- Forwarded message -- From: jim holtman

[R] Outputting csv file from dataframe with columns in a particular order

I have a dataframe with columns "ID",'date","estimate","actual" (but not necessarily in that order - I do a merge somewhere and that somehow messes up the order of the columns). How can I output it to a csv file with the columns in the order that I want? Thanks. _

Re: [R] debug biglm response error on bigglm model

Thank you, Greg. The issue was in the simulation logic, where one of the values was not changing correctly for some iterations... On Jan 10, 3:20 pm, Greg Snow wrote: > Not sure, but one possible candidate problem is that in your simulations one > iteration ended up with fewer levels of a facto

Re: [R] extracting more information from optim in R?

I have no experience with writing C code, but if I have such problems in R code, I add a line to my function which prints the values to the console: eg: fr <- function(x) { ## Rosenbrock Banana function x1 <- x[1] x2 <- x[2] cat (paste(x1, x2, "\n")) 100 * (x2 - x1 * x1)^2 +

Re: [R] Sum by column

David and Josh, Thanks very much for your help, it is much appreciated. Peter On 12 Jan 2011, at 14:28, David Winsemius wrote: There are two functions you need to become familiar with: ?tapply ?ave If you wanted these summed values to be placed in another column of the same dataframe, you w

[R] Bootstrapping to Correct Standard Errors in Two-Stage Least Square Estimation

Dear friends I want to estimate an equation using two-stage least square but suspect that the model suffers from autocorrelation. Can someone please advise how to implement bootstrapping method in order to calculate the correct standard errors in R? Thank you. Kind regards Thanaset -- View t

[R] Basic Stars Plot - help ..

Hi there Rers I am trying a very basic stars plot: x<-matrix(c(1,4,3,1.1,2,3,4,3,1,1,5,2), ncol = 3, byrow = TRUE, >> dimnames=list(c("a","b","c","d"),c("x","y","z"))) > > stars(x, draw.segments = TRUE, radius=TRUE) > > Can anyone explain what I am seeing there - EACH of my plots should have 3 co

[R] How to define values for color distribution in the package gplots and function heatmap.2

Hi, This question is about the package gplots and the function heatmap.2. I'm not a programmer and I did not understand the answers I found when I googled it. Thank you in advance. A similar question was asked with title "heatmap color distribution" in 2005 and got the answer to use breaks. When

Re: [R] Outputting csv file from dataframe with columns in a particular order

Hi! Let's say your data.frame is called df and that you want column 1, then column 4, then 3 and then 2: df <- data.frame(ID=LETTERS[1:5], date=rnorm(5), estimate=rnorm(5), actual=rnorm(5)) write.csv(df[c(1,4,3,2)], file="df.csv") HTH, Ivan Le 1/12/2011 16:16, analys...@hotmail.com a écrit :

Re: [R] snowfall

You forgot to load the required packages on the client nodes by sfLibrary(fCalendar) sfLibrary(fractalrock) and you really should not tryCatch without evaluating the errors for yourself. Best wishes, Uwe Ligges On 12.01.2011 09:47, Santosh Srinivas wrote: Hello, Just wondering why I am u

Re: [R] Outputting csv file from dataframe with columns in a particular order

On 2011-01-12 07:16, analys...@hotmail.com wrote: I have a dataframe with columns "ID",'date","estimate","actual" (but not necessarily in that order - I do a merge somewhere and that somehow messes up the order of the columns). How can I output it to a csv file with the columns in the order that

[R] Odp: Outputting csv file from dataframe with columns in a particular order

Hi r-help-boun...@r-project.org napsal dne 12.01.2011 16:16:16: > I have a dataframe with columns "ID",'date","estimate","actual" (but > not necessarily in that order - I do a merge somewhere and that > somehow messes up the order of the columns). If you have datafreme with column order a, b, c,

Re: [R] Sum by column

Or with ddply : library(plyr) dat <- structure(list(ED = c(21.809467, 36.229566, 51.861284, 11.36232, 27.264634, 12.261986, 46.519313, 7.815376, 2.810428, 13.478372, 35.670182, 27.128715, 19.010294, 15.475368, 18.597983, 29.292615, 6.749846, 14.981488, 14.93511, 14.93511, 21.040785, 8.271615, 12.94

Re: [R] graphics: 3D regression plane

On Jan 12, 2011, at 6:10 AM, Federico Bonofiglio wrote: Hello Masters, wishing you all a great 2011 I was also going to ask if anyone knows a quick and efficient way to plot a regression plane (z~x*y). There are many. There are limitations to using the ?? operator in that it only brings

Re: [R] Basic Stars Plot - help ..

On 12.01.2011 15:53, JP wrote: Hi there Rers I am trying a very basic stars plot: x<-matrix(c(1,4,3,1.1,2,3,4,3,1,1,5,2), ncol = 3, byrow = TRUE, dimnames=list(c("a","b","c","d"),c("x","y","z"))) stars(x, draw.segments = TRUE, radius=TRUE) Can anyone explain what I am seeing there - EAC

[R] adonis, amova and haplotype frequency

Dear All, I'd like to perform adonis (from the vegan package) rather than amova (in ade4) on some haplotype data, as I have crossed factors. Is there a simple way to tweak the source to allow weights (haplotype frequencies) in a similar way to amova? Best Simon __

Re: [R] R not recognized in command line

Im sorry for the late reply. The output to echo %PATH% : C:\GTK\bin; C:\Program Files\MiKTeX 2.8\miktex\bin ;C:\Windows\system32 ;C:\Windows ;C:\Windows\System32\Wbem; C:\Windows\System32\WindowsPowerShell\v1.0\; C:\Program Files\MATLAB\R2008a\bin; C:\Program Files\MATLAB\R2008a\bin\win32; C:\Pr

Re: [R] Formatted output with alternating format at different rows

thx a lot Jim. sprintf solved my problem. Ray On Mon, Jan 3, 2011 at 3:40 PM, jim holtman wrote: > 'sprintf' if your friend: > >> dummy3 = c(1.1, 2.2, 3.3) >> dummy4 = c(4.4, 5.5, 6.6, 7.7) >> dummy2 = c(8.8, 9.9) >> >> cat(sprintf("%5.1f%6.2f%7.3f\n", dummy3[1], dummy3[2], dummy3[3])) >  1.1  

Re: [R] R not recognized in command line

> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf Of Aaditya Nanduri > Sent: Wednesday, January 12, 2011 8:44 AM > To: Uwe Ligges > Cc: r-help@r-project.org; spec...@stat.berkeley.edu > Subject: Re: [R] R not recognized in command li

Re: [R] R not recognized in command line

On 12.01.2011 18:13, Daniel Nordlund wrote: -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Aaditya Nanduri Sent: Wednesday, January 12, 2011 8:44 AM To: Uwe Ligges Cc: r-help@r-project.org; spec...@stat.berkeley.edu Subject: Re:

[R] Require

I think that the "quietly" argument in "require" isn't working > require('JumboShrimp', quietly=TRUE) Warning in library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE, : there is no package called 'JumboShrimp' > By the way, the behavior is the same with options(warn

Re: [R] Require

On 12.01.2011 18:53, Gene Leynes wrote: I think that the "quietly" argument in "require" isn't working require('JumboShrimp', quietly=TRUE) Warning in library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE, : there is no package called 'JumboShrimp' ?require

Re: [R] Require

Gene Leynes writes: > I think that the "quietly" argument in "require" isn't working > > > require('JumboShrimp', quietly=TRUE) > Warning in library(package, lib.loc = lib.loc, character.only = TRUE, > logical.return = TRUE, : > there is no package called 'JumboShrimp' > > Isn't quietly mean

[R] Don´t know what test i have to use

Hello, I´m starting with my PhD and I have to stop because i got a little knowledge in R and statistics. I´ve got a model of this kind: binary response variable: prevalence of infection (0/1) 3 categorical independent variables: sex, month and name of the area I was trying with a full model li

[R] Metafor vs Meta vs Spreadsheet: wrong numbers

Hello, I experimented the Metafor and Meta packages in the scope of replacing Excel for meta-analysis. I performed the first working example provided in Michael Borenstein's book "Introduction to Meta-Analysis" with Excel, Metafor and Meta. The numbers given by my spreadsheet, which I vali

Re: [R] R not recognized in command line

Although it could easily be "user error", I never got Rpy or Rpy2 working an any sort of reliable way. However I did learn a couple of things about the Windows PATH First, (as others have mentioned) it's easiest to modify the PATH through the Windows GUI that comes up when you right click "My Com

Re: [R] A question on dummy variable

Thanks Gabor and other for their input. I admit that I must have placed some reproducible codes on what I wanted. However it was actually in my mind however I restrained because it was not any R related query rather a general Statistics related. Here I am using dummy variables in ***Time series co

Re: [R] Require

I read the help first, and read it again now, but I still don't see why the warning was generated. The help seems to state clearly that the user can suppress warnings "most often, no errors/warnings are printed if package loading fails." If the package doesn't exist, then it would fail to loa

Re: [R] Don´t know what test i have to use

On Jan 12, 2011, at 12:51 PM, gaiarrido wrote: Hello, I´m starting with my PhD and I have to stop because i got a little knowledge in R and statistics. I´ve got a model of this kind: binary response variable: prevalence of infection (0/1) 3 categorical independent variables: sex, month and

Re: [R] Don´t know what test i have to use

Hi, That is basically correct. You can specify the link as logit (see my example), but that is the default so you do not strictly need to in this case. II would encourage you to keep your variables (prevalencia, edad, sexo, mes) stored in a data frame, in which case you would add the data = argu

Re: [R] Integrate and subdivisions limit

> Dear all, > I have some issues with integrate in R thus I would like to request > your help. I am trying to calculate the integral of f(x)*g(x). > The f(x) is a step function while g(x) is a polynomial. > If f(x) (step function) changes its value only few times (5 or 6 'steps') > everything is ca

[R] 2d plot with modification of plotting symbol to indicate third dimension.

I would like to plot 3-dimensional data on a two-dimensional scatter-plot. Is there a way I can automatically modify the plot symbol (e.g. changing size or color) to indicate the value of a third variable? E.g. How can I plot weight vs. age and indicate the value of muscle mass for each value wei

Re: [R] 2d plot with modification of plotting symbol to indicate third dimension.

You don't give an example, but in general you can use a vector for cex with the values proportional to the third variable. Same goes for color: col can be a vector, not just a single value. This has been discussed before on-list, and fairly recently. Sarah On Wed, Jan 12, 2011 at 2:19 PM, John

[R] metafor/ meta-regression

Hi I have tryed to do the meta-regression in metafor package, but I would like to get the standardized coefficients for each variable, however in command:   Ø  res<-rma.uni (yi, vi, method="REML", mods=~cota+DL+uso+gadiente+idade, data= turbidez)   I just have the coefficients no standardized (

Re: [R] Don´t know what test i have to use

Thanks very much both. I´m starting playing with it, i was a little afaid because it was part of my job, but now i've found it very funny. Josh, I've got just data for 3 representatives months, and it's not a priori rejectable that could be differences in the ratio of changes along the months bet

Re: [R] data frame subset too slow

Sorry for the late response. I was away for vacation and was unable to keep on working on the codes. Anyway, I was unable to provide *str* of that specific data since they are all in a big package with lots of inputs/outputs. Quickly gazing through the code, I narrowed them down (and made a ba

Re: [R] 2d plot with modification of plotting symbol to indicate third dimension.

Look at the symbols function for some options of doing what you suggest (you can also do a search for "bubble plot" for a couple of other implementations). If you want to go a bit further than what symbols does for you then look at the my.symbols function in the TeachingDemos package. -- Greg

[R] aggredating date data

I tried a date by date forecast of a time series and it seems to be too wild. How can I aggregate the date into weeks or months as required? Thanks. The input looks like ID datadate("-MM-DD") value_for_day -- ---- -- --

Re: [R] Outputting csv file from dataframe with columns in a particular order

Thanks to all who responded. On Jan 12, 10:34 am, Peter Ehlers wrote: > On 2011-01-12 07:16, analys...@hotmail.com wrote: > > > I have a dataframe with columns "ID",'date","estimate","actual" (but > > not necessarily in that order - I do a merge somewhere and that > > somehow messes up the order

[R] help in calculating ar on ranked vector

I was using ar(stats) to calculate autoregressive coefficient. It works on vector z, but it will not work on vector rz <-rank (z, ties.method="average").  What did I miss? Any info will be greatly appreciated.  TIA [[alternative HTML version deleted]] __

Re: [R] Help with Data Transformation

Hi John, Thank you for your patience. I was away for a State certification exam yesterday, so am just getting back to this. Reading through you response I believe I wasn't clear enough about what I'm trying to do. Your description seems to rearrange the matrix without grouping the analytica

[R] navigating in lists

Dear list members, I am stuck with navigating in a rather complicated list object. In general I would need a solution to access all first (or other) elements of the different sublists in one list: test=list(a=list(1,2),b=list(3,4),c=list(5,6)) like: test[[1:3]][[1]] which should result in

Re: [R] plot: skip a range of axis

You could always create a new vector, something like Xprime<- if x<-0 x else x-2 #not valid R code Thus mapping +1 to -1, and shifting everything else down. Fixing the x-tick labels is left as a homework problem :-) I'm assuming from your description that there are no y-values correspondi

Re: [R] 2d plot with modification of plotting symbol to indicate third dimension.

Google on "R Graph Gallery" to find examples with code. You can also almost certainly RSiteSearch() on appropriate keys to find a pre-existing function (which someone may provideyou on the list). Also: ?symbols -- Bert On Wed, Jan 12, 2011 at 11:19 AM, John Sorkin wrote: > I would like to plot

Re: [R] navigating in lists

> sapply(test, '[[', 1) a b c 1 3 5 -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of Jannis > Sent: Wednesday

Re: [R] Help with Data Transformation

Hi: This seems like a problem that is well suited for the cast() function in package reshape2. Here's a toy example: library(reshape2) df <- data.frame(idnum = rep(101:104, c(3, 2, 4, 3)), lab = unlist(sapply(c(3, 2, 4, 3), function(x) sample(LETTERS[1:6], x))),

Re: [R] Don´t know what test i have to use

... But I would think that month should be treated as a cyclical quantity, not as a factor with 12 independent levels, e.g. by transforming month to sin( 2*pi*monthNumber/12) . This assumes 1 year periodicity, which might not be right, of course. Time series methods could obviously be relevant he

[R] How to disable using "enter" key to exit the browser in debugging mode

Dear R, How can I disable using "enter" key to exit the browser() in debug mode? I would love to have this option because it is so annoying to jump out of the debugging mode unexpectedly when I don't want to. I guess some of us have encouraged at least one of these situations, 1, Accidentally pre

[R] RNetCDF: retrieving variable names and units

Dear List, does anybody has experience with the RNetCDF package? I manage to open a connection and copy data from a ncdf file but would need a way to automatically retrieve variable names (ideally all of them from one file) and units from the file. Any ideas? Jannis __

Re: [R] RNetCDF: retrieving variable names and units

Hi Jannis, although I don't know how you'd do that with RNetCDF, with the ncdf package it's pretty easy: ncid = open.ncdf( 'file.nc' ) nvars = ncid$nvars for( ivar in 1:nvars ) print(paste("var number",ivar,"is named", ncid$var[[ivar]]$name, "and has units", ncid$var[[ivar]]$units )) Regards

Re: [R] RNetCDF: retrieving variable names and units

There are a number of functions in the package to inquire about the file contents. See library(help = RNetCDF). For example: library(RNetCDF) nc <- open.nc("file.nc") var.inq.nc(nc, 0) $id [1] 0 $name [1] "longitude_U" $type [1] "NC_DOUBLE" $ndims [1] 1 $dimids [1] 1 $natts [1] 0 You can

[R] syntax for extending a line in a script??

Hello, A hopefully simple question. I use 'R' through emacs, but I suspect the following would occur with any manner of text editor: - my editor has a normally quite handy feature where it will automatically indent to the appropriate level when I start a new line. However, this occa

[R] speed up subsetting with certain conditions

Hi folks, I am working on a project that requires subsetting of a found file based on some known file. The known file contains several lines like below: chr132375463237547rs523104280+ chr132375493237550rs520975820+ chr245133264513327rs2976928

Re: [R] syntax for extending a line in a script??

On 1/12/2011 2:46 PM, Mike Williamson wrote: Hello, A hopefully simple question. I use 'R' through emacs, but I suspect the following would occur with any manner of text editor: - my editor has a normally quite handy feature where it will automatically indent to the appropriate le

Re: [R] speed up subsetting with certain conditions

On 1/12/2011 2:52 PM, Duke wrote: Hi folks, I am working on a project that requires subsetting of a found file based on some known file. The known file contains several lines like below: chr132375463237547rs523104280+ chr132375493237550rs520975820+ chr

Re: [R] syntax for extending a line in a script??

On Jan 12, 2011, at 5:46 PM, Mike Williamson wrote: Hello, A hopefully simple question. I use 'R' through emacs, but I suspect the following would occur with any manner of text editor: - my editor has a normally quite handy feature where it will automatically indent to the appropri

Re: [R] aggredating date data

?aggregate I would transfer your date character sting into a date object (as.POSIXct) and then extract month or week numbers (?format) from this vector and use them as indices for the aggregate function. There may be more elegant ways though HTH Jannis --- analys...@hotmail.com schrieb a

[R] Grouped bars in barplot

Dear all, I am trying to make a barplot with clustered pairs of bars, using class=numeric data and the following command: barplot(c(bline_precip[10,9], bline_runoff[10,9], cccma_precip[10,9], cccma_runoff[10,9], csiro_precip[10,9], csiro_runoff[10,9], ipsl_precip[10,9], ipsl_runoff[10,9], m

Re: [R] navigating in lists

Or, if for some reason the lists differ in length... test=list(a=list(1,2),b=list(3,4),c=list(5,6,7)) picker <- function(x, i) { if(length(x)>=i) x[[i]] else NA } pick <- function(list,i) { sapply(list, function(x) picker(x, i)) } > pick(test, 1) a b c 1 3 5 > pick(test, 2) a b c 2 4 6

[R] Multivariate autoregressive models with lasso penalization

I wish to estimate sparse causal networks from simulated time series data. Although there's some discussion about this problem in the literature (at least a few authors have used lasso and l(1,2) regularization to enforce sparsity in multivariate autoregressive models, e.g., http://user.cs.tu-berli

Re: [R] speed up subsetting with certain conditions

On 1/12/11 6:12 PM, Martin Morgan wrote: The Bioconductor project has many tools for dealing with sequence-related data. With the data k <- read.table(textConnection( "chr132375463237547rs523104280+ chr132375493237550rs520975820+ chr24513326451332

Re: [R] syntax for extending a line in a script??

Thanks to Brian & David! For some reason, I'd thought to use "{", but not "(". I guess I can chalk that up to a slow brain. Interestingly (and I didn't bother to figure out why), the "{" didn't work for me. But I tried a few cases, and "(" always seems to work, so far.

Re: [R] How to disable using "enter" key to exit the browser in debugging mode

That also drives me crazy! I don't have that problem when I use the StatEt plug-in for Eclipse. Of course, using a new IDE is a big undertaking, but I can assure you: it's worth it! This is just one small benefit. On Wed, Jan 12, 2011 at 4:00 PM, Feng Li wrote: > Dear R, > > How can I disable

[R] 2d plot with modification of plotting symbol to indicate third dimension.

I would like to plot 3-dimensional data on a two-dimensional scatter-plot. Is there a way I can automatically modify the plot symbol (e.g. changing size or color) to indicate the value of a third variable? E.g. How can I plot weight vs. age and indicate the value of muscle mass for each value wei

[R] Openbugs and rbugs on mac with wine

Hello list, I’ve been trying to get OpenBUGS running on my mac using the wine emulator. I can run Openbugs just fine by doing: wine ~/OpenBUGS312/OpenBUGS.exe In the terminal, so OpenBUGS works. When I try to run the schools example using rbugs(), the OpenBUGS process starts in wine, but i

[R] Openbugs and rbugs on mac with wine

Hello list, I’ve been trying to get OpenBUGS running on my mac using the wine emulator. I can run Openbugs just fine by doing: wine ~/OpenBUGS312/OpenBUGS.exe In the terminal, so OpenBUGS works. When I try to run the schools example using rbugs(), the OpenBUGS process starts in wine, but it ju

Re: [R] Grouped bars in barplot

Tena koe Steve Convert your data into a matrix: dataMat <- matrix(c(bline_precip[10,9], bline_runoff[10,9], cccma_precip[10,9], cccma_runoff[10,9], csiro_precip[10,9], csiro_runoff[10,9], ipsl_precip[10,9], ipsl_runoff[10,9],

[R] easy loop question

Hi everyone, I am new in R and programming. I have tried to remove the values out of range in some variables using a loop: 1) var <- names(est8vo[, 77:83]) # I got the variable names > var [1] "p16.1" "p16.2" "p16.3" "p16.4" "p16.5" "p16.6" "p16.7" for (i in 1:7) { var.i <- var[i] est8

[R] Rotated, Right-Justified Labels for Shortened Tick Marks

Hello R-help, I'm trying to make a fairly simple plot axis that goes something like this: plot(-10:10,-10:10, yaxt='n') axis(side=2, las=1, hadj=1, tck=-.01, cex.axis=.6) ...but as you can see, the labels are not close enough to the y-axis (where I want them... to save space for publication).

Re: [R] aggredating date data

I like the zoo package, and there are several helpful examples. library(zoo) You can easily convert your data into a zoo object using I was actually just doing this using this function: LoadReturnData=function(x){ ret = read.csv(x) ret = zoo(ret[ , -1], as.Date(ret[ , 1])) colnames(ret

[R] What does the shell() command do?

Dear R community, I am trying to understand what the shell() function does. An example is: xfile <- shell(paste("dir/b ", paste(directory.folder,"file.name",sep="")),intern=T) I'm afraid I wasn't able to completely understand the explanation under the Help files. Thanks for your help! Leanne

[R] standard errors in johansen test

Dear all, I have a question. How to get the standard errors of alpha and beta when using "ca.jo" to test cointergration? In the paper by Bernhard Pfaff and Kronberg im Taunus “VAR, SVAR and SVEC Models: Implementation Within R Package” pp.24-25. The standard errors are listed on the table

[R] unicode&pdf font problem

Dear List, I would like to print a plot into pdf. The problem is that the character \U0171 is replaced by a simple 'u' (i.e. without accents) in the pdf file. Example: # this works fine plot(1,type="n") text(1,1,"print \U0171") # this fails pdf("trial.pdf") plot(1,type="n") text(1,1,"print \U01

[R] question about svm(e1071)

Dear all, I executed svm calculation using e1071 library with a microarray data (http://www.iu.a.u-tokyo.ac.jp/~kadota/R/data_Singh_RMA_3274.txt). Then, I shuffled the data samples and executed svm calculation again. The results of 2 calculation were different (in SV, coefs and weights). I attac

[R] Repeating value occurence

How can achieve this in R using seq, or rep function c(-1,0,1,0,-1,0,1,0,-1,0) The range value is between-1 and 1, and I want it such that there could be n number of points between -1 and 1 Anyone? Please help Thanks Rusty [[alternative HTML version deleted]] _

Re: [R] Help with Data Transformation - RESOLVED

Hi Dennis, SOLVED!!! My thanks to both you, John, and others who chimed in. Took a little more digging before finally working, but it is working! Here's a little more of what I did and the ultimate resolutions: I installed the reshape2 library from CRAN. Executed > (.packages()) [1] "reshape2"

Re: [R] easy loop question

On Jan 12, 2011, at 10:54 PM, Sebastián Daza wrote: Hi everyone, I am new in R and programming. I have tried to remove the values out of range in some variables using a loop: 1) var <- names(est8vo[, 77:83]) # I got the variable names > var [1] "p16.1" "p16.2" "p16.3" "p16.4" "p16.5"

  1   2   >