similar to: problem selecting 330 to 30 longitude

i use library(ncdf) to read this file as follow library(ncdf) sst.nc = open.ncdf(title) lonall = get.var.ncdf(sst.nc,'lon') latall = get.var.ncdf(sst.nc,'lat') precip = get.var.ncdf(sst.nc,'pre') close(sst.nc) if i use this method my pc freeze and not respond until i restart it is there -- View this message in context:

Dear RUsers, I guess this is an easy question for someone a little familiar with programming...(which I am not)... I've got 2 colummns, one shows just dates(SST_date, Class 'Date' num), the other one shows the SeaSurfaceTemperature (SST, num) at that certain date. SST_date SST 2008-01-01 22.2 2008-01-02 21.8 2008-01-03 22.8 2008-01-04 22.9 2008-01-05 23.1 2008-01-06 23.2 ...

Hi, I am fitting a weighted least square regression and trying to compute SSE,SST and SSReg but I am not getting SST = SSReg + SSE and I dont know what I am coding wrong. Can you help please? xnam <-colnames(X) # colnames Design Matrix fmla1 <- as.formula(paste("Y ~",paste(xnam, collapse=

Dear R Users, I was hoping for some help with a recurrent error message in lmer. I am trying to model the effect of temperature on metabolic rate in animals (response = int.length) at different temperatures (mean.sst), with repeated measurements on the same individuals (random effect = female). Ideally I would make a random slope and intercept model where the rate can change differently with

Hello I am a new user of R . and I ve problem with R and netcdf . I succed installation , I could use all examples . But when I take my netcf it is different . I want to do statistic on this kind of file . 1) first calculate mean . my data is like that through ncdump -h test.nc netcdf test { dimensions: lat = 301 ; lon = 401 ; time = UNLIMITED ; // (80 currently)

Hi- I'm wrestling with some of my data apparently not being called into a GLM or an LM. I'm looking at factors affecting fish annual catch rates (ie. CPUE) over 30 years. Two of the factors I'm using are sea surface temperature and sea surface temperature anomaly. A small sample of my data is below: CPUE Year Vessel_ID Base_Port Boat_Lgth Planing SST Anomaly 0.127

Hi, I'd like to compare SST data for year 2000 with SST for year 2001. I managed to get filled contour plots showing monthly SST for both years (12 plots for each year, 24 plots in total). In order to compare year 2000 and year 2001, however, I'd need to have the same legend scale (same color bar) in both years plots. Thus, my question is: is there a way to set the legend scale in

Hi everyone ! So here's my question of the day ! I need to make a decision on whether or not to go to a voip solution or configure an existing pbx (norstar) that my company has available. We are a small startup. I'm wanting a solution that will support up to about 200 people, with direct dial-in capability, up to about 30 concurrent phone calls and good voice quality. Right now I have

Hi there - I'm having some echo problems on my snom 320 phones. Anybody experience this before ? I don't have any issues with the sipura 841s I have though. Any help is greatly appreciated. Thanks ! Nora Lavelle -------------- next part -------------- An HTML attachment was scrubbed... URL:

Sorry for all the newbie questions. I really appreciate everyone's help today. Okay I've got outgoing and incoming calls working with no echo. yay! Now I'm having an issue with SIP extension to extension calling. Any time I dial another extension it goes right into voice mail. My extensions.conf is pretty small and rough but, here's what I have right now. Most of it was taken

Hi, I am looking at the effects of two explanatory variables on chlorophyll. The data are an annual time-series (so are autocorrelated) and the relationships are non-linear. I want to account for autocorrelation in my model. The model I am trying to use is this: Library(mgcv) gam1 <-gam(Chl~s(wintersecchi)+s(SST),family=gaussian, na.action=na.omit, correlation=corAR1(form =~

Hi, I have been working with R for the last year and using the UKFSST package to look at satellite tag track data and SST information. Fpr those not familiar, the package uses the positions estimated by the satellite tags themselves and the associated SST data from servers (in this case, from the University of Hawaii and the NOAA Coastwatch website) for the time preiod to estimate a reasonable

Okay everyone - I'm moving away from using sipura 841 phones. I'm starting to test with Polycom IP 501 phones. We plan to upgrade our server to a dual t1 but, for now we have a digium tdm400P with 4 analog lines coming into it. So my question is will upgrading the IP phones with my existing digium tdm400 card be enough to satisfy my users ? or is it really a combo deal needing to

Dear colleagues, I am trying a glm.nb for the distribution of a plant species with 93 environmental variables. I execute the instruction and I get the following message: "Error: C stack usage is too close to the limit". How can I increase the memory of R? Your sincerely, Nora. _________________________________________________________________ [[alternative HTML version deleted]]

Hi, I have a problem with the appearance of legend in my plots. If I specify the legend positions in characters like "topright"..etc, it appears, if i specify it in terms of coordinates like "-1, 1" .. etc, it does not appear. Can anyone help me with this? script - x.date <- as.Date(paste(year, month, day, sep="-")) ts1.n.e3 <- ts(data.nemr.e3[,3]) z1.n.e3

Dear colleagues in R, Has anybody implemented the 1) (Goodman & Kruskal) lambda or the 2) (Thiel's) uncertainty coefficient statistics (in the asymmetric and symmetric forms), or is anyone aware that they might reside in some package? A search in the R archives does indicate that they are (somehow) part of the CoCo package, but I would rather not start learning how to transform my

Hello All, I am looking for some help with was is probably very simple: I want the y-axis title '$B&D(B13C$B"s(B VPDB' and similarly for nitrogen on a secondary axis. I am able to get the delta notation and per mil notation separately but one messes up the other with the inclusion of quotations: ylab=expression("{delta}^13*C \211 VPDB") gives the per mil just fine but

Hi everyone, My background is not math and I am trying to figure out exactly what equation to use to map a response variable in GIS based on the coefficients obtained from the GLM and the values of the independent variables in each grid cell of my study area. Most specifically, I want to know how to incorporate the Poisson family and log link function in the equation. I would really appreciate if

Hi, I seem to be having a problem adding the following two raster objects together - one is a rasterLayer, the other is a rasterBrick. The extent, resolution, and origin are the same, so according to my understand it should work. The objects look like so: > obs.clim class : RasterLayer dimensions : 60, 200, 12000 (nrow, ncol, ncell) resolution : 0.5, 0.5 (x, y) extent : -70,

Hi, I've extracted netcdf data in R, and am trying to output it in a matrix format. I can make the matrix, and I have the data, but I can't get the matrix to include the data... I'm not very experienced at using R, so please forgive the crudeness of these codes... # For the matrix wind_dir_data<-matrix( NA, 17, 365 , byrow = FALSE ) wind_speed_data<-matrix( NA, 17, 365,