[R] geotiff calculations

2007-11-14 Thread Laura Poggio
Dear list, I have to compare two digital elevation models in raster format (geotiff). I then have to calculate the differences in altitude for each cell and make some statistics (basic as mean, median, std, range but also more advanced as RMSE) on that. I do not know very much how to proceed: 1)

[R] R function, sink() and empty file

2008-08-10 Thread Laura Poggio
Dear all, I wrote a simple script in order to put together some functions and method to be executed on various files I am trying to have to possibility to call the script changing few parameters in order to use the different files. I succeeded partly using the function method. However in my script

[R] cluster.stats

2008-06-13 Thread Laura Poggio
Dear list, I just tried to use the function cluster.stat in the package fpc. I just have a couple of questions about the syntax: cluster.stats(d,clustering,alt.clustering=NULL, silhouette=TRUE,G2=FALSE,G3=FALSE) 1) the distance object (d) is an object obtained by the function dist() on my own

Re: [R] cluster.stats

2008-06-14 Thread Laura Poggio
Thank you very much for your answer. I tried to run the function on my data and now I am getting this message of error Error in as.dist(dmat[clustering == i, clustering == i]) : (subscript) logical subscript too long Below the code I am using (version2.7.0 of R with all packages updated): data

Re: [R] cluster.stats

2008-06-14 Thread Laura Poggio
) whether data is what you expect it to be. Christian On Sat, 14 Jun 2008, Laura Poggio wrote: Thank you very much for your answer. I tried to run the function on my data and now I am getting this message of error Error in as.dist(dmat[clustering == i, clustering == i]) : (subscript

Re: [R] cluster.stats

2008-06-14 Thread Laura Poggio
know ... I simply try to see if something was changing with different structure of data Christian On Sat, 14 Jun 2008, Laura Poggio wrote: I am sorry I did not provide enough information. I am not using img later, but data that is data.frame. I wrote that img is a image just to explain

Re: [R] cluster.stats

2008-06-14 Thread Laura Poggio
. Sorry, Christian On Sat, 14 Jun 2008, Laura Poggio wrote: Thank. See below. Laura 2008/6/14 Christian Hennig [EMAIL PROTECTED]: What does str(ddata) give? Class 'dist' atomic [1:130816] 69.2 117.1 145.6 179.9 195.6 ... dcent doesn't make sense as input

[R] variable as part of file name

2008-07-02 Thread Laura Poggio
Dear all, sorry for this very basic question, but I did not find any good example yet. I would like to set up a variable that can be recall later to substitute a part of a file name. As example: var_filename = as.name(aaa) jpeg(var_filename.jpg) plot() dev.off() It would be very useful in

Re: [R] variable as part of file name

2008-07-02 Thread Laura Poggio
Thank you perfect! it worked also with the GDAL option. Laura 2008/7/2 Richard Pearson [EMAIL PROTECTED]: Laura Does jpeg(paste(var_filename, jpg, sep=.)) do what you want? Regards Richard Laura Poggio wrote: Dear all, sorry for this very basic question, but I did not find any good

[R] x axe values

2008-11-07 Thread Laura Poggio
Dear list, I have to draw a simple plot. On y axe some numerical values that correspond to various categories on axe x. The table I am reading looks like: cat Obj1 Obj2 Obj3 max 23 27 34 ave 21 25 32 min 19 23 30 In order to avoid that the first column is reordered alphabetically I used: (found

Re: [R] R-square in robust regression

2008-11-13 Thread Laura Poggio
to the Proportion of variation in response(s) explained by model(s) computed by S-Plus. There are several options. I have had good results using wle.lm() in package wle and lmRob() in package robust. The second option is perhaps closest to what you want. Regards, Mark. Laura POggio wrote: I

Re: [R] R-square in robust regression

2008-11-13 Thread Laura Poggio
? Thank you very much in advance. Laura Poggio - Date: Mon, 20 Oct 2008 06:15:49 +0100 (BST) From: Prof Brian Ripley [EMAIL PROTECTED] Subject: Re: [R] R-square in robust regression To: PARKERSO [EMAIL PROTECTED] Cc: r

Re: [R] R-square in robust regression

2008-11-13 Thread Laura Poggio
== Laura Poggio [EMAIL PROTECTED] on Thu, 13 Nov 2008 10:43:14 + writes: LP yes thank you! it is perfect. LP I was using lmrob in package robustbase and it did not have that option in LP the summary. Yes lmRob() from robust is from a company which -- often being excellent