Re: [R] Efficient way to convert covariance to Euclidian distance matrix

2013-11-02 Thread Jari Oksanen
Rolf Turner r.turner at auckland.ac.nz writes: On 10/31/13 23:14, Takatsugu Kobayashi wrote: I am struggling to come up with an efficient vectorized way to convert 20Kx20K covariance matrix to a Euclidian distance matrix as a surrogate for dissimilarity matrix. Hopefully I can apply

[R] ascii-grid export

2013-11-02 Thread Enzo Cocca
hi, I want to export in ascii-grid file a semivariogram map that I have did with gstat library. How can I make it? Somebody can help me? thanks! enzo -- Enzo Cocca (PhD Candidate) Research Fellow Università di Napoli L'Orientale mail: enzo@gmail.com cell: +393495087014

Re: [R] mapping data to a geographic map of Europe

2013-11-02 Thread Markus Kainu
Dear Claudia, I have written a tutorial for this task here; http://www.markuskainu.fi/r-tutorial/eurostat/spatial.html It uses the shapefile provided by Gisco and plots data from Eurostat. Packages rgdal and maptools are used to process the data and SpatialPolygonDataFrame is fortified into

Re: [R] ascii-grid export

2013-11-02 Thread Gergely Daróczi
If I understand what you are up to, my pander package might help. You mean something like this? I used the example form the package manual. library(pander) data(meuse) coordinates(meuse) = ~x+y pander(variogram(log(zinc)~1, meuse)) np dist

Re: [R] R Packages in Mac

2013-11-02 Thread Jim Silverton
Hello everyone, I am interested in finding out why some packages don't work on Mac like the package 'vcd' Is there a fix? -- Thanks, Jim. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] Path Analysis

2013-11-02 Thread Sarah Rogers
Hello, I have just started to work on a path analysis (see attached image for the diagram), but I have encountered an error message. This is the code I have used: cov_matrix-var(xdata) library(sem) model.xdata-specifyModel() x1 - y2, xy12, NA x2 - y1, xy21, NA y1 - y2, yy12, NA y2 - y3,

[R] how to compute correlation between times (with buffer before and after times)?

2013-11-02 Thread Yla Savh
Dear forum members:   I have two series of times to correlate. Additionally, I have to allow for a buffer before and after the times so that they can be within 20 seconds of each other and not seem discordant. Example of my data: Time 1 Time 2 9:26 9:26 10:31 10:29 9:22  9:35 9:35 9:38

Re: [R] R Packages in Mac

2013-11-02 Thread Prof Brian Ripley
On 02/11/2013 17:03, Jim Silverton wrote: Hello everyone, I am interested in finding out why some packages don't work on Mac like the package 'vcd' Is there a fix? What does don't work mean? That one does: see http://cran.r-project.org/web/packages/vcd/index.html and the results page.

Re: [R] R Packages in Mac

2013-11-02 Thread Ista Zahn
The vcd package works for me on a mac. Please be more specific, describe what you did, and what happened that you found surprising. Also include the output of sessionInfo() Best, Ista On Sat, Nov 2, 2013 at 1:03 PM, Jim Silverton jim.silver...@gmail.com wrote: Hello everyone, I am interested

Re: [R] Path Analysis

2013-11-02 Thread John Fox
Dear Sarah, It's generally a good idea to include a reproducible example if you want to get help with a problem, but in this case it's a safe bet that the problem is that the model you specified has no variance or covariance parameters for the variables x1 and x2, which, I assume, you mean to

[R] How to obtain nonparametric baseline hazard estimates in the gamma frailty model?

2013-11-02 Thread Y
I can easily get the other parameter estimates by using coxph() but don't know how to get the baseline hazard estimates from it. Does anyone know how to obtain nonparametric baseline hazard estimates in the gamma frailty model? Thanks, YH [[alternative HTML version deleted]]

Re: [R] (gam) formula: Why different results for terms being factor vs. numeric?

2013-11-02 Thread Marius Hofert
Dear Bert, Thanks for helping. Your questions 'answers' why I get the expected behavior if 'group' is a factor. My question was why I don't get the expected behavior if 'group' is not a factor. From a theoretical (non-programming) point of view, there is no difference in a factor with two

Re: [R] ascii-grid export

2013-11-02 Thread Barry Rowlingson
Or do you mean you want to write the gridded output of an interpolation you did (perhaps using kriging) in gstat as an ESRI ASCII Grid file for reading into a GIS? If so, you can probably do it with writeGDAL from the rgdal package, or writeRaster from the raster package. I don't really know

Re: [R] Installing RCurl: 'configure' exists but is not executable

2013-11-02 Thread Michael Hannon
That seems OK. Are you trying to install the RCurl package from source? IIRC, the configure script is commonly used in GNU-style, compile-from-source installations. -- Mike On Fri, Nov 1, 2013 at 8:38 PM, Rainer Schuermann rainer.schuerm...@gmx.net wrote: # ls -l `which curl-config`

Re: [R] Package(s) for making waffle plot-like figures?

2013-11-02 Thread Jim Lemon
On 11/02/2013 10:35 AM, Zhao Jin wrote: Dear all, I am trying to make a series of waffle plot-like figures for my data to visualize the ratios of amino acid residues at each position. For each one of 37 positions, there may be one to four different amino acid residues. So the data consist of

Re: [R] plot time series data in wide format

2013-11-02 Thread Jim Lemon
On 11/02/2013 11:08 AM, Gary Dong wrote: Dear R users, I wonder if there is a way that I can plot a time series data which is in a wide format like this: CITY_NAME 2000Q12000Q2 2000Q32000Q4 2001Q1 2001Q2 2001Q3 2001Q4 2002Q1 2002Q2 CITY1