similar to: climatological standard deviation

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

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: >

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:

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

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 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

Dear R Users, I am using R on windows. how to use 'points' function to plot two curves with errbar I am doing like: x.val <- as.integer(names(co2mean)) errbar(x.val, co2mean, co2mean + co2sd, co2mean - co2sd, xaxt='n', col=1, xlab=NA,ylab=NA)# obs error bar lines(x.val, co2mean, col=1, lwd=2) errbar(x.val, co2tm3.month.mean, co2tm3.month.mean + co2sd.tm3, co2tm3.month.mean -

Hello R Users, I want to convert date (yr, mo, day, hr, min, sec) to decimal date, for example: Date (in six columns): yr mo dy hr min sec 1993 02 13 05 52 00 Converted to Decimal date : 1993.3542 How to write a small code in R so I can convert six column date to decimal date Many thanks, Yogesh -- Yogesh K. Tiwari (Dr.rer.nat), Scientist, Indian Institute

Dear R Users, I have a CO2 time series. I want to fit this series seasonal cycle and trend with fourth harmonic function, and then compute residuals. I am doing something like: file<-read.csv("co2data.csv") names(file) attach(file) fit<-lm(co2~1+time+I(time^2)+sin(2*pi*time)+cos(2*pi*time)+sin(4*pi*time)+cos(4*pi*time)+

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:

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

Hello R Users, 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: ######### Latitute: 89.5 167 0 0 0 0 0 182 0 0 0 0 0 0 0 0 0 0 0 0 200

Hi, I am using R on Windows XP. I am using 'arrows' funtion to plot the variance as error bar, BUT error bar goes only one side of the data point, I need to plot the error bar on both side of the data point (plot is attached), I am using following commands to plot, plot(file3$lat,file3$STotwoKm,pch=21,cex=2.5,ylim=c(-0.2,2.5),xlim=c(-50,50),xlab=NA,ylab=NA, col=1,

Dear Sir I have successfuly downloaded R package. Kindly let me know where i will get _/*user mannual*/_ for the same. awaiting a prompt reply with regards Pankaj @#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@#@ Pankaj Kumar Research Scholar Climatology & Hydrometeorology Divison Indian Institute Of Tropical Meteorology Homi Bhabha Road, Pune-411008 India Phone No. +91-20-5893600 Ext. 361

Hello everybody, Which one is a better choice 1. Gateway device with FXO <-> SIP ( example Addpac http://www.addpac.com/addpac_eng2/addpac_product_view_detail.php?class_id=19&item_id=59 ) 2. Digium (Wildcard TDM400P) 3. Sangoma (A200 Analog FXO/FXS) All i need is to put asterisk in place with 4-8 incomming lines (ordinary POTS ). With IVR, Voice mail and International Call via SIP.

Hi all , I'm running Debian Etch . I just finished configuring SAMBA as PDC to authenticate against LDAP server which works. The system in question uses default debian etch packages. As My Linix/unix accounts can authenticate against it. The LDAP works. I Used the default shipped smbldap-populate script to setup SAMBA. Everything seems to work as Anonymous User or as user

GNU ld version 2.20.51.0.2-5.36.el6 20100205 Copyright 2009 Free Software Foundation, Inc. This program is free software; you may redistribute it under the terms of the GNU General Public License version 3 or (at your option) a later version. This program has absolutely no warranty. From: Kostya Serebryany [mailto:kcc at google.com] Sent: Monday, August 19, 2013 5:36 AM To: Sergey Matveev Cc:

Hello R Users, How to make a variable equidistance with time i.e. how to interpolate a variable if it is not sampled at equal time interval. Many thanks, Regards, Yogesh [[alternative HTML version deleted]]

Hi, I want 12 ticks at axis 1 and want to write Jan-Dec on each. something like: axis(1, at=1:12, labels=c('J','F','M','A','M','J','J','A','S','O','N','D')) I could omit default ticks but now how to control ticks. plot(file$time, file$ch4*1000, ylim=c(1500,1700), xaxt='n', xlab= NA,

Hi Daniel, To my knowledge, there is no other implementation available for Steensgaard AA except poolalloc. I have tried Steensgaard AA under LLVM 2.7 poolalloc, but there are lots of bugs which need to be fixed. It will be a good idea to use previous implementation, but you will have to modify (a lot) it with respect to the latest version of LLVM poolalloc. Thanks, -Yogesh -------------- next