similar to: Handle missing values

Hi everyone As I couldn't succeed with manual installation of Rmpi I decided to start again from the beginning. I removed R and MPICH in my Ubuntu Hardy installation. Then, to avoid any dependencies problems I have installed MPICH and R from synaptic, not from sources. But now I can't install Rmpi. An error message appears when trying to install Rmpi, you can find in

Hello everyone, When looking for information about clustering of spatial data in R I was directed towards DBSCAN. I've read some docs about it and theb new questions have arisen. DBSCAN requires some parameters, one of them is "distance". As my data are three dimensional, longitude, latitude and temperature, which "distance" should I use? which dimension is related to

Hi R-users I'm a beginner with R and statistics, so I need some help to start my data analysis. I've been reading some docs and tutorials on R and cluster analysis. I've got a large dataset (102000 points) with values of longitude, latitude and temperature and want to see if I can find groups (clusters). Following some tutorials I can look for principal components but get an error

Hi everyone I am trying to do some time series analysis with daily temperature data (40 years). I have created a zoo object and ts object but can't apply stl function. It says the series is not periodic or has less than two periods. I've searched through google and found a lot of messages about this problem but not a solution/example to look for trend and seasonal component of a

Hi everyone I want to install Rmpi to use R in parallel mode in a Linux cluster (Ubuntu, Hardy Heron). It seems to be properly installed but a problem appears when loading Rmpi library. R version 2.11.1 (2010-05-31) > library("Rmpi") Error: package 'Rmpi' was built before R 2.10.0: please re-install it Should I remove R-2.11 and install R-2.10? I have tried to

Hola a todos Estoy intentando extraer un subconjunto de datos de una serie temporal, concretamente quiero extraer los meses de verano de varios años pero no encuentro la manera. Los datos tienen este formato: FECHA;H_SOLAR;DIR_M;DIR_S;VEL_M;VEL_S;VEL_X;U;V;TEMP_M;HR;BAT;PRECIP;RAD;UVA;UVB;FOG;GRID;

Hola a todos Estoy tratando de analizaruna serie temporal con datos diarios de temperaturas (40 años). He creado un objeto zoo (con ayuda de la lista, gracias) sobre el que encuentro la regresión lineal. He probado también a crear un objeto ts a partir del zoo. El problema que encuentro es que nose puede aplicar la función stl para hallar la componente estacional y la tendencia. Rdice que la

Hola a todos Estoy intentando hacer un análisis rápido de una serie temporal de datos diarios pero me encuentro con algunos problemas. Me gustaría en primera instancia hacer una regresión lineal pero no encuentro la forma. Tras leer los datos diarios creo un objeto de la clase zoo y sobre éste no puedo utilizar lm(). He leído algo sobre dynlm pero no encuentro la forma. Se agradece

R-help - I have this set of aggregated tables (sample data below via dput()). And I would like to have delayValue as the column variables with the "temp" (temp1, temp2, temp3) values as the row variables. However I would like to have the temp variables *aggregated into single rows* so that I have the frequency ("Freq" | counts) of each time each "delayValue" occurs

Hi everybody: Could anybody help how I solve the next problem?. I'm doing interpolation maps of tropospheric ozone of my region, and after create it using IDW, and kriging methods, I want from shapefiles (*.shx, *.shp, *.dbf, *.sbx & *.sbn ) create contour over the interpolation maps. Could anybody tell me how do it?. Thanks in advance. ###################################### Antonio

Dear R-collegues, I'm getting an error message (Error in round) when summarising a glm.nb model, and when using anova.negbin (in R 1.3.1 for windows): > m.nb <- glm.nb(tax ~ areal) > m.bn Call: glm.nb(formula = tax ~ areal, init.theta = 5.08829537115498, link = log) Coefficients: (Intercept) areal 3.03146 0.03182 Degrees of Freedom: 283 Total (i.e. Null); 282

Dear R-users, I'm sure it must be a specific function or a better way to convert matrix to x,y,z coordinates (and viceversa), than my function below (it works). Any help? m2coord <- function(m) { k <- nrow(m)*ncol(m) aa <- data.frame(r=1:k, c=1:k, v=1:k) k <- 0 for (i in 1:nrow(m)) for (j in 1:ncol(m)) { k <- k+1 aa$f[k]=i; aa$c[k]=j; aa$v[k]=m[i,j] } aa } Juli

I don't quite understand the difference between the two methods for performing a chi-squared test on contingency tables: summary(table()) and chisq.test() They may different results. E.g.: aa <- gl(2, 10) bb <- as.factor(c(1,2,2,2,1,2,1,2,2,2,1,2,2,2,1,1,1,2,1,1)) aa <- c(aa, aa) bb <- c(bb, bb) table(aa, bb) summary(table(aa, bb)) chisq.test(aa, bb) Could somebody give me

Hi, Thanks for showing me how to use RODBC and odbcConnect. Now works nicely. The question I've got now is: noms <- list.files(pattern=".DBF") # removing extension names: noms <- sapply(noms, function(x) as.character(strsplit(x,".DBF")) , USE.NAMES=F) for (i in 1:length(noms)) { s <- sqlFetch(bdades, noms[i]) # etc. } But it seems that sqlFetch()

Dear colleges, I'm trying to combine a barplot and a plot in a single figure as follows: data <- 1:6 t <- barplot(data, axes=F) par(new= T) plot(t, data, type="b") However, as you can see in the example, the dots of the second plot do not fall in the midpoint of the bars in the first. Any trick for setting the 2 plots at the same scale? I have unsuccessfully tried: plot(t,

R-users, Can an offset term be included in a Poisson model? I get an error message when trying that: >r3o <- glm(tax ~ areal + offset(o), family=poisson) Error in (if (is.empty.model(mt)) glm.fit.null else glm.fit)(x = X, y = Y, : inner loop 1; can't correct step size In addition: Warning message: Step size truncated due to divergence in: (if (is.empty.model(mt))

Hi everybody: I just generate interpolation maps with differents methods, like IDW and kriging, and now i want to compare the predict values versus real, and i don't who is the commant to do it. I want for example, if I pass from console the coordinates of a place, return me the predict value. Could anybody tell me how should do?. Thanks in advance --

Hi everybody: I just generate different maps of interpolation with different methods but with the same package -(gstat, i used the krige command). And now for evaluate the best method i want to compare the predict values versus experimental, but i don't know how do it. Could anybody tell me how should do it?. Thanks in advance. -- ######################################## Antoni Vi??dez

Dear R collegues, I'm trying to understand what's AIC in R (ver. 1.3.1), and I'm getting a different answer if I look at the AIC(of the fitted model) or the aic in the summary( of the fitted model). Is this correct? Can somebody explain me the difference between the two values? or Is the AIC criterion not appropiated for Poisson models? R session: > t1 <- glm(tax ~ areal,

Dear R collegues, I'm trying to understand what's AIC in R (ver. 1.3.1), and I'm getting a different answer if I look at the AIC(of the fitted model) or the aic in the summary( of the fitted model). Is this correct? Can somebody explain me the difference between the two values? or Is the AIC criterion not appropiated for Poisson models? R session: > t1 <- glm(tax ~ areal,