Re: [R] plotting rasters - no plot, no error

ubject: Re: [R] plotting rasters - no plot, no error > To: pisican...@hotmail.com; r-help@r-project.org > > On Wed, 15 Apr 2015 at 01:16 Michael Sumner > mailto:mdsum...@gmail.com>> wrote: > On Wed, 15 Apr 2015 at 00:50 Monica Pisica > mailto:pisican...@hotmail.com>> w

[R] plotting rasters - no plot, no error

Hi, I have the current version of R installed on 2 different Windows computers. On one i can plot using the plot function a raster in geotif format, on the other it plots only the axis, an empty color bar and no raster what so ever without generating any errors or warnings. So i suspect som

Re: [R] extracting slot "coords" from a Polygons class object

1]]@Polygons[[1]]@coords > > or > > coordinates(pt1.cpoly@polygons[[1]]@Polygons[[1]]) > >> class(pt1.cpoly) > [1] "SpatialPolygons" > attr(,"package") > [1] "sp" > > So see > ?"SpatialPolygons-class" > for details.

[R] extracting slot "coords" from a Polygons class object

Hi, I am struggling to extract the polygon vertices from a list of an object class "Polygons", specifically the slot "coords". I have a point, and i "draw" a buffer around with gBuffer, i am "extracting" the polygon form the SpatialPolygons class and i end up with a list of 1 one objec

Re: [R] colors in violin plot, ggplot2

;notchwidth = .5, alpha = 0.5, colour = "grey50") + > scale_colour_manual(values = c("orange", "brown", "forestgreen")) + > scale_fill_manual(values = c("orange", "brown", "forestgreen")) > >

[R] colors in violin plot, ggplot2

Hi, I am using ggplot and geom_violin to build a violin plot of some with only 2 categories. All is good except that I cannot set up the colors I want or the violin plots. Either I have same color for both my categories or colors from probably rainbow(2), which are red and blue. What if I wan

[R] problem with user defined panel function in xyplot

Hi everyone,   I am trying to do a horizonplot using my own time series data. I know that there is a horizonplot function in latticeExtra, but on closer examination i think that the graph itself is slightly wrong (it displays some regions as triangles and i think they should be trapezoids, an

[R] Plotted circle does not go through desired points - very long email with code

Hi,   I am trying to solve a problem related to Poincare circles ( for more info http://www.ms.uky.edu/~droyster/courses/spring08/math6118/Classnotes/Chapter09.pdf). In a nutshell i am trying to replicate the method in the above pdf section 9.2.1. that explains in broad terms how to draw the

Re: [R] x, y for point of intersection

on suggested by Don is faster, I haven't checked that. > > Regards, Hans Werner > > > 2011/11/29 Monica Pisica > > > > Hi again, > > > > Working with my real data and not with the little example i sent to the > > list i discovered that segm_distanc

Re: [R] x, y for point of intersection

landt " > wrote: > > >If it's a one off, the identify() function might be of help -- if you > >need something algorithmic it's harder due to floating point stuff and > >sampling frequencies. Let me know if that's the case. > > > >Michael >

[R] x, y for point of intersection

Message: 99 Date: Wed, 23 Nov 2011 08:11:22 + From: Hans W Borchers To: Subject: Re: [R] x, y for point of intersection Message-ID: Content-Type: text/plain; charset="utf-8" Monica Pisica hotmail.com> writes: > Hi everyone, > > I am trying to get a

Re: [R] x, y for point of intersection

4 -0500 > To: pisican...@hotmail.com > > If it's a one off, the identify() function might be of help -- if you need > something algorithmic it's harder due to floating point stuff and sampling > frequencies. Let me know if that's the case. > > Michael >

[R] x, y for point of intersection

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

Re: [R] How to write a shapefile with projection

Hi Patrick, Thanks for letting me know. I mostly use rgdal to read and write rasters so until now i kind of ignore other functionality. Unfortunately i supposed that a package dedicated to shapefiles would be the answer and had the functionality i needed. But rgdal does a nice job in saving my

[R] How to write a shapefile with projection - problem solved

Hi, Sorry i have put such a detailed question to the list about writing a shapefile with projection. I realized that if i use writeOGR from rgdal and not the other write shapefile functions i can get a shapefile with projection recognized by ArcGIS. The command is (in case anybody wonders):  

[R] How to write a shapefile with projection

Hi,   I am trying to write a shapefile with projection. I have my data in a data.frame called try and consists in xy coordinates and a numerical attribute value z1.   Libraries loaded are: sp, rgdal, raster, maptools   head(try) x       y                    z1 1 610237.1    

Re: [R] Getting the values out of histogram (lattice)

Hi everybody, Thank you so much for all the explanations. Now it is much more clear to me. And sorry for the delay in answer but i move office and i was not on line for some time (not that you are really interested in that, though ;-)). I should have sent a dummy example, my real data is a lit

[R] Getting the values out of histogram (lattice)

Hi,   I have a relatively big dataset and I want to construct some histograms using the histogram function in lattice. One thing I am interested in is to look at differences between density and percent. I know I can use the hist function but it seems that this function gives sometimes some wr

Re: [R] matrix manipulations

u use table(): > > > table(d1$group,d1$cat) > cat1 cat2 cat3 cat4 cat5 cat6 > land 0 17 0 0 19 23 > water 0 15 0 0 10 16 > > table(d2$group,d2$cat) > cat1 cat2 cat3 cat4 cat5 cat6 > land 14 0 18 22 0 0 > water 19 0 16 11 0 0 > > Hope this helps. > - Phil Spector > Statisti

[R] matrix manipulations

Hi, I am having some difficulties with matrix operations. It is a little hard to explain it so please bear with me. I have a very large data set, large enough that it needs to be split in parts in order to deal with. I can work things on these "parts" but the problem lies in adding together th

Re: [R] Custom ticks on x axis when dates are involved - many thanks!

: pisican...@hotmail.com > CC: r-help@r-project.org > Subject: Re: [R] Custom ticks on x axis when dates are involved > > On 11/25/2010 06:27 AM, Monica Pisica wrote: > > > > ... > > Now the graph looks very close to what i want, but i know that my ticks > > actually

[R] Custom ticks on x axis when dates are involved

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 ca

Re: [R] Error: could not find function "extract" in package raster

ges: > [1] stats graphics grDevices utils datasets methods base > > other attached packages: > [1] raster_1.6-10 sp_0.9-72 > > loaded via a namespace (and not attached): > [1] grid_2.12.0 lattice_0.19-13 tools_2.12.0 > > On Mon, Nov 8, 2010 at 8:08 AM, Mo

[R] Error: could not find function "extract" in package raster

Hi, I would like to use the function "extract" from package raster and i get the following error: m01e <- extract(marsh01, p) Error: could not find function "extract" marsh01 is a raster object and p is an intersectExtent object that is not null. The package is installed and loaded, i use a

Re: [R] intersection of 2 density curves

2C 2009=2C at 10:01 AM=2C Rainer M Krug wrote: > >> On Thu=2C Oct 8=2C 2009 at 3:48 PM=2C Monica Pisica >> wrote: >> >>> >>> >>> Hi=2C >>> >>> I would like to find out the coordinates of the intersection points >>> of 2 >

Re: [R] intersection of 2 density curves

again=2C =20 Monica > Date: Thu=2C 8 Oct 2009 16:01:36 +0200 > Subject: Re: [R] intersection of 2 density curves > From: r.m.k...@gmail.com > To: pisican...@hotmail.com > CC: r-h...@stat.math.ethz.ch > > > > On Thu=2C Oct 8=2C 2009 at 3:48 PM=2C Monica P

[R] intersection of 2 density curves

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

Re: [R] adding color bar to a graph

> CC: r-h...@stat.math.ethz.ch > Subject: Re: [R] adding color bar to a graph > > Monica Pisica wrote: >> Hi everybody, >> >> I am wondering how i can add a "stretch" color bar / legend to a graph that >> uses colorBrewer to define the colors in it. I will

Re: [R] split strings

wn so far: > > sub(".tif", "", basename(metr_list), fixed = TRUE) > > It does not allow file names with .tif in the middle > of them since it will delete the first occurrence rather > than the last but such a situation is highly unlikely. > > > On Tue, May

Re: [R] split strings

; From: waclaw.marcin.kusnierc...@idi.ntnu.no > To: pisican...@hotmail.com > CC: r-help@r-project.org > Subject: Re: [R] split strings > > Monica Pisica wrote: >> Hi everybody, >> >> I have a vector of characters and i would like to extract certain parts. My >> vector is named

[R] split strings

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.t

Re: [R] statistical significance of accuracy increase in classification

Thanks. I will surely try this as well. Monica > Date: Thu, 26 Feb 2009 08:14:31 -0500 > Subject: Re: [R] statistical significance of accuracy increase in > classification > From: mxk...@gmail.com > To: pisican...@hotmail.com > CC: r-help@r-project.org > > > Do you know about any good refe

Re: [R] statistical significance of accuracy increase in classification

Hi Stefan, Thanks so much. This will help me i am sure. These past 2 days i was away on a trip so please excuse my delayed answer. Monica > From: stefan.ev...@uos.de > To: pisican...@hotmail.com > Subject: Re: [R] statistical significance of accuracy increase in > classification > Date

Re: [R] statistical significance of accuracy increase in classification

Max, thanks for the reply. Yes, the models are done outside R (i will see what i can do to run some of them inside R in the future ) and the sampling is extremely skewed. But we use as truth or reference results from a field exercise where people actually went and gave detailed descriptio

Re: [R] statistical significance of accuracy increase in classification

Hi again, Looking more into test statistics i realized that maybe i can use the power.prop.test to see if the difference between the 2 accuracies are zero or not. Do you have any comments about that? Also, should i considered kappa statistics also a kind of proportion and use the same test? I

[R] statistical significance of accuracy increase in classification

Hi everyone, I would like to test for the statistical significance(for what it worth ...) in increasing classification accuracy and kappa statistics from different land classifications. The classifications were done using other software (like eCognition and See5), but the results were "sampled

Re: [R] joining "one-to-many"

ry merge(t1, t2) > > > On Tue, Feb 17, 2009 at 9:33 AM, Monica Pisica wrote: >> >> 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 usin

[R] joining "one-to-many"

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 se

Re: [R] Aggregrate function

ompact version of the same idea: > > do.call(rbind,by(xveg,xveg['loc'],function(x)x[x$tot == max(x$tot),])) > > - Phil Spector > Statistical Computing Facility > Department of Statistics > UC Berkeley > spec...@stat.berkeley.edu > > > On Thu, 12 Feb 2

Re: [R] Aggregrate function

max(x$tot), ])) > > loc sp tot > L1 L1 b 60 > L2 L2 e 30 > L3 L3 b 68 > > -Christos > >> -Original Message- >> From: r-help-boun...@r-project.org >> [mailto:r-help-boun...@r-project.org] On Behalf Of Monica Pisica >> Sent: Thursday, Februar

[R] Aggregrate function

Hi, I have to recognize that i don't fully understand the aggregate function, but i think it should help me with what i want to do. xveg is a data.frame with location, species, and total for the species. Each location is repeated, once for every species present at that location. For each lo

Re: [R] maptools: Test if point is in polygon

Hi, You may want to check the package sp. There is one function called "point.in.polygo" that might do what you want. Good luck, Monica - >Message: 32 >Date: Thu, 5 Feb 2009 14:51:36 + >From: baptiste auguie >Subject: Re: [R] maptools: Test if point is in polygo

[R] quantiles for sorted pairs of data

Hi everybody, Suppose I have continuous measurements of an energy waveform that is sampled discretely for different heights every 0.5m. Let's say I want to find out the height for which I have equal amount of energy above and below. My colleague did the following: a. calculate the cumulative

Re: [R] Spatial join - optimizing code

ch.and.dist <- function(point) { > sqdists <- colSums((point - ypos)^2) > ind <- which.min(sqdists) > c(ind, sqrt(sqdists[ind])) > } > match <- sapply(xposl, get.match.and.dist) > cbind(xpos, mindist=match[2,], yinfo[match[1,],]) > } > > It's marginally fa

[R] Spatial join – optimizing code

Hi, Few days ago I have asked about spatial join on the minimum distance between 2 sets of points with coordinates and attributes in 2 different data frames. Simon Knapp sent code to do it when calculating distance on a sphere using lat, long coordinates and I've change his code to use Euclidi

[R] Join data by minimum distance

I am wondering if there is a function which will do a join between 2 data.frames by minimum distance, as it is done in ArcGIS for example. For people who are not familiar with ArcGIS here it is an explanation: Suppose you have a data.frame with x, y, coordinates called track, and a second data

[R] Adding a legend to R graph device with several plots (no to individual plots!)

Hi Nelson, I don't know if you got your answer or not but here it is my low-tech solution for what is worth. I like to do a layout as I want my graphs to be and add some space where I "plot" an invisible graph with dummy data and add the legend there. In this case I decided to put the legend o

[R] Non-constant variance and non-Gaussian errors

Hi Paul, Take a look at gam() from package mgcv (gam = generalized additive models), maybe this will help you. GAMs can work with other distributions as well. Generalized additive models consist of a random component, an additive component, and a link function relating these two components. Th

Re: [R] savePlot() does not save plot with format set

Did you try: savePlot("test.bmp", type = "bmp") Monica Message: 118 Date: Tue, 26 Aug 2008 10:29:13 +0100 From: "Luis Ridao Cruz" Subject: [R] savePlot() does not save plot with format set To: Message-ID: Content-Type: t

Re: [R] GAM-binomial logit link

de the formula for discovering possible problems.Fabrizio Cipollini ---- Message: 54Date: Thu, 21 Aug 2008 20:09:08 +From: Monica Pisica <[EMAIL PROTECTED]>Subject: Re: [R] GAM-binomial logit linkTo: Message-ID: <[EMAIL PROTECTED]>

Re: [R] GAM-binomial logit link

Hi, I am not sure it is the best to use a binomial distribution for a continuous bounded variable. A beta distribution would be more appropriate, although I don't know how to define one for the gam() function. On the other hand beta distribution is closely linked to the gamma distribution so m

Re: [R] : bquote inside legend()

nt > > val <- substitute("value" == a, list(a = a)) > legend(6, 5, do.call("expression", list("Legend", val)), bty = 'n') > > ? > > Your first attempt is a character vector and your second fixed up quote > and substitute. But more

Re: [R] : bquote inside legend() - TANKS!

t; "), title="Legend" ,bty = 'n') > > val <- substitute("value" ==.(a), list(a = a)) > legend(6, 5, do.call("expression",list("Legend", val)), bty = 'n') > > # --- > # Code end > > On Mon, Aug 11,

[R] : bquote inside legend()

Hi, I have a graph and I would like to write some values inside the legend that were saved in a variable. Please revise the code below in which I've wrote 2 different legends, but I am not happy with either of them. What I want is a legend with tile "Legend" and underneath a line with a name a

[R] SOLVED - memory problem - failing to load rgl in R 2.7.1 patched

Hi, Sorry for that . Today R did load both ca and rgl packages with 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

[R] memory problem - failing to load rgl in R 2.7.1 patched

> Hi, > > yesterday i had the surprise not to be able to load the package "ca" on R > 2.7.0 saying that cannot find required package rgl although it was there. So > today i've upgraded to 7.2.1. patched and i got the following error: > >> local({pkg <- select.list(sort(.packages(all.available =

[R] failing to load package in R 2.7.1 patched

Hi, yesterday i had the surprise not to be able to load the package "ca" on R 2.7.0 saying that cannot find required package rgl although it was there. So today i've upgraded to 7.2.1. patched and i got the following error: > local({pkg <- select.list(sort(.packages(all.available = TRUE))) +

Re: [R] Warning message in if else statement

CTED]> Subject: Re: > [R] Warning message in if else statement> To: [EMAIL PROTECTED]> > hi monica: > i don't understand the statement below. could that be related > to your > problem ?> > if (pval==0) p_value <- "<- pval> > > > > > On

Re: [R] Warning message in if else statement

. If I use this function as it is …. No problems, if I incorporate inside a different function so to be used internally I get the warning. Any idea why? Thanks, Monica > Date: Wed, 23 Jul 2008 12:24:55 -0500 > From: [EMAIL PROTECTED] > To: [EMAIL PROTECTED] > CC: r-help@r-project.org &

[R] Warning message in if else statement

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 con

Re: [R] saving plot both as jpg and pdf

1) 408-8111> > > > > -----Original Message-> > From: [EMAIL PROTECTED]> > [mailto:[EMAIL PROTECTED] On Behalf Of Monica Pisica> > Sent: Tuesday, July 22, 2008 7:50 AM> > To: r-help@r-project.org> > Subject: [R] saving plot both as jpg and pdf> >

[R] saving plot both as jpg and pdf

Hi, I want to save a plot automatically as a pdf and jpg, and if I open the pdf device first and jpeg second only the jpeg file saves correctly …. If I do reverse, only the pdf file saves correctly. pdf('E:/my_graphs/test1.pdf', bg = "white") jpeg('E:/my_graphs/test1.jpg', quality = 100, bg =

[R] Comparing differences in AUC from 2 different models

Hi, I would like to compare differences in AUC from 2 different models, glm and gam for predicting presence / absence. I know that in theory the model with a higher AUC is better, but what I am interested in is if statistically the increase in AUC from the glm model to the gam model is signifi

Re: [R] PCA on image data

Hi, You have the values of the principal component in "scores" (for each "location" where you have a RGB reading) and the eigenvectors in "loadings", see ?princomp So your first pca component would be: comp1 <- pca$scores[,1] Now you can transform this in the matrix you need to display your

Re: [R] tiff()-bug (was re:Preparing high quality figures

Hi, I don't know if this matters but it worked for me with no problems …. > sessionInfo() R version 2.7.0 (2008-04-22) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_Unit

Re: [R] loop with files

Hi, I am not sure this is what you want …. But alist <- list.files(path = "whatever path you have for your files", full.names = TRUE) Now you have a list with the names of your files …. Form each name you can make a variable name and assign that file …. fname <- substring(alist[i], first, l

[R] PCA analysis

Monna, The way i do it is to re-create the biplot for the PCA . I am attaching my code (i am sure this can be done even easier . but this works as well) where i am using pca() function from labdsv and my data is called veg1. library (labdsv) pca.1<-pca(veg1,cor=TRUE) # The scores are

Re: [R] Pros and Cons of R

em ;-) Thanks so much for your replies, Monica> Date: Thu, 22 May 2008 12:14:11 -0400> From: [EMAIL PROTECTED]> To: [EMAIL PROTECTED]> CC: r-help@r-project.org> Subject: Re: [R] Pros and Cons of R> > Monica Pisica wrote:> > > > Cons:> > > > -

Re: [R] Pros and Cons of R

Lenovo ThinkPad and was able to allocate memory beyond 2G. Have not done much else otherwise but it seems to work just fine.> > H> > -Original Message-> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Monica Pisica> Sent: Thursday, May 22, 2008 9:00 AM>

Re: [R] Pros and Cons of R

rate how to do familiar things using R.> > > Charles Annis, P.E.> > [EMAIL PROTECTED]> phone: 561-352-9699> eFax: 614-455-3265> http://www.StatisticalEngineering.com> > > -Original Message-> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On> B

[R] Pros and Cons of R

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://ad

Re: [R] applying cor.test to a (m, n) matrix

AIL PROTECTED]> wrote: have a look also at function rcor.test() from package ltm.Best,DimitrisDimitris RizopoulosBiostatistical CentreSchool of Public HealthCatholic University of LeuvenAddress: Kapucijnenvoer 35, Leuven, BelgiumTel: +32/(0)16/336899Fax: +32/(0)16/337015Web: http://med.kuleuven

Re: [R] applying cor.test to a (m, n) matrix - SUMMARY

Hi again, I've got few very good options from the list and since they were not posted to the list, I will provide a summary. Thank you very much to all who answered and I hope this summary will benefit others interested in solving similar problems like that. Yasir Kaheil re-wrote my original

[R] applying cor.test to a (m, n) matrix

Hi everybody, I would like to apply cor.test to a matrix with m rows and n columns and get the results in a list of matrices , one matrix for p.val, one for the statistic, one for the correlation and 2 for upper and lower confidence intervals, something analog with cor() applied to a matrix.

Re: [R] Correspondence and detrended correspondence analysis

; CC: r-help@r-project.org> Date: Sat, 19 Apr 2008 11:18:25 +0100> > On Fri, 2008-04-18 at 18:45 +, Monica Pisica wrote:> > Hi,> > > > I hope someone knows the answer to this or has a real good reference> > about it (I am using Legendre & Legendre, Numer

[R] Correspondence and detrended correspondence analysis

Hi, I hope someone knows the answer to this or has a real good reference about it (I am using Legendre & Legendre, Numerical Ecology, 1998) My data is a data.frame with locations as rows and vegetation assemblages / species as columns. I've done a PCA, a correspondance analysis (CA) using

Re: [R] unexpected GAM result - at least for me!

ll the help, Monica > Date: Mon, 31 Mar 2008 09:30:01 -0400 > From: [EMAIL PROTECTED] > To: [EMAIL PROTECTED] > CC: r-help@r-project.org > Subject: Re: [R] unexpected GAM result - at least for me! > > On 3/31/2008 9:01 AM, Monica Pisica wrote: >> Thanks Duncan. >

Re: [R] unexpected GAM result - at least for me!

predict presence / absence. Do you think it is still a valid result? Thanks again, Monica> Date: Mon, 31 Mar 2008 08:47:48 -0400> From: [EMAIL PROTECTED]> To: [EMAIL PROTECTED]> CC: r-help@r-project.org> Subject: Re: [R] unexpected GAM result - at least for me!> > On 3/31

[R] unexpected GAM result - at least for me!

Hi I am afraid i am not understanding something very fundamental and does not matter how much i am looking into the book "Generalized Additive Models" of S. Wood i still don't understand my result. I am trying to model presence / absence (presence = 1, absence = 0) of a species using

Re: [R] dreaded p-val for d^2 of a glm / gam

say, 1000 random permutations of your response variable, collect > anova(glm.D93t, glm.D93)[2, "Deviance"] in a vector, and then find out > how extreme the deviance you actually got is relative to this > permutation distribution. > > Hope this helps. > Spencer Graves> p

[R] dreaded p-val for d^2 of a glm / gam

OK, I really dread to ask that much more that I know some discussion about p-values and if they are relevant for regressions were already on the list. I know to get p-val of regression coefficients - this is not a problem. But unfortunately one editor of a journal where i would like to pu

Re: [R] a more elegant way to get percentages?

c", "c", "c", "d", "d", "d"), > val = c(5, 5, 15, 5, 20, 5, 10, 5, 15, 10)) > > Try this: > > x$percent1 <- ave(x$val, x$locat, FUN = function(x) 100*x/sum(x)) > > On Thu, Mar 13, 2008 at 9:36 AM, Monica Pisica wrote: >

[R] a more elegant way to get percentages?

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

[R] How to Read?

Well, i am not sure what do you want how to do the graphs you showed or how to interpret them?? If you want to know how to do them, look at package ade4, especially functions s.corcircle and scatter.dudi, but you can use the function plot as well and choose the right variables since PCA

Re: [R] newbie (me) needs to model distribution as two overlapping gaussians

take a look at: Du, 2002, Master Thesis, http://www.math.mcmaster.ca/peter/mix/Rmix.pdf Macdonald, P., 2003, RMIX routine for R, http://www.math.mcmaster.ca/peter/mix/mix.html I don't think this package was actually posted on CRAN (the mix package on CRAN is a different one as far as i reme

Re: [R] How to search for packages - wrap up!

Hi everybody, I have to recognize that my post certainly shows my lack of skills in really navigating the R web page i am surprised that only one or two persons wrote me about the "Task Views" - it is what i was after - although maybe too general for my lazy taste - but hei - it is ther

Re: [R] How to search for packages

eb 2008 10:03:38 -0600 > From: [EMAIL PROTECTED] > To: [EMAIL PROTECTED] > Subject: Re: [R] How to search for packages > CC: r-help@r-project.org > > On Feb 4, 2008 8:34 AM, Monica Pisica wrote: >> >> Hi everybody, >> >> >> I know this might be ve

Re: [R] How to search for packages

Hi Jonathan, I didn't know about your web page - now it is in my favorites. I am sure it will help a lot. Google does and doesn't do the job - sometimes ;-) But R and packages are an open source and i realize that people already put a lot of good effort into something which is essentially fre

[R] How to search for packages

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 package

Re: [R] as.numeric rounds up - Thanks!

.org> > options(digits=10)> > > b> > or just> > print(b, 10)> > On 01/02/2008, Monica Pisica <[EMAIL PROTECTED]> wrote:> >> >> > Hi list,> >> > It seems that as.numeric is rounding up my numbers that are in a character for

[R] as.numeric rounds up

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

Re: [R] titles with superscript and variale value

titles with superscript and variale value> CC: r-help@r-project.org> > Try this> > plot(rnorm(100), main=bquote(Family ~ Gaussian ~ ";" ~ D^2 == .(d2)))> > > On 27/01/2008, Monica Pisica <[EMAIL PROTECTED]> wrote:> >> >> > Hi everyone,>

[R] titles with superscript and variale value

Hi everyone, I am trying to write a title for a plot which has a superscript like D^2 and a value of a variable stored in that variable. i am not sure if i am clear so i will try an example: Suppose i want my title to be like: Family: Gaussian; D^2 = 0.45 where i have the value 0.45

[R] modifing barplot

Hi, I would like to modify how barplot plots. For reasons only my supervisor comprehends i need to do the crisscross lines in the bar thicker. If i use lwd=2 the y-axis is thicker, (i use the option beside = TRUE) but nothing else is affected by this parameter - it seems. I looked at th

[R] FW: Question about making active a previously opened graphic device - SOLVED

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: [EMAIL PROTECTED]: [EMAIL PROTECTED]: [R] Question about making active a previously opened graphic deviceDate: Fri,

[R] Question about making active a previously opened graphic device

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 hav

[R] error in my selection

Hi everybody, Well, I really hoped that until now i know how to do a selection . but obviously i don't. So risking your laughing at me - here you go. I have a table with lots of vegetation data, and one column specifies in which region the data is. There are only 2 regions, region 1 and

[R] formatting ftable

Hi all, I am using ftable and xtabs to get total counts of species per location per hight category (understory, midstory and overstory). i can save the data in almost the format i would like it to be - but not quite. Here it is an example of what i get: x - is a dataframe with the following col

[R] formating ftable

Hi all, I am using ftable and xtabs to get total counts of species per location per hight category (understory, midstory and overstory). i can save the data in almost the format i would like it to be - but not quite. Here it is an example of what i get: x - is a dataframe with the following

Re: [R] cor.test formula

Hi, Thanks. I usually don't fail to google though .. I really appreciate your answer. It is exactly what i needed. Monica> Date: Mon, 17 Dec 2007 10:50:30 -0500> From: [EMAIL PROTECTED]> Subject: Re: [R] cor.test formula> To: [EMAIL PROTECTED]> CC: [EMAIL PROTECT

  1   2   >