search for: pisicandru

Hi everybody, I have a vector of characters and i would like to extract certain parts. My vector is named metr_list: [1] "F:/Naval_Live_Oaks/2005/data//BE.tif" [2] "F:/Naval_Live_Oaks/2005/data//CH.tif" [3] "F:/Naval_Live_Oaks/2005/data//CRR.tif" [4] "F:/Naval_Live_Oaks/2005/data//HOME.tif" And i would like to extract BE, CH, CRR, and HOME in a

Hi everybody, I have a vector of characters and i would like to extract certain parts. My vector is named metr_list: [1] "F:/Naval_Live_Oaks/2005/data//BE.tif" [2] "F:/Naval_Live_Oaks/2005/data//CH.tif" [3] "F:/Naval_Live_Oaks/2005/data//CRR.tif" [4] "F:/Naval_Live_Oaks/2005/data//HOME.tif" And i would like to extract BE, CH, CRR, and HOME in a

Hi, I am using an if else statement inside a function ?. If I use that function I have no problems ?. If I use the function with the if else statement inside a second function I get the following waring: Warning message: In if (pval == 0) p_value <- "< 2.2e-16" else p_value <- pval : the condition has length> 1 and only the first element will be used Using the second

Hi everyone, ? I am trying to get a point of intersection between a polyline and a straight line ?.. and get the x and y coordinates of this point. For exemplification consider this: ? ? set.seed(123) ? k1 <-rnorm(100, mean=1.77, sd=3.33) ?k1 <- sort(k1) q1 <- rnorm(100, mean=2.37, sd=0.74) q1 <- sort(q1, decreasing = TRUE) plot(k1, q1, xlim <- c((min(k1)-5),

Hi everybody, I know this might be very off topic and it took me quite a while to up my courage to post this?. But I remember a thread some time ago about how we can find the packages we need to do specific tasks in R if we don?t know before hand which ones actually do it. Now all the packages are listed alphabetically on the web site. Since I am not very advanced in writing my own functions I

Hi, I would like to find out the coordinates of the intersection points of 2 density curves. I did a search but i didn't get any significant results. I really hope some of you have some ideas. here it is an example: set.seed(123) x1 <- rnorm(100, 1, 1) x2 <- rnorm(100, 0, 1) d1 <- density(x1) d2 <- density(x2) plot(d1, col = 2) lines(d2, col = 4) Now i would like to know

Dear all, I'm using a binomial distribution with a logit link function to fit a GAM model. I have 2 questions about it. First i am not sure if i've chosen the most adequate distribution. I don't have presence/absence data (0/1) but I do have a rate which values vary between 0 and 1. This means the response variable is continuous even if within a limited interval. Should i use

Hi everyone, I have a vector of strings, each string made up by different number of words. I want to get a new vector which has only the first word of each string in the first vector. I came up with this: str <- c('aaa bbb', 'cc', 'd eee aa', 'mmm o n') str1 <- rep(1, length(str)) for (i in 1:length(str)) { str1[i] <- strsplit(str, "

Hi, I have a set of irregular time series and i want to produce a simple plot, with dates on x axis and attribute value of y axis. This is simple enough but my x axis is divided automatically by ticks every 5 years. I would like to have a tick every year at January 1st. I am not sure how i can do that - i end up with something very close to what i want, but it is clunky and not very correct. I

Hi, I am doing a very informal presentation for my office about R capabilities to deal with and analyze spatial data, display data and maps, and connections with GIS. I've used in my presentation info from the CRAN, the spatial Task view, and the more striking graphics examples from http://addictedtor.free.fr/graphiques/thumbs.php and NCEAS

Hi, I am trying to get percentages in a more elegant way. I have a data.frame with locations and values (counts) of species at that location. Each location is repeated for each species i have values for and i would like to get percentages of each species at that location. I am not sure if i am clear in my explanations so i will paste my code below: ##################### > x locat val 1

Hi, I have a matrix of data which i can vizualize as an image - for example. I would like to save this image as a geotiff file or at a tiff file with a world file which holds the projection of my data (ultimately the data represent a map of some sort). I know i can save the data as an ESRI grid, but i am not interested in that. I wonder if anybody knows about any code which will help me do

Hi, I wrote a panel function called panelwhite.corplot. If i use this function in a coplot is working fine, but if i use same function with xyplot i get the error: Error using packet 1: plot.new has not been called yet ..... and so on filling all my plots. So with: coplot(be~ch|dbh1, envmetr1, panel=panelwhite.corplot) I get the expected result with: xyplot(be~ch | dbh1,

Hi, I have opened 2 different graphic devices on screen with windows(3) and Windows(4) (I already have had an other graphic on screen opened in device 2). Now device 3 is inactive and device 4 is active. How i can make device 3 active to plot in device 3 and after make device 4 active to have a different plot in it?? For what i am doing it is important to first open the graphic devices -

Hi again, It seems that dev.set() will solve my problem. It took me quite a bit of time, search and frustration until i stumbled on dev.set ;-) Thanks anyway, Monica From: pisicandru@hotmail.comTo: r-help@r-project.orgSubject: [R] Question about making active a previously opened graphic deviceDate: Fri, 18 Jan 2008 15:02:40 +0000 Hi, I have opened 2 different graphic devices on screen with windows(3) and Windows(4) (I already have had an other graphic on screen opened in devi...

Hi list, It seems that as.numeric is rounding up my numbers that are in a character format. Example: a [1] "776554.45" "776985.31" "776076.03" "776092.01" "776151.42" "776276.97" b <- as.numeric(a) b [1] 776554.4 776985.3 776076.0 776092.0 776151.4 776277.0 I've tried as.numeric(a,2) and as.numeric(a, digits = 2) but it

...ith no problems ... the single thing i've done since Friday was to put "defrag" on drive C, i even didn't re-start the computer ...actually it was re-started several time with same weird problem ... so i suppose something was segmented badly on the hard drive. Monica > From: pisicandru at hotmail.com > To: r-help at r-project.org > Subject: [R] memory problem - failing to load rgl in R 2.7.1 patched > Date: Mon, 4 Aug 2008 12:27:31 +0000 > > >> Hi, >> >> yesterday i had the surprise not to be able to load the package "ca" on R 2.7.0 sayi...

Hi, I am sure there is a reason but ...... why larger decimal numbers get rounded to the nearest integer? Example: a <- 3308000.5 a [1] 3308001 I would like my numbers to be decimals .... since they do represent coordinates and i don't want them rounded .... how can i keep them as they are? Thanks, Monica _________________________________________________________________

...so you can actually learn the language also. I have written about the ease of learning R , if you begin with these two packages first and are working in a commercial data environment on http://decisionstats.com/2008/learning-r-easily-two-guis/ On Thu, May 22, 2008 at 9:30 PM, Monica Pisica <pisicandru@hotmail.com> wrote: > > Hi, > > I am doing a very informal presentation for my office about R capabilities > to deal with and analyze spatial data, display data and maps, and > connections with GIS. I've used in my presentation info from the CRAN, the > spatial Task vie...

Hello list, I am wondering if a joining "one-to-many" can be done a little bit easier. I tried merge function but I was not able to do it, so I end up using for and if. Suppose you have a table with locations, each location repeated several times, and some attributes at that location. The second table has the same locations, but only once with a different set of attributes. I would