similar to: Thornthwaite Method in R

Hola eRReros. Os lo explico con un ejemplo: # Cargamos los paquetes y el ejemplo install.packages("SPEI") library("SPEI") data(wichita) # los primeros 12 datos head(wichita,12) # mi subset de los primeros 12 datos meu<-wichita[c(1:12),] meu # como veis los valores de TMED son iguales en ambos dataframes. # ahora viene el problema # calculamos

Hi Experts, I am newbie in R. Could you please share your idea to create a stack with no value or zero value? Nrow=1,Ncol=1, Ncell=1, I aim to work with looping, I want to make a first layer with zero values to add other layers further. Thanks. Best Regards, ..................Anup KhanalNorwegian Institute of science and Technology (NTNU)Trondheim, NorwayMob:(+47) 45174313

Good afternoon, I am a new in R. I have to work with large climate data.I am not able to read the netcdf file. Can anyone try with this file attached ? Best Regards, ..................Anup KhanalNorwegian Institute of science and Technology (NTNU)Trondheim, NorwayMob:(+47) 45174313

Hi all, I have gone through some basics of R in eclipse environment. What I need to work is more than basic level so I could not find out the solution by myself. could anyone help me ? My case, I have climate data (precipitation) in nc format with single layer for each time interval in a day. That means for January 1, 2013, I have 8 separate files each containing 3 hrs cumulative precipitation

Dear exports,I have created a hypsometric curve (area-elevation curve) for my watershed by using simple command hypsometric(X,main="Hypsometric Curve", xlab="Relative Area above Elevation, (a/A)", ylab="Relative Elevation, (h/H)", col="blue")It plots the hypsometric curve in "RGraphics window", My question is how can I export values which is used

Hi Experts, Could I ask one question ? I have plotted hypsographic curve for one watershed. I want to see the values of curve as I have to create elevation-volume curve from it. Could you please explain me how to see data that is used in curve ? Best Regards,..................Anup KhanalNorwegian Institute of science and Technology (NTNU)Trondheim, NorwayMob:(+47) 45174313

*Dear **Madam**/**Sir,* *What would happen if we irrigate the degraded land using the desalinated saline water in Minqin Oasis, Gansu Province, China? Can we prevent the disappearing of Minqin Qasis?* * * *Background:* *1. **The Minqin Oasis is surrounded by two deserts, Badain Jaran Desert and Tengger Desert. * *2. **Water is the root of ecological

I'd appreciate advice on the following: I've written an R function that uses 3 vectors (temperature, precipitation and potential evapotranspiration) for a given site, calculates a water budget (which implies few (<4) iterations), and, from thresults of this water budget, calculates a number of bioclimatic indices. Now I want to calculate these indices for a large number of sites

Dear R users, I'm glad to announce the new version of water package (0.5). As this is my first message to the list, I want to add that this package provides tools to estimate actual evapotranspiration from surface energy balance models. Right now you can run the well-know METRIC model using it. This model allows to estimate the energy balance using landsat (7 or 5) images and a weather

Dear R users, I'm glad to announce the new version of water package (0.5). As this is my first message to the list, I want to add that this package provides tools to estimate actual evapotranspiration from surface energy balance models. Right now you can run the well-know METRIC model using it. This model allows to estimate the energy balance using landsat (7 or 5) images and a weather

Hi System: Linux Open SuSE 10.2 I upgraded to R 2.6.0 and wanted to upgrade my packages as well. I tried to use the checkBuild=TRUE option in update.packages, it asked me if I wanted to update the package and I answered with yes, and finally ig tave me the error message unknown parameter "checkBuild" in function install.packages. My question: is the option still necessary or is

Dear list, I have to look examine hourly time - series and would like to plot variable section of them using plot.zoo. Hourly time series data which looks like this: YYYY MM DD HH P-uk P-kor P-SME EPOT EREA RO R1 R2 RGES S-SNO SI SSM SUZ SLZ 2003 1 1 1 0.385 0.456 0.021 0.000 0.000 0.000 0.013 0.223 0.235 0.01 0.38

Hi everybody. I'm looking for the way to increase the depth value in 0.01 for each index group. Easier to explain with this example: >my_data=read.table("clipboard", header=TRUE) Depth s_name index 3852 Site_1 144 3852 Site_1 144 3852 Site_1 144 3852 site_A 145 3852 site_A 145 3852 site_A 145 3852 site_A 145 3852 site_B 147 3852 site_B 147 3852 site_B 147 3852 site_B 147 54962

Dear list members, I am new member and fairly new into R world! I hope what I have is not beyond the purpose of this list. I did first search for similar experimental designs without success. I want to perform an ANOVA analysis using the aov() function. I am not 100% sure that I have it right. If anyone can help me, that will be greatly appreciated. My design is not balanced for any of the

Dear all, I am trying to transform polynomial coefficients from orthogonal form to the standard power basis. There's poly.transform in S-plus. Does anybody know how to do that in R ? I've found question about that in the archives of R-help but no real answer. Example : I'm doing polynomial regression of percentage of one insect in a community on altitude, precipitations,

marked by the genius, the courage, the originality, and the written. Not only is Tradition silent, but even Fable herself does not attempt to tell the story of the origin of her population. Here stand In accordance with these principles, we observe that, while the most tear her away, and then, maddened by her resistance and the sight of the sea and of the land by evaporation, falls again, under

marked by the genius, the courage, the originality, and the written. Not only is Tradition silent, but even Fable herself does not attempt to tell the story of the origin of her population. Here stand In accordance with these principles, we observe that, while the most tear her away, and then, maddened by her resistance and the sight of the sea and of the land by evaporation, falls again, under

Hello there, I’m computing the total value of an order from the price of the order items using a “for” loop and the “ifelse” function. I do this on a large dataframe (close to 1m lines). The computation of this function is painfully slow: in 1min only about 90 rows are calculated. The computation time taken for a given number of rows increases with the size of the dataset, see the example with

Dear Dennis, Your question assumes that people know both SAS PROC MIXED and R nlme. Only a limited number of people do. Add the mathematical formulation of the model. That will increase the number of people that can help you. Adding the number of levels in each categorical variable and the number of observation per group is useful too. Best regards, ir. Thierry Onkelinx Instituut voor natuur-

I am trying to reproduce some old SAS PROC MIXED code using R and nlme. The data consists of emission readings from vehicles and fuel properties. All variables are real numbers except "study" and "vehicle", which are character. Unfortunately, since the data are confidential, I'm unable to provide an example. The original SAS v6.12 code is provided below: