[R] Partial loop?

My dear R buddies, I'm writing a loop program like this: for(i in 1:n){ for(j in 1:i){ ... } } I wonder if there is any simple apply()-like function to make the loop a little bit easier and faster. Thanks a lot. Best wishes, -- 彭河森 Hesen Peng http://hesen.peng.googlepages.com/

[R] compressing the plot's white space

I need to fit a graph into a column of a 2-column paper. I found that just specifying width and height parameters (3.2in x 3.5in) to plot doesn't decrease the fonts of the main title, axis titles, and labeling numbers, and tick sizes. So I have to add cex to all labels and titles and manage

Re: [R] 50993 point distance matrix, too big to as.matrix, looking for another way to calculate point-level summary

Hi, If you are only interested in row means, you can work the distance matrix at the c level. You might like to adapt this post: http://tolstoy.newcastle.edu.au/R/e6/devel/09/04/1378.html Romain On 06/26/2009 09:40 PM, leif olson wrote: Hello, Im working on a 50933 point count bird

Re: [R] Partial loop?

Finally I ended up doing this: temp - expand.grid(1:n,1:n) temp-temp[temp[,1]temp[,2]] apply(temp,1, ... ) and it seems much faster :) On Sat, Jun 27, 2009 at 3:26 PM, Patrick Burnspbu...@pburns.seanet.com wrote: See 'The R Inferno'. Patrick Burns patr...@burns-stat.com +44 (0)20 8525

[R] Compute correlation matrix for panel data with specific ordering

Hello All, I have a panel date - here a small-scale example: df - data.frame(cbind(rep(c(AUT,BEL,DEN,GER),4),cbind(rep(c(1999,2000,2001,2002),4)),sample(10,16,replace=T))) names(df) - c(country,year,x) SORT - c(GER,BEL,DEN,AUT) I need to compute the correlation between countries in the

Re: [R] questions about meta-analysis

Le samedi 27 juin 2009 à 13:02 +0800, sdzhangping a écrit : Dear R users: In the example of meta-analysis (cochrane, package rmeta), I can not found the p-value of Test for overall effect, and some other indices (Z, I, weight and et al). How can I get the these indices listed?

Re: [R] Compute correlation matrix for panel data with specific ordering

Serguei Kaniovski Serguei.Kaniovski at wifo.ac.at writes: df - data.frame(cbind(rep(c(AUT,BEL,DEN,GER),4), cbind(rep(c(1999,2000,2001,2002),4)),sample(10,16,replace=T))) names(df) - c(country,year,x) SORT - c(GER,BEL,DEN,AUT) I need to compute the correlation between countries in the

Re: [R] How to automatically placing a legend in empty space

Jason Rupert jasonkrupert at yahoo.com writes: At one point I believe I heard of an R package that would automatically find the most empty space in a plot, and then that answer could then be used to intelligently place a legend. Have a look a some of the plotting function in package Hmisc.

[R] Creating Excel-Charts via RDCOMClient

Hi R-users! I'm trying to create an easy Excel chart using the package RDCOMClient. The following example is working fine: #--- library(RDCOMClient) xlLocationAsObject -2 xlXYScatterSmoothNoMarkers - 73 ex -

Re: [R] questions about meta-analysis

The object a has log.estimates and selog.estimates: str(a) List of 10 $ logOR : num [1:7] -0.548 -1.804 -1.404 -0.357 -1.055 ... $ selogOR : num [1:7] 0.22 1.11 0.611 0.372 0.715 ... $ logMH : num -0.632 $ selogMH : num 0.16 snipped remainder of output. You can also look at the

[R] Regression; how to get t-values for all parameters estimates

Dear all, Even after a couple of hours looking at old messages I still haven't found a solution for my problem. I'm trying to fit an additive linear regression model with 2 effects, both fixed, to some dataset. The function contrasts(effectA) - contr.sum can gaurantee that the coefficients per

Re: [R] How do I define the method for gcheckboxgroup in gWidgets?

On Thu, Jun 25, 2009 at 8:29 AM, Bryan Hanson han...@depauw.edu wrote: Hi All... I¹m trying to build a small demo using gWidgets which permits interactive scaling and selection among different things to plot. I can get the widgets for scaling to work just fine. I am using gcheckboxgroup

Re: [R] parallel R?

For cross-validation, the caret package was designed to easily go between sequential and parallel processing (using nws, mpi or anything else). See the last examples in ?train. Max On Jun 26, 2009, at 11:28 AM, Michael comtech@gmail.com wrote: I guess when we move to Amazon AWS,

Re: [R] How do I define the method for gcheckboxgroup in gWidgets?

Thanks Michael... I was working by analogy to the gbuttons, so I was trying to ³add² the gcheckboxgroup, which is apparently not necessary (due to, I guess, the intrinsic differences between the widgets). The index thing I was just screwed up on! I have it working now. Nice package. Bryan

Re: [R] compressing the plot's white space

On 27/06/2009 3:06 AM, Alexy Khrabrov wrote: I need to fit a graph into a column of a 2-column paper. I found that just specifying width and height parameters (3.2in x 3.5in) to plot doesn't decrease the fonts of the main title, axis titles, and labeling numbers, and tick sizes. I normally

Re: [R] Creating Excel-Charts via RDCOMClient

The statconn tools allow you configure your application with the Excel client on one machine and the R server on another machine. In (D)COM, the D stands for distributed COM. Dr. Michael Wolf wrote: Hi R-users! I'm trying to create an easy Excel chart using the package RDCOMClient. The

Re: [R] questions about meta-analysis

At 06:02 27/06/2009, sdzhangping wrote: Dear R users: In the example of meta-analysis (cochrane, package rmeta), I can not found the p-value of Test for overall effect, and some other indices (Z, I, weight and et al). How can I get the these indices listed? library(rmeta) data(cochrane)

Re: [R] changing default arguments of a function and return the modified function as a result

see also `%but%.character` in the operators package. rnorm %but% list( mean = 3 ) function (n, mean = 3, sd = 1) .Internal(rnorm(n, mean, sd)) environment: namespace:stats baptiste [[alternative HTML version deleted]] __

Re: [R] How to read a specific dataset, not the entire data, from HDF5?

Thanks, Felix. You are watching new questions here. I will give a try to rgdal and let you know the result. BTW, it is surprising to me that there seems no R package handling HDF5 like PyTables does for Python community, considering HDF5 is as popular as NetCDF as a file format for scientific

Re: [R] panel.text and saving to pdf

If it is not the device then it's probably just the usage of an unexpected unit somewhere - as reported in another reply by another helper to your original question on R-help. Best, Uwe Ligges willem vervoort wrote: Hi Uwe, I should send these messages from my member e-mail number Thanks

Re: [R] Re gression; how to get t-values for all parameters estimates

gzf200 wrote: Even after a couple of hours looking at old messages I still haven't found a solution for my problem. I'm trying to fit an additive linear regression model with 2 effects, both fixed, to some dataset. The function contrasts(effectA) - contr.sum can gaurantee that the

[R] Looking for Rdonlp2

Dear All, I tried to download the package Rdonlp2 from the address given at CRAN Task View: Optimization and Mathematical Programming: http://arumat.net/Rdonlp2/ However, this link seems to be dead. Any ideas? Thanks in advance, Paul __

Re: [R] a plot of stacked boxes

check out 'map.market' in the 'portfolio' package. On Fri, Jun 26, 2009 at 11:27 PM, Osman Al-Radi osman.al.r...@gmail.comwrote: Dear Richard and David, Thanks for this reference. I looked into vcd and mosaic plot, it is a nice plot for investigating associations between two or more

[R] gradient fill of a grid.polygon

Following up on my previous post. I've managed to have the function return a gList rather than plot everything directly, but I get a rather obscure error message when I try to wrap the grobs in a gTree with a rotated viewport, Error in x$children[[i]] : attempt to select less than one element

[R] Bayesian Model Averaging

Two questions pertaining to the iBMA.glm procedure: 1) Does prior.param=1 force variables into the model as it does in S-Plus? If not, how can I force the first p variables into every model? 2) Is there some way in which to suppress the intercept as in the case of origin models?

[R] data formatting

All, I have three columns of data: id, date, hab. I am trying to set up a matrix that has the id as the rows, date as columns, and the hab value as the data values. Each id/date combination can only have one hab value. I would like for it to look something like this date 1, date 2,

Re: [R] Looking for Rdonlp2

Not sure what word in my prior message triggered the list-filter but this is what I got after a bit of cutting and pasting on the link offered: This directory does not contain anymore the source of DONLP2. The code is available for some users from its author, Peter Spellucci,

Re: [R] data formatting

Is this what you are after (uses the reshape package) # test data x - expand.grid(id=1:5, date=seq(as.Date('2009-01-01'), by='1 day', length=5)) x$hab - seq(nrow(x)) # add hab head(x) id date hab 1 1 2009-01-01 1 2 2 2009-01-01 2 3 3 2009-01-01 3 4 4 2009-01-01 4 5 5

Re: [R] how to interpret coefficients for a natural spline smooth function in a GLM

I have not seen a reply to this question, so I will offer a comment; someone who knows more than I may correct or add to my comments. There are many different kinds of splines. Perhaps the most common are B-splines, which sum to 1 inside their range of definition and are 0

[R] Defining class as a member of another class

Hi, When I define a new class (through setClass), members defined in representation argument doesn't seem to like a class. For example, if I do the following: setClass(NotWork,representation=(x=zoo)) It seems to me that representation members will take in only primitive type to R. Is there any

Re: [R] Defining class as a member of another class

Try this: setClass(zoo) [1] zoo setClass(Work,representation=(x=zoo)) [1] Work On Sat, Jun 27, 2009 at 10:01 PM, R_help Helprhelp...@gmail.com wrote: Hi, When I define a new class (through setClass), members defined in representation argument doesn't seem to like a class. For example, if

Re: [R] Defining class as a member of another class

Thank you Gabor. It works. I read a bit more and understood what you're doing. One more question, I want to know more about .Data property of a class. I know that it defines type of the class when asking typeof. But I don't know much how to use them. Would you mind pointing me to some reference

Re: [R] How to automatically placing a legend in empty space

I really apprecaite the suggestion. I just downloaded the Hmisc.pdf document and looked through it. Do you have any further thoughts about a specific function within the Hmisc package? I looked at largest.empty, but not sure how that is applicable to my data points. x = seq(0, 1000, by

[R] applying a function to a pair of components for each row of a list

Hi, I have a set of (x,y) coordinate pairs that are stored as a list my_list $x [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 $y [1] -8.0866819 -7.3876052 -6.6849311 -5.9837693 -5.2967432 -4.6525466 [7] -4.0999453 -3.6556190 -3.3076102 -3.0360780 -2.8220465

Re: [R] applying a function to a pair of components for each row of a list

Hi Kavitha, I must confess you that I not understood well what you are looking for. But.. mylist-list(x=1:25, y=runif(n=25)) plot(mylist, type=n) points(mylist, type=p, col=mylist$x) Hth, miltinho On Sun, Jun 28, 2009 at 1:03 AM, Kavitha Venkatesan kavitha.venkate...@gmail.com wrote: Hi,

Re: [R] applying a function to a pair of components for each row of a list

Milton, Thanks, the answer was in fact as simple as you pointed out. I was thinking more complicated than needed! Kavitha On Sun, Jun 28, 2009 at 1:11 AM, milton ruser milton.ru...@gmail.comwrote: Hi Kavitha, I must confess you that I not understood well what you are looking for. But..

[R] multiple regression w/ no intercept; strange results

I am writing some software to do multiple regression and am using r to benchmark the results. The results are squaring up nicely for the with-intercept case but not for the no-intercept case. I am not sure what R is doing to get the statistics for the 0 intercept case. For example, I would