similar to: paquete SPEI función thornthwaite

Hi,Has anybody worked with thornthwaite method to calculate the daily evapotranspiration ? I have used the following code to calculate monthly evaporation but what I need is daily evapotranspiration. Epot <- thornthwaite(Tave, lat, na.rm = FALSE) # Evaluates evapotranspiration for each month. Do you have any code which works for daily timestep ? Its more technical question, may be some

Dear R users, I am trying to produce some plots in a postscript file, but I am experiencing some issues. I open the device with ----------------------------- setPS() postscript (file='gs_mcmc_dust.ps',width=5*3,height=5*3,horizontal = FALSE, paper = "special",family = "ComputerModern",encoding="TeXtext.enc")#, onefile = FALSE -----------------------------

Dear list, This is a question about overdispersion and the ML estimates of the parameters returned by the glm.poisson.disp (L. Scrucca) and glm.nb (Venables and Ripley) functions. Both appear to assume a negative binomial distribution for the response variable. Paul and Banerjee (1998) developed C(alpha) tests for "interaction and main effects, in an unbalanced two-way layout of counts

Hello I have configured my samba server at my RedHat AS3 and I use nss_ldap to access to my W2k domain. Everything works fine. Now i have upgraded my server to RedHat EL4, and all my Windows users when the username have Uppercases in username can't acces to my samba shares, and all other users works. I modify my log level to 10 and when any user try access to samba server the log

Hola eRReros, Necesito extraer datos de un dataframe según uno de los factores. La cosa sería algo así: df.nuevo <- iris[iris$Species=="setosa"|iris$Species=="virginica",] Que me crea un df nuevo solo con las filas en que Species es setosa o virginica. El problema es que mi factor de interés tiene 22 niveles y a veces he de seleccionar 13, 10, 8 niveles o números así.

Hola a todos, Estoy trabajando con modelos lineales generalizados (GLM), en particular con las funciones glm y anova.glm de R. Tengo una pregunta que es más bien técnica en el sentido estadístico (y un poco Off Topic y probablemente un tanto naive). No tengo claro si es correcto decir análisis de varianza o de devianza al utilizar GLMs, veo que en el help de anova.glm dice análisis de

Hola a todos. Quería saber si conocéis alguna librería que haga esto: SPSS permite enlazarse con PowerPoint de manera que si quiero modificar un gráfico del ppt, con solo pinchar en el gráfico me aparece los datos del spss y al modificar estos cambia la gráfica directamente sin tener que hacer la gráfica de nuevo. Hay algo parecido en R???? Saludos José Miguel Contreras [[alternative

Hola eRReros, Estoy trabajando con datos geográficos unos en formato raster y otro en formato shape. Lo que quiero es hacer las estadísticas (sobretodo la media) de los píxeles del raster que caen dentro de algunos de los polígonos del shape. Lo he intentado con sextante pero se flipa con los valores perdidos, nuestros queridos NA. Y he pensado: estando el paquete raster en R seguro que se puede

Estimados Estoy realizando algunos cambios, migrando desde windows 7 64 bit a linux SUSE 12.02. No encontré inconvenientes hasta que: Error en plot.new() : márgenes de figura muy grandes Este error se da en Rstudio, ¿alguna idea? En windows sabía que usaba la consola en varias ventanas, pero en linux ¿...?. ¿Hay alguna configuración que desconozco? Desde ya muchas gracias. Javier Marcuzzi

I'm setting up a Samba/LDAP PDC (samba-3.0.13 / openldap-2.1.30) and I'm trying to join a machine (called "melisa") to the domain. When I try to join the domain, I type in the Administrator (which maps to root) credentials. and after a few seconds I get a windows error that says: The following error occurred while attempting to joing the domain "PCA-USERS":

Hello, I'm trying to read a file that has ascci data in it. The format of the data is in this order: test_number(integer) 4 double readings. Simple,ok, but for life of me I can't get R to put it into the internal format I need. What I want to do is to make an array called dat and dim it 10000 rows by 5 columns. That seems simple enough, I start out by: dat <-

I found the need of converting a matrix into a list of its columns (for use with do.call), and was surprised there was no method as.list.matrix, it could easily be a part of as.list default I wrote my.as.list.matrix <- function(mat) { if(!is.matrix(mat))stop("Argument must be a matrix") n <- NCOL(mat) res <- vector(mode="list", length=n)

Buenos dìas compañeros. Estoy trabajando con la funciòn FMA del paquete FactoMiner, pero no me permite ver las salidas gràficas de la funciòn con el siguiente mensaje: Error en plot.new() : márgenes de figura muy grandes Mis variables continuas estàn estandarizadas y no he podido conguir resultados. Agradezco cualquier ayuda. Cordialmente, -- Luis Alberto López González [[alternative HTML

Hola a todos. Un problema dos pasos. Quiero desde RStudio trabajar con dos versiones de R (pongamos R2.13 y R2.15) desde Linux, Ubuntu 10.04 (lucid). Ya tengo instalada la R.2.15.1. Por defecto trabajo con ella. Pasos que creo he de dar: 1º. Instalar dos versiones de R en la misma maquina. ¿Como se hace?. Me estoy metiendo con el PATH y tal... hay algún protocolo clarito por ahí? 2º.

Hola, ¿qué tal? Llevo un rato dándole vueltas y no encuentro la solución , seguro que es una tontería pero no acabo de verlo. Tengo un data frame de la siguiente forma. TARJETA FECHA 1 2010006207318 2011-02-01 2 2010006207318 2011-02-01 3 2010006209657 2011-02-01 4 2010006209657 2011-02-01 5 2010006209657 2011-02-01 6 2010006209657 2011-02-01 7

So I have 2 sets of data - a training data set and a test data set. I've been doing the analysis on the training data set and then using predict and feeding the test data through that. There are 114 rows in the training data and 117 in the test data and 1024 columns in both. It's actually the same set of data split into two. The rows are made of 5 different numbers. They do represent

Hi, I am a new user of R using R 2.8.1 in windows 2003. I have a csv file with single column which contain the 30,000 students names. There were typo errors while entering this student names. The actual list of names is < 1000. However we dont have that list for keyword search. I am interested in grouping/cluster these names as those which are similar letter to letter. Are there any

Andrej Borsenkow wrote: > Current aside, the February2000 was last modified on February, 2nd. Last > message there is http://us1.samba.org/listproc/samba/February2000/0042.html > The same is true for samba-technical, samba-cvs (did not check others) > Or is it a problem of my proxy? I don't think so, I see the same thing. Samba web-folks, could you check this out, please?

Hola a todos. En linux uso rkward o gedit con el plugin Rgedit, pero me gustaría saber qué editores se pueden usar en windows. He probado Tinn-R , ¿qué usáis vosotros? Gracias

Hola a todos: Estoy intentado realizar un modelo senoidal de unos datos de radiación solar con el fin de afrontar el relleno de la serie y aplicar la prueba de Thom para verificar su homogeneidad [0]. De momento me encuentro con los siguientes problemas: 1- ¿Existe la prueba de Thom en R? ¿O debo crearme mi propia función? 2- Para la realización del modelo senoidal estoy siguiendo los pasos