similar to: harmonic function fiting? how to do

Dear R Users, If a variable, say CO2(ppm), is varying with time. Then how to calculate CO2 (ppm) growth rate /a-1 I have CO2 time series (1991-2000), as: time, year, month, day, hour, min, sec, lat, long, height, CO2 1991.476722 1991 6 24 0 5 0 -38.93 145.15 4270 353.680 1991.476741 1991 6 24 0 15 0 -39.20 145.22 4270 353.950 1991.476747 1991 6 24 0 18 0 -39.43 145.28 4270 353.510

Dear 'R' Users, I have CO2 timeseries which I want to fit with fourth harmonic function and then would like to compute residuals. I tried with two option but not sure which one is correct, kindly any one can help correcting me: # #I would like to fit a model with Fourth harmonic function as:

Hello R Users, I am using R on Windows, I am ploting CO2 (350,380) variying with year(1993,2003) I want to over-plot rainfall (10,250) varying with year(1993,2003) on axis-4 axis-1=year(1993,2003), axis-2=CO2(350,380) axis-3=None axis-4=rainfall(10,250) Kindly help how to over-plot another variable on axis-4. Many thanks, Regards, Yogesh [[alternative HTML version deleted]]

Dear R Users, I have a time series of CO2 (1993-2002) x-axis as time y-axis as CO2 I would like to over plot a CO2 average line at y-axis on the time series plot. How to do it. Kindly help, Thanks, Regards, Yogesh [[alternative HTML version deleted]]

Hello, When I plot any data with simple plot command in R, for example :- plot(time,co2,ylim=c(350,380),xlim=c(1993,2003),xlab=NA,ylab=NA,type="p",col=5) Then the first value of x-axis(350) and y-axis(1993) never starts from origin, always they sifted from the origin. Is there any command that I can correct this in the ploted figure and both the axis values start from origin.

Hello 'R' Users, I have a monthly mean CO2 necdf data file defined on 1x1 lat by lon coordinate. I want to sub-sample this variable CO2 on the coordinates of another ASCII data file. The coordinates of another ASCII data file are as: -24.01 152.06 -18.58 150.19 -13.46 148.35 -8.29 147.03 -3.14 146.19 1.53 145.59 7.08 145.33 12.25 145.02 17.46 144.31 22.44 142.35 27.53 141.26 33.04

Dear R users, I am comparing two data sets (CO2 observation vs. CO2 simulation, during 1993-2002). In order to do it I am calculating Root-Mean-Square(RMS) difference with following formula: > sqrt(sum((observed_residual - simulated_residual)^2)/n) # 'n' is number of observations Residuals are computed by fitting a harmonic function on both the data:

Good day: I am trying to use readcsvIts("nwr_data_qc.txt",informat=its.format("%Y%m%d%h%M %Y"),header=TRUE,sep="",skip=0,row.names=NULL,as.is=TRUE,dec=".") to read in a file (nwr_data_qc.txt) that looks like this: Time Y M D H Min CO2 2000.18790 2000. 3. 9. 18. 30. 373.60 2000.20156 2000. 3. 14. 18. 30. 373.34

Hello R Users, I am new to R and I have simple problem for R users. I have CO2 observations defined on time axis(yr,mo,day,hr,min,sec). (DATA ATTACHED HERE) First I want to convert time axis as one axis as 'hour' on regular interval as 1 hour. Say 00 hrs to 24hrs(jan1), 25hrs to 48hrs(jan2) and so on. Then I want to interpolate CO2 at every hour. Kindly anybody can help, Many thanks,

Dear R Users, I have CO2 mixing ratio(ppm) time series data during 1991-2000. I would like to calculate CO2 growth rate /ppm.a-1 (is it derivative dt/dco2 ??) Kindly can any one advise how to calculate above. My data file looks like; time, year, month, day, hour, min, sec, lat, long, height, CO2 1991.476722 1991 6 24 0 5 0 -38.93 145.15 4270 353.680 1991.476741 1991 6 24 0 15 0 -39.20 145.22

Sorry for re-posting the question, I did not get any reply. Kindly reply please if any one can. ### Hello R users, I have a montly time series over a several year period. It's easy to compute a monthly climatology (12 values), Now, I would like to calculate the corresponding standard deviation, ie the 12 values calculated from the january values, february values, etc. What's the best way

Dear R Users, I am using R 2.14.1 on windows How to read netcdf files in R ? Which packeges do we need to install for this, and what commands are used for reading netcdf files. Thanks, Best Regards, Yogesh Tiwari -- Yogesh K. Tiwari (Dr.rer.nat), Scientist, Centre for Climate Change Research, Indian Institute of Tropical Meteorology, Homi Bhabha Road, Pashan, Pune-411008 INDIA Phone:

Dear R Users, I did not get any reply on my question so I am re-asking. This time I am giving sample data: 1 60.3162 -13.5993 -0.4353 46.0938 0.1877 -0.194E-07 2 60.3713 -13.5992 -0.4423 46.1241 0.2057 -0.231E-06 3 60.3430 -13.5981 -1.6163 44.9048 0.2237 -0.270E-06 4 60.3227 -13.5970 -2.6258 43.8785 0.2213

Hello R users, I have a montly time series over a several year period. It's easy to compute a monthly climatology (12 values), Now, I would like to calculate the corresponding standard deviation, ie the 12 values calculated from the january values, february values, etc. What's the best way for such a calculation ? Regards, Yogesh -- Yogesh K. Tiwari (Dr.rer.nat), Scientist, Indian

Hello R Users, I am new to R. I have two data sets i) TOMS aerosol optical depth(AOD) and ii) TOMS ozone(O3). > > AOD data is on 1x1 grid and O3 data is on 5x5 grid. > > First I want to read AOD and O3 as it is and then I want to regrid AOD on > 5x5 grid as O3. > > Reading is first problem. > > FIRST PROBLEM READING AOD: > > AOD data is in following format: >

Hello R Users, How to convert decimal date to date as YYYY.mm.dd.hr.min.sec For example, I have decimal date in one column , and want to convert and write it in equivalent date(YYYY.mm.dd.hr.min.sec) in another next six columns. 1979.000000 1979.020833 1979.041667 1979.062500 Is it possible in R ? Kindly help, Regards, Yogesh -- Dr. Yogesh K. Tiwari, Scientist, Indian Institute

Hi r-help forum, I have been collecting contact data (with proximity logger collars) between a few different species of animal. All animals wear the collars, and any contact between the animals should be detected and recorded by both collars. However, this isn't always the case and more contacts may be recorded on one collar of the two. This is fine, it depends on battery life and other

I am having trouble fitting correlation structures within nlme. I would like to fit corCAR1, corGaus and corExp correlation structures to my data. I either get the error "step halving reduced below minimum in pnls step" or alternatively R crashes. My dataset is similar to the CO2 example in the nlme package. The one major difference is that in my case the 'conc' steps are

Dear all, When analyzing data from a climate change experiment using linear mixed-effects models, I recently came across a situation where: - the summary(model) showed a significant difference between the levels of a two-level factor, - while the anova(model) showed no significance for that factor (see below). My question now is: Is the anova.lme() approach correct for that model? And why does

I need some help understanding how on of the example data sets is formatted in the basic R installation. If I load the Mona Loa CO2 data, with the command: > data(co2) I can view the data with: > co2 And the data are in the form of 11 rows labeled as years (1994-2004) and 12 columns labeled (Jan - Dec). This structure appears to be a dataframe, however, if I type the command