[R] R for Windows 32-bit mode versus 64-bit mode

2013-02-06 Thread Roy Mendelssohn - NOAA Federal
Hi All: We are developing (actually updating) a piece of code that runs in R. Due to some limitations of some the libraries we use, at the present time we would be limited to the 32-bit version of R on Windows. Does anyone have a feel (or even real knowledge) if most people these days run

Re: [R] R for Windows 32-bit mode versus 64-bit mode

2013-02-07 Thread Roy Mendelssohn - NOAA Federal
Thanks. -Roy M. On Feb 6, 2013, at 11:55 PM, Prof Brian Ripley rip...@stats.ox.ac.uk wrote: On 06/02/2013 22:47, Roy Mendelssohn - NOAA Federal wrote: Hi All: We are developing (actually updating) a piece of code that runs in R. Due to some limitations of some the libraries we use

Re: [R] help with population matrix

2013-02-18 Thread Roy Mendelssohn - NOAA Federal
See if this helps at all: http://cran.r-project.org/web/packages/popbio/popbio.pdf -Roy M. On Feb 18, 2013, at 8:49 AM, John Kane jrkrid...@inbox.com wrote: Currently your code does not seem to make any sense ===comments in line John Kane Kingston ON Canada -Original

Re: [R] EEG data for time series

2013-03-10 Thread Roy Mendelssohn - NOAA Federal
Hi Erin: On Mar 9, 2013, at 7:36 PM, Erin Hodgess erinm.hodg...@gmail.com wrote: Dear R People: I have a data set with EEG data. There are 128 measurements per second for 16 locations. What is the best way to handle these series, please? Thanks, Erin I don't know about best, but

Re: [R] Classification of Multivariate Time Series

2013-05-27 Thread Roy Mendelssohn - NOAA Federal
Look at: State - Space Discrimination and Clustering of. Atmospheric Time Series Data. Based on Kullback Information Measures. Thomas Bengtsson If you Google the topic, there are host of other papers too, but the one meshes with exiting star-space methods. -Roy On May 27, 2013, at 4:34 AM,

Re: [R] mean

2013-08-30 Thread Roy Mendelssohn - NOAA Federal
Use is.numeric(list) on the list to see if it is a list of numbers or if it is list of characters. If it is a list of characters (which could be the case if you read it in using read.csv or the like) and it makes sense to convert to numeric, then do something like: list-as.numeric(list)

Re: [R] Having problems with doing anything with my data sets

2013-09-03 Thread Roy Mendelssohn - NOAA Federal
Hi: I know that it is frustrating when something like what you describe occurs, but realize that it is also frustrating for anyone on the list to try and help you, because there is no information to go on. What would help is information like: 1. Your OS and version 2. How you did the

Re: [R] Hdf files Download

2013-10-22 Thread Roy Mendelssohn - NOAA Federal
Your script didn't make it. There are restrictions on the mail list for attached files. You might try putting the script directly into the email. -Roy On Oct 22, 2013, at 9:40 AM, Tereza Smejkalova terka...@gmail.com wrote: I am triing one more time since my first message was rejected by

[R] State-space model with long tails

2012-12-12 Thread Roy Mendelssohn - NOAA Federal
Hi All: I have a dataset that when estimated as a dynamic state-space model, appears that the trend term would be better modeled with long-tails, such as a t-distribution, rather than Gaussian. I have looked at the manuals for the packages dlm, KFAS and sspir, and if I can estimate such a

Re: [R] funtion equivalent of jitter to move figures on device

2013-03-12 Thread Roy Mendelssohn - NOAA Federal
Hi Michael: I am not totally certain what you are trying to do, but Hadley Wickham and a student have been working on a package that allows a variety of plots to be added to maps. Examples can be found in this paper in Environmetrics : Glyph-maps for visually exploring temporal patterns in

Re: [R] EDCR build from source with namespace problem

2013-04-19 Thread Roy Mendelssohn - NOAA Federal
Several changes in R and in Java broke the code in EDCR. We have in beta a new version that fixes these problems. I will be certain to let you know when it is ready. -Roy M. On Apr 19, 2013, at 4:11 AM, ltdm luis.tito-de-mor...@ird.fr wrote: Hi R community, I am running R version

[R] ggmap (or ggplot2 with mapping question)

2014-12-26 Thread Roy Mendelssohn - NOAA Federal
HI All: Okay I am feeling particular brain dead today. i have the following map created so far using ggmap: library(ggmap) myLocation - c(-178, -28, -72, 48) myMap - get_map(location=myLocation, source=google, maptype=satellite, crop=FALSE) I want to color in some lat-lon 5 degree

Re: [R] ggmap (or ggplot2 with mapping question)

2014-12-26 Thread Roy Mendelssohn - NOAA Federal
Thanks. The map specs are stored in myMap, it can be displayed just by typing at the prompt: myMap Thanks for the geom.rect suggestion. I will look at that. -Roy On Dec 26, 2014, at 4:16 PM, David Winsemius dwinsem...@comcast.net wrote: On Dec 26, 2014, at 3:26 PM, Roy Mendelssohn

Re: [R] ggmap (or ggplot2 with mapping question)

2014-12-26 Thread Roy Mendelssohn - NOAA Federal
Winsemius dwinsem...@comcast.net wrote: On Dec 26, 2014, at 3:26 PM, Roy Mendelssohn - NOAA Federal roy.mendelss...@noaa.gov wrote: HI All: Okay I am feeling particular brain dead today. i have the following map created so far using ggmap: library(ggmap) myLocation - c(-178, -28, -72

Re: [R] Reading *.nc4 data

2015-01-26 Thread Roy Mendelssohn - NOAA Federal
I don’t think the ncdf library open netcdf4 files, which your file would seem to be. You need the library ncdf4 for that. For variety of reasons perhaps others can explain, ncdf4 for Windows is not on CRAN. You can download it from: http://cirrus.ucsd.edu/~pierce/ncdf -Roy On Jan 26, 2015,

Re: [R] Cube of Matrices or list of Matrices

2015-01-19 Thread Roy Mendelssohn - NOAA Federal
I believe what Karim is after is often referred to as a “ragged array”. For disk storage, such structures have been added to netcdf4 for things like subsurface profiles with a different number of depths. This blog might be of interest:

[R] ggplot2, geom_point with some symbols unfilled

2015-01-05 Thread Roy Mendelssohn - NOAA Federal
Hi All: I have a data frame the first part of which is the following: date lonlat mean chlcol missing 1 2003-04-23 203.899 19.664 0.0730 0.15176637 0 2 2003-04-24 204.151 19.821NaNNaN 1 3 2003-04-30 203.919 20.351 0.0740

Re: [R] ggplot2, geom_point with some symbols unfilled

2015-01-05 Thread Roy Mendelssohn - NOAA Federal
)) + geom_point(aes(colour=chlcol, shape=factor(missing)),size=3) + scale_colour_gradient(low=lightblue, high=red, na.value=black) + scale_shape_manual(values=c(19,1)) + theme_bw() Best, Ista On Mon, Jan 5, 2015 at 6:24 PM, Roy Mendelssohn - NOAA Federal roy.mendelss...@noaa.gov wrote: Hi

Re: [R] Conversion of Matlab code to an R code

2015-03-23 Thread Roy Mendelssohn - NOAA Federal
On Mar 23, 2015, at 8:10 AM, Abhinaba Roy abhinabaro...@gmail.com wrote: Hi, Can a Matlab code be converted to R code? I am finding it difficult to do so. Could you please help me out with it. Your help will be highly appreciated. snip If you mean something that can do an automatic

Re: [R] Lubridate and NameSpace

2015-02-28 Thread Roy Mendelssohn - NOAA Federal
Not certain but perhaps: lubridate::%m+% there is a subtle difference between what the two and three colons do. -Roy On Feb 28, 2015, at 9:25 AM, Glenn Schultz glennmschu...@me.com wrote: Hello All, I am working on a package very near completion - roxygenizing it now. The namespace

Re: [R] lag1.plot {astsa} vs. lag.plot {stats}

2015-04-20 Thread Roy Mendelssohn - NOAA Federal
snip By the way, the tsa3 issues page that you reference above...it's indicates the problems with existing time series functions as the reason for developing corrected functsion in astsa/tsa3. But the actual documentation for these corrected functions are extremely sparse. Is there

[R] GSVD and qr questions

2015-05-04 Thread Roy Mendelssohn - NOAA Federal
Hi All: I posted some of this earlier to the r-sig-mac because at that time I got a Mac specific error, but now I have a more general question. For a canonical variate analysis (CVA) used in some system identification routines, I need to calculate a generalized svd (GSVD). On the web there

Re: [R] GSVD and qr questions

2015-05-04 Thread Roy Mendelssohn - NOAA Federal
completion. Does this provide what you are looking for? Rich On Mon, May 4, 2015 at 11:56 PM, Roy Mendelssohn - NOAA Federal roy.mendelss...@noaa.gov wrote: Hi All: I posted some of this earlier to the r-sig-mac because at that time I got a Mac specific error, but now I have a more

[R] Changing layout in grid.arrange

2015-05-08 Thread Roy Mendelssohn - NOAA Federal
Hi All: I am doing something very similar to the the example in the grid.arrange package: require(ggplot2) plots = lapply(1:10, function(.x) qplot(1:10,rnorm(10), main=paste(plot,.x))) require(gridExtra) do.call(grid.arrange, plots) If you run this, the layout is 4 rows and 3 columns with

Re: [R] Changing layout in grid.arrange

2015-05-08 Thread Roy Mendelssohn - NOAA Federal
(mar=c(0,0,0,0)) plot(0.5,0.5,xlim=c(0,1),ylim=c(0,1), axes=FALSE,xlab=,ylab=,type=n) box() text(0.5,0.5,scr) } close.screen(all=TRUE) Jim On Sat, May 9, 2015 at 12:53 AM, Roy Mendelssohn - NOAA Federal roy.mendelss...@noaa.gov wrote: Hi All: I am doing something very similar

[R] xtractomatic package

2015-04-18 Thread Roy Mendelssohn - NOAA Federal
xtractomatic R package for accessing environmental data xtractomatic is an R package developed to subset and extract satellite and other oceanographic related data from a remote server. The program can extract data for a moving point in time along a user-supplied set of longitude, latitude and

Re: [R] lag1.plot {astsa} vs. lag.plot {stats}

2015-04-17 Thread Roy Mendelssohn - NOAA Federal
Not certain which plot you are looking at, but my guess is the answer is contained somewhere here: http://www.stat.pitt.edu/stoffer/tsa3/Rissues.htm in particular perhaps issues 4-5. -Roy M. On Apr 17, 2015, at 7:30 PM, Paul Domaskis paul.domas...@gmail.com wrote: I'm following

[R] ggpairs question

2015-06-10 Thread Roy Mendelssohn - NOAA Federal
Hi All: I am not totally grasping GGally, in particular how to customize ggpairs (or maybe I need to use a different part of GGally for what I want). Toy example: dput(junk) structure(list(state1 = c(-1.78772815343945, 0.347049897473099, -0.18529149576171, 1.79047951942822,

[R] Simulating a time series with a given spectrum

2015-05-27 Thread Roy Mendelssohn - NOAA Federal
Hi All: Is there a routine in R that allows me to simulate a time series that has a given spectrum? I have looked at the R Time Series Task view, and have done a web search also, including an article to appear in the handbook of statistics on time series in R, but I don’t see anything

Re: [R] Simulating a time series with a given spectrum

2015-05-27 Thread Roy Mendelssohn - NOAA Federal
Thanks, I will look at that. -Roy On May 27, 2015, at 9:15 AM, Prof Brian Ripley rip...@stats.ox.ac.uk wrote: On 27/05/2015 16:54, Roy Mendelssohn - NOAA Federal wrote: Hi All: Is there a routine in R that allows me to simulate a time series that has a given spectrum? I have looked

Re: [R] opendap accessibility on Windows

2015-07-06 Thread Roy Mendelssohn - NOAA Federal
Hi Christina: My memory is that Pierce had problems getting the netcdf4 library to work with OpeNDAP included on Windows. OPeNDAP does work with ncdf4 on a Mac or a Linux box. So one option is to run on a Linux VM on Windows. There may be more direct options, depending on what datasets you

Re: [R] About manipulating NetCDF files in ncdf package.

2015-05-20 Thread Roy Mendelssohn - NOAA Federal
Hi Thanh: You are confounding several issues here, and are providing only incomplete information about your data. The confounded issues are: 1. How to write a netcdf with a given coordinate system. 2. How to calculate a monthly average or monthly subsample from a larger sample. I would

Re: [R] Environmental Data Connector v1.3

2015-08-03 Thread Roy Mendelssohn - NOAA Federal
Hi Robert: I didn’t see this until Dan sent me something offline. I apologize for the problem. Yes the function should work as is, and a lot of the EDC is Java. I have forwarded your email to the people who did the coding. But in the meantime can you do two things for me to help us in the

Re: [R] Environmental Data Connector v1.3

2015-08-03 Thread Roy Mendelssohn - NOAA Federal
On Aug 3, 2015, at 8:12 AM, Robert in SA ri.william...@outlook.com wrote: Hello. I have successfully installed EDC v1.3 on linux ubuntu 14.04. I am running a 64bit machine with R 3.2.1 via Rstudio. I have tried example1 - EDC.get(1) after loading the ncdf and EDCR libraries, from both the

Re: [R] Maxent Jarfile

2015-07-09 Thread Roy Mendelssohn - NOAA Federal
Message From: Roy Mendelssohn - NOAA Federal Sent: Thursday 9 July 2015 18:48 To: annemariefische...@gmail.com Cc: David Winsemius; John Kane; R. Mailing List Subject: Re: [R] Maxent Jarfile ** The contents of this message do not reflect any position of the U.S. Government

Re: [R] Trellis Plots: translating lattice xyplot() to ggplot()

2015-07-10 Thread Roy Mendelssohn - NOAA Federal
Don’t know for certain but might this help: http://journal.r-project.org/archive/2015-1/murrell.pdf From the latest issue of R Journal. Abstract The gridGraphics package provides a function, grid.echo(), that can be used to convert a plot drawn with the graphics package to a visually

Re: [R] Maxent Jarfile

2015-07-09 Thread Roy Mendelssohn - NOAA Federal
Have you been able to run the example that comes with dismo::maxent such as at the bottom of: http://www.inside-r.org/packages/cran/dismo/docs/maxent -Roy On Jul 9, 2015, at 9:38 AM, annemariefische...@gmail.com wrote: Hi John, I have attempted to add links to the files to the thread

Re: [R] TSclust multivariate time series clustering

2015-08-26 Thread Roy Mendelssohn - NOAA Federal
I suggest you look at this paper: http://www.jstatsoft.org/v62/i01/paper -Roy On Aug 26, 2015, at 1:02 AM, GloriaSolano sol...@europrofessionals.org wrote: Hello, I am trying to cluster multivariate time series with the R package TSclust. I have a dataset of 45 companies with 10 years

Re: [R] system() command not working

2016-06-04 Thread Roy Mendelssohn - NOAA Federal
Hi John: When El Capitan first came out there was a discussion in the R-SIg-Mac list about environmental variables not being passed down to applications (not just R abut in general). I believe a work around was suggested, but I would search the archives for that. So what is happening,

Re: [R] Request for help

2016-06-03 Thread Roy Mendelssohn - NOAA Federal
Hi All: > On Jun 3, 2016, at 11:33 AM, jlu...@ria.buffalo.edu wrote: > > There is a video tutorial on the RStudio web site showing how to create R > packages within RStudio. Hadley Wickham also has a book on creating R > packages. > And I would add that Hadley has kindly put the book

Re: [R] How to get bathymetry data using R

2016-06-11 Thread Roy Mendelssohn - NOAA Federal
Use the “rerddap” package that accesses our ERDDAP server. You can see the data available at: http://upwell.pfeg.noaa.gov/erddap In the search box type in “bathymetry”.You can subset as you want, using “rerddap” you get back a netcdf file that is already read into your R workspace, or

Re: [R] R editor for Mac

2016-01-20 Thread Roy Mendelssohn - NOAA Federal
Hi: Both the default Mac OS X installation of R (with the GUI) as well as RStudio, have very nice editors for coding in R. -Roy > On Jan 20, 2016, at 10:22 AM, Christofer Bogaso > wrote: > > Hi, > > Could you please suggest a good R editor for Mac OS X

Re: [R] Reshaping an array - how does it work in R

2016-03-18 Thread Roy Mendelssohn - NOAA Federal
o. > > HTH, > Denes > > > > On 03/18/2016 10:28 PM, Roy Mendelssohn - NOAA Federal wrote: >> Hi All: >> >> I am working with a very large array. if noLat is the number of latitudes, >> noLon the number of longitudes and noTime the number of time p

Re: [R] Reshaping an array - how does it work in R

2016-03-18 Thread Roy Mendelssohn - NOAA Federal
> On Mar 18, 2016, at 2:56 PM, Bert Gunter wrote: > > However copying may occur anyway as part of R's semantics. Others will > have to help you on that, as the details here are beyond me. > > Cheers, > Bert Hi Bert: Thanks for your response. The only part I was

Re: [R] Reshaping an array - how does it work in R

2016-03-18 Thread Roy Mendelssohn - NOAA Federal
8, 2016 at 3:15 PM, Roy Mendelssohn - NOAA Federal > <roy.mendelss...@noaa.gov> wrote: >> Thanks. That is what I needed to know. I don’t want to play around with >> some of the other suggestions, as I don’t totally understand what they do, >> and don’t want

[R] Reshaping an array - how does it work in R

2016-03-19 Thread Roy Mendelssohn - NOAA Federal
Hi All: I am working with a very large array. if noLat is the number of latitudes, noLon the number of longitudes and noTime the number of time periods, the array is of the form: myData[noLat, no Lon, noTime]. It is read in this way because that is how it is stored in a (series) of netcdf

Re: [R] Reshaping an array - how does it work in R

2016-03-19 Thread Roy Mendelssohn - NOAA Federal
> On Mar 19, 2016, at 8:18 AM, Henrik Bengtsson <henrik.bengts...@gmail.com> > wrote: > > On Fri, Mar 18, 2016 at 8:28 PM, Roy Mendelssohn - NOAA Federal > <roy.mendelss...@noaa.gov> wrote: >> Hi Henrik: >> >> I want to do want in oceanography is ca

Re: [R] Memory usage in prcomp

2016-03-22 Thread Roy Mendelssohn - NOAA Federal
> On Mar 22, 2016, at 10:00 AM, Martin Maechler <maech...@stat.math.ethz.ch> > wrote: > >>>>>> Roy Mendelssohn <- NOAA Federal <roy.mendelss...@noaa.gov>> >>>>>>on Tue, 22 Mar 2016 07:42:10 -0700 writes: > >> Hi

Re: [R] Reshaping an array - how does it work in R

2016-03-22 Thread Roy Mendelssohn - NOAA Federal
uot;outside", >> notably Tomas Kalibera. >> >> And hence:*NO* such strange workarounds are needed in this specific case: >> >> > Workaround: use data.table::setattr or bit::setattr to modify the >> > dimensions in place (i.e., without mak

[R] Memory usage in prcomp

2016-03-22 Thread Roy Mendelssohn - NOAA Federal
Hi All: I am running prcomp on a very large array, roughly [50, 3650]. The array itself is 16GB. I am running on a Unix machine and am running “top” at the same time and am quite surprised to see that the application memory usage is 76GB. I have the “tol” set very high (.8) so that it

Re: [R] Reshaping an array - how does it work in R

2016-03-21 Thread Roy Mendelssohn - NOAA Federal
Thanks for the info, but I will stay with regular R. Work -arounds for what I want to do just took some thought and programming, I just didn’t know if R copied the array or just manipulated indices, and given the size of the array I am memory limited. This gets into the old thing of

Re: [R] Unexpected values obtained when reading in data using ncdf and ncdf4

2016-04-22 Thread Roy Mendelssohn - NOAA Federal
Hi Louise: If Dave can’t figure it out, I can give a look also. A couple of things I would suggest: 1. Don’t use the name “data” in the nc_open command, that is a reserved command in R and you never know what problems that can cause. 2. You are doing calculations to get set the start and

Re: [R] Matrix multiplications

2016-05-21 Thread Roy Mendelssohn - NOAA Federal
> str(t(y-X %*% b)) num [1, 1:10] 0.595 -1.7538 -0.0498 -1.651 -0.6328 ... > str((y-X %*% b)) num [1:10, 1] 0.595 -1.7538 -0.0498 -1.651 -0.6328 … -Roy > On May 21, 2016, at 12:00 PM, george brida wrote: > > Dear R users: > > I have written the following lines : >

Re: [R] about netcdf files

2016-07-26 Thread Roy Mendelssohn - NOAA Federal
count > = c(1, 1, -1)) where the -1 value tells it to get all the values along the time dimension. You now have that slice of the data stored in myPrecip. I hope this gets you started, but i strongly urge you to read through the link above. -Roy > > On Tue, Jul 26, 2016 at 2

Re: [R] about netcdf files

2016-07-26 Thread Roy Mendelssohn - NOAA Federal
Hi Lily: I doubt the mail-list would pass through the netcdf file. Instead, could you do the following, and post the results: library(ncdf4 pre1 = nc_open('sample_precip_daily.nc') str(pre1) nc_close(pre1) I have a feeling you haven't worked much with netcdf files. I will try to find a

Re: [R] C/C++/Fortran Rolling Window Regressions

2016-07-21 Thread Roy Mendelssohn - NOAA Federal
I have no idea which method produces the fastest results, but the package KFAS has a function to do recursive regressions using the Kalman filter. One difference is that it is not, as far as a I can telll, a moving window (so past data are being dropped), just a recursively computed

Re: [R] about netcdf files

2016-08-01 Thread Roy Mendelssohn - NOAA Federal
Hi Lily: If you download the vignette to my xtractomatic package (http://coastwatch.pfeg.noaa.gov/xtracto/index.html) there are any number of examples using ggplot2 to make maps from netcdf data, HTH, -Roy > On Aug 1, 2016, at 10:35 AM, lily li wrote: > > Hi all, > >

Re: [R] a package-building opinion question, please

2016-07-08 Thread Roy Mendelssohn - NOAA Federal
Hi Erin: Everyone's tastes differ, but when I started out knowing nothing about packages, I found Hadley's guide (thank you Hadley!!!) indispensable: http://r-pkgs.had.co.nz/intro.html For obvious reasons, his guide really is geared to built-in features of RStudio, though I believe

Re: [R] dplyr : row total for all groups in dplyr summarise

2016-07-06 Thread Roy Mendelssohn - NOAA Federal
Thanks muchly. I hate the smart quotes! -Roy > On Jul 6, 2016, at 11:42 AM, David Winsemius wrote: > > >> On Jul 6, 2016, at 9:45 AM, rmendelss gmail wrote: >> >> >>> On Jul 6, 2016, at 9:36 AM, David Winsemius wrote:

Re: [R] netcdf data precision or least significant digit

2016-07-07 Thread Roy Mendelssohn - NOAA Federal
Hi Ismail: Can you point me to a particular netcdf file you are working with. I would like to play with it for awhile. I am pretty certain the scale factor is 0.01 and what you are seeing in rounding error (or mor precisely I should say problems with representations of floating point

Re: [R] netcdf data precision or least significant digit

2016-07-07 Thread Roy Mendelssohn - NOAA Federal
ically, their absolute values must be equal. > > 1- http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis2.html > 2- > ftp://ftp.cdc.noaa.gov/Datasets/ncep.reanalysis2.dailyavgs/pressure/uwnd.2015.nc > > > >> On 08 Jul 2016, at 02:27, Roy Mendelssohn - NOAA Fed

Re: [R] netcdf data precision or least significant digit

2016-07-07 Thread Roy Mendelssohn - NOAA Federal
kage but author said that it is nothing > to do with ncdf4 package. > > Please, forgive me for taking your time. > > >> On 08 Jul 2016, at 03:21, Roy Mendelssohn - NOAA Federal >> <roy.mendelss...@noaa.gov> wrote: >> >> After looking at the file, doing an extrac

[R] Update to the xtractomatic package

2016-08-08 Thread Roy Mendelssohn - NOAA Federal
xtractomatic is an R package developed to subset and extract satellite and other oceanographic related data from a remote server. The program can extract data for a moving point in time along a user-supplied set of longitude, latitude and time points; in a 3D bounding box; or within a polygon

Re: [R] open a zip file

2016-08-07 Thread Roy Mendelssohn - NOAA Federal
If I break it into parts, I find that the "GET" fails. > Year <- format(Sys.Date(), "%Y") >Month <- format(Sys.Date(), "%m") > junk <- paste("https://mbsdisclosure.fanniemae.com/disclosure-docs/monthly/;, +"mbs", +as.character(Month), +

Re: [R] Extracting matrix from netCDF file using ncdf4 package

2016-07-02 Thread Roy Mendelssohn - NOAA Federal
Sending this to Hemant a second time as i forgot to reply to list. Hi Hemant: Well technically the code you give below shouldn’t work, because “start” and “count” are suppose to be of the same dimensions as the variables. I guess Pierce’s code must be very forgiving if that is working. One

[R] testers sought for xtractomatic-like package

2017-02-03 Thread Roy Mendelssohn - NOAA Federal
rerddapXtracto is an R package developed to subset and extract satellite and other oceanographic related data from a remote ERDDAP server. The program can extract data for a moving point in time along a user-supplied set of longitude, latitude and time points; in a 3D bounding box; or within a

[R] plotly example that highlights a line

2017-02-28 Thread Roy Mendelssohn - NOAA Federal
Hi All: In searching online, I have found examples of using plotly with ggplot2 graphics, say using geom_line, where there are multiple lines and by selecting the "factor" in the legend makes the particular line disappear or reappear (see https://plot.ly/ggplot2/). I am wondering if

[R] Variable String formation

2016-09-22 Thread Roy Mendelssohn - NOAA Federal
Hi All: I am trying to write code to create a string to be executed as a command. The string will be of the form: "param <- param[,rev(seq_len(dataYLen)),,drop = FALSE]" Now just creating that string is simple enough. Where the problem arises is the array param could be 2, 3, or 4

Re: [R] OPeNDAP access / OPeNDAP subsetting with R

2016-09-30 Thread Roy Mendelssohn - NOAA Federal
t;) > x<-nc_open(gribfile) > > > nc_open doesn't work. > > which command should I use? > > with regards > -Deb > > > On Tue, Sep 27, 2016 at 9:30 PM, Michael Sumner <mdsum...@gmail.com> wrote: > Opendap won't work on Windows CRAN build of ncdf4,

Re: [R] OPeNDAP access / OPeNDAP subsetting with R

2016-09-27 Thread Roy Mendelssohn - NOAA Federal
Look at the package ncdf4. You can use an OPeNDAP URL in place of the file name to perform subsets., -Roy > On Sep 27, 2016, at 9:06 AM, Debasish Pai Mazumder wrote: > > Hi all, > > I would like to access and subset following OpeNDAP files. > server: >

Re: [R] OPeNDAP access / OPeNDAP subsetting with R

2016-09-27 Thread Roy Mendelssohn - NOAA Federal
it > doesn't work. > the OPeNDAP link of the data > http://nomads.ncdc.noaa.gov/thredds/dodsC/modeldata/cfsv2_forecast_ts_9mon/2014/201404/20140403/2014040312/ > > datafile: > tmax.01.2014040312.daily.grb2 > > Thanks > -Deb > > On Tue, Sep 27, 2016 at 11:51

[R] xtractomatic package is now available in CRAN

2016-10-12 Thread Roy Mendelssohn - NOAA Federal
I am pleased to announce that the xtractomatic package is now available from CRAN. xtractomatic is an R package developed to subset and extract satellite and other oceanographic related data from a remote server. The program can extract data for a moving point in time along a user-supplied set

Re: [R] Error in R_nc4_open

2017-01-11 Thread Roy Mendelssohn - NOAA Federal
> On Jan 11, 2017, at 8:39 PM, Debasish Pai Mazumder wrote: > > Thanks so much Roy. It works. > Thanks Jeff for all your help. > As a part of NCAR Command Language help group, I was only concern about the > first response I received from this help group which will

Re: [R] Error in R_nc4_open

2017-01-11 Thread Roy Mendelssohn - NOAA Federal
Try replacing http with https. > gribfile <- > "https://nomads.ncdc.noaa.gov/thredds/dodsC/modeldata/cfsv2_forecast_ts_9mon/2011/201104/20110401/2011040100/tmax.01.2011040100.daily.grb2; > nc <- nc_open(gribfile) > str(nc) List of 14 $ filename : chr

Re: [R] Lubdridate: subset based on hour and minute

2017-01-02 Thread Roy Mendelssohn - NOAA Federal
> Perhaps you're showing me the way and I'm missing it - how would I > subset to only 1030 and 1100, excluding 1000? It seems I would need to > say, give me all time greater than 10:00, but the hours and minutes > are in separate slots, which is throwing me off. > > Thanks ag

Re: [R] Lubdridate: subset based on hour and minute

2017-01-02 Thread Roy Mendelssohn - NOAA Federal
Hi Joe: See below. > On Jan 2, 2017, at 12:01 PM, Joe Ceradini wrote: > > Hi folks, > > I must be missing something obvious/painfully simple here > > How do I subset a time vector based on hours AND minutes? So, in this > example, I want all time greater than 10:00,

Re: [R] Lubdridate: subset based on hour and minute

2017-01-02 Thread Roy Mendelssohn - NOAA Federal
somehow couldn't find an example like that via google. > > Joe > > On Mon, Jan 2, 2017 at 1:35 PM, Roy Mendelssohn - NOAA Federal > <roy.mendelss...@noaa.gov> wrote: >>> test > hm("10:00") >> [1] FALSE FALSE TRUE TRUE >>> test[test > hm(

[R] xtractomatic v3.2.0

2017-01-04 Thread Roy Mendelssohn - NOAA Federal
Hi All: I am pleased to announce that xtractomatic v3.2.0 is now available on CRAN. The changes in this version will be invisible to the user - the major changes are the use of https instead of http, and some changes in the vignette so that multiple attempts are made to download the data

Re: [R] Kalman filter for a time series

2017-07-30 Thread Roy Mendelssohn - NOAA Federal
> structSSM Is no longer part of KFAS. All you needed to do was: library(KFAS) ?KFAS and you would have seen that if you went to the index. A structural state space model is now built up from its components, much like in LM. Look at; ?SSModel -Roy > On Jul 29, 2017, at 9:26 PM, Staff

Re: [R] Accessing Pointers

2017-06-22 Thread Roy Mendelssohn - NOAA Federal
Hi Lawrence: > On Jun 22, 2017, at 4:26 PM, David Winsemius wrote: > > > >> is pointing to in the following line of code. Need some help. >> >> #install.packages('xml2') >> library('xml2') >> pg1 <- read_html("www.msn.com") > > Error: 'www.msn.com' does not exist in

[R] Xtractomatic version 3.3.2 now available on CRAN

2017-05-23 Thread Roy Mendelssohn - NOAA Federal
The xtractomatic package version 3.3.2 is now available on CRAN. Besides the improvements listed below, this release fixes a problem caused by an update to the Apache Tomcat used by the ERDDAP server, that broke an important function in the package. Many thanks to the CRAN maintainers for

Re: [R] Reversing one dimension of an array, in a generalized case

2017-06-01 Thread Roy Mendelssohn - NOAA Federal
From: Ismail SEZEN [mailto:sezenism...@gmail.com] > Sent: Thursday, June 1, 2017 3:35 PM > To: Roy Mendelssohn - NOAA Federal <roy.mendelss...@noaa.gov> > Cc: David L Carlson <dcarl...@tamu.edu>; R-help <r-help@r-project.org> > Subject: Re: [R] Reversing one dimensio

[R] Reversing one dimension of an array, in a generalized case

2017-06-01 Thread Roy Mendelssohn - NOAA Federal
Hi All: I have been looking for an elegant way to do the following, but haven't found it, I have never had a good understanding of any of the "apply" functions. A simplified idea is I have an array, say: junk(5, 10, 3) where (5, 10, 3) give the dimension sizes, and I want to reverse the

Re: [R] Reversing one dimension of an array, in a generalized case

2017-06-01 Thread Roy Mendelssohn - NOAA Federal
ep coming along > and sticking things into it." > -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) > > > On Thu, Jun 1, 2017 at 9:51 AM, Roy Mendelssohn - NOAA Federal > <roy.mendelss...@noaa.gov> wrote: >> Hi All: >> >> I h

Re: [R] Reversing one dimension of an array, in a generalized case

2017-06-01 Thread Roy Mendelssohn - NOAA Federal
ean relative difference: 0.5855162" > > David C > > -Original Message- > From: Bert Gunter [mailto:bgunter.4...@gmail.com] > Sent: Thursday, June 1, 2017 2:00 PM > To: David L Carlson <dcarl...@tamu.edu> > Cc: Roy Mendelssohn - NOAA Federal <roy.mendel

Re: [R] Reversing one dimension of an array, in a generalized case

2017-06-01 Thread Roy Mendelssohn - NOAA Federal
latitudes in the array is irrelevant, as long as the mapping is correct. -Roy > On Jun 1, 2017, at 1:35 PM, Ismail SEZEN <sezenism...@gmail.com> wrote: > >> >> On 1 Jun 2017, at 22:42, Roy Mendelssohn - NOAA Federal >> <roy.mendelss...@noaa.gov

[R] RMarkdown question

2017-08-29 Thread Roy Mendelssohn - NOAA Federal
Hi All: In creating a R Notebook I know that in the text I can link to a (sub) section by using the command: [Header 1](#anchor) and putting the appropriate anchor name at the appropriate header. But can the same be done for code chunks, if the code chunk is named? What I want to do

Re: [R] RMarkdown question

2017-08-29 Thread Roy Mendelssohn - NOAA Federal
hen you can jump to this > via a link like "see [this code chunk](#foo)". > > Regards, > Yihui > -- > https://yihui.name > > > On Tue, Aug 29, 2017 at 1:30 PM, Roy Mendelssohn - NOAA Federal > <roy.mendelss...@noaa.gov> wrote: >> Hi All: >>

Re: [R] withr::set_makevars

2017-09-09 Thread Roy Mendelssohn - NOAA Federal
gt; TIBCO Software > wdunlap tibco.com > > On Wed, Sep 6, 2017 at 3:41 PM, Roy Mendelssohn - NOAA Federal > <roy.mendelss...@noaa.gov> wrote: > Hi All; > > This problem has come about from trying to learn some of the review practices > recommend by rOpensci. One of them is

Re: [R] withr::set_makevars

2017-09-06 Thread Roy Mendelssohn - NOAA Federal
nction(e)conditionMessage(e})) > > Then run your test and see what file set_makevars is complaining about and > what in the file might cause trouble for set_makevars. > > > Bill Dunlap > TIBCO Software > wdunlap tibco.com <http://tibco.com/> >

[R] withr::set_makevars

2017-09-06 Thread Roy Mendelssohn - NOAA Federal
Hi All; This problem has come about from trying to learn some of the review practices recommend by rOpensci. One of them is to use the package goodpractice. After installing goodpractice, it kept failing on my own packages which are under development, and I was concerned something was funny

Re: [R] Extracting subset from netCDF file using lat/lon and converting into .csv in R

2017-08-28 Thread Roy Mendelssohn - NOAA Federal
Two questions: 1. Is the order of the dimensions shown what is shown if you look at str(ncin) - I mean shown at the end where it describes the variable and its dimensions? 2. Is you problem how to subset the netcdf file, how to write to the .csv file, or both? -Roy > On Aug 28, 2017, at

[R] ggridges help

2017-10-17 Thread Roy Mendelssohn - NOAA Federal
Hi All: I am just not understanding ggridges. The data I have are time series at different depths in the ocean. I want to make a joy plot of the time series by depth. If I was just doing a ggplot2 line plot I would be doing: ggplot(plotFrame, aes(x = time, y = cycle, group = depth)) +

Re: [R] ggridges help

2017-10-17 Thread Roy Mendelssohn - NOAA Federal
d your question is...? > ... and the code you tried that didn't work was? > > Bert > > > On Oct 17, 2017 12:22 PM, "Roy Mendelssohn - NOAA Federal" > <roy.mendelss...@noaa.gov> wrote: > Hi All: > > I am just not understanding ggridges. The data I have are tim

Re: [R] ggridges help

2017-10-17 Thread Roy Mendelssohn - NOAA Federal
line(fill="red", min_height=-0.25) > > > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > On Tue, Oct 17, 2017 at 12:43 PM, Roy Mendelssohn - NOAA Federal > <roy.mendelss...@noaa.gov> wrote: > I have tried: > > ggplot(plotFrame, aes(x = time, y = cycl

Re: [R] ggridges help

2017-10-17 Thread Roy Mendelssohn - NOAA Federal
lap <wdun...@tibco.com> wrote: > > > > Does the following work for you? > > > >ggplot2::ggplot(plotFrame, aes(x = time, y = depth, height = cycle, > > group = depth)) + ggridges::geom_ridgeline(fill="red", min_height=-0.25) > > > > &g

[R] dygraphs, multiple graphs and shiny

2017-10-18 Thread Roy Mendelssohn - NOAA Federal
Hi All: This is really getting into the weeds, but I am hoping someone will have a solution. I am trying to use dygrahs for R, within Shiny. The situation arises when I am combining a number of dygraphs into one plot. If I am just in an RNotebook, if you look at:

Re: [R] dygraphs, multiple graphs and shiny

2017-10-18 Thread Roy Mendelssohn - NOAA Federal
raphs::renderDygraph({ > res = list() > res[[1]] <- dygraph(lungDeaths[, 1], group = 'lungs') %>% > dyRangeSelector() > res[[2]] <- dygraph(lungDeaths[, 1], group = 'lungs') %>% > dyRangeSelector() >res <- htmltools::tagList(res) >

[R] httr::content without message

2018-01-02 Thread Roy Mendelssohn - NOAA Federal
Hi All: I am using httr to download files form a service, in this case a .csv file. When I use httr::content on the result, I get a message. Since this will be in a package. I want to suppress the message, but haven't figured out how to do so. The following should reproduce the result:

[R] Development versions of xtractomatic and rerddapXtracto

2017-12-22 Thread Roy Mendelssohn - NOAA Federal
If you are a user of the R package "xtractomatic", I have a new development version available, as well as a test version of the package "rerddapXtracto". The biggest changes are functions that can take the output of any of the data download functions and quickly map the data. These

  1   2   >