similar to: Problem installing Rmpi

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

Hi everyone I am new to R and have a question about missing values. I am trying to do a cluster analysis of monthly temperatures and my data are 14 columns with spatial coordinates (lat,lon) and 12 monthly values: /lat - lon - temp1 - //temp2 - temp3 - .... - //temp12/ If I omit missing values (my missing values are 99.00) with /mydata <- na.omit(mydata)/ every row with a

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

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

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

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

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

Hello R users It's some time I am playing with a dataset to do some cluster analysis. The data set consists of 14 columns being geographical coordinates and monthly temperatures in annual files latitutde - longitude - temperature 1 -..... - temperature 12 I have some missing values in some cases, maybe there are 8 monthly valid values at some points with four non valid. I don't want to

Hello everyone I am trying to use CLARA method for finding clusters in my spatial surface temperature data and noticed one problem. My data are in the form lat,lon,temperature. I extract lat,lon and cluster number for each point in the dataset. When I plotted a map of cluster numbers I found empty areas in the map. The point is that the number of points that were assigned a cluster number are

HI R users With clara function I get a data frame (maybe this is not the exact word, I'm new to R) with the following variables: > names(myclara) [1] "sample" "medoids" "i.med" "clustering" "objective" [6] "clusinfo" "diss" "call" "silinfo" "data" I want to

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

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,

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

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

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

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,