Re: [R] Announce: edtdbg, integrating R's debug() with your text editor

Hello On 12/7/09, Norm Matloff wrote: > Please give me your comments. I probably will upload to CRAN at some > point, possibly after an ESS guru contributes the ESS code. :-) > Do you think that the package can be made to work with Geany? Thank you Liviu _

[R] Announce: edtdbg, integrating R's debug() with your text editor

I've just developed edtdbg, a small package that integrates R's debug() with one's text editor. Excerpt from the README file: Goal The debug() function in R is primitive. My goal was to make it more usable by integrating it with one's text editor. Hence I wrote the package here

Re: [R] [ggplot2] Wind rose orientation

Hi Tom, The problem was that the bars have width 22.5, so setting the limits to c(0, 360) chopped off the last bar. To get around that, you needed to add half a bar width (11.25) to the start and end limits, and then rotate the polar coordinate system so that 0 lines up again. ggplot(wind.data,

Re: [R] Label individual points in lattice by and ID

Sorry. My previous code didn't quite display my problem correctly. To make my point more clear, I want to use my "groups" as the label for the point, or rather the point itself. The key becomes illegible with many groups. x<-as.factor(1:5) y<-rnorm(1:10) f1<-c("a","b") f2<-c("x","y") dat<-expa

Re: [R] hi there!

Hi Dave, May be a good idea you (1) use subjects avoiding "help", "Hi there", "problem", etc. The list is archived for future search, and nobody with problems with MCMC, Metropholis Hasting algorithm will search "Hi there"; (2) how about minimum reproducible problem? If you are needing use this ki

Re: [R] 4D plotting

Hi Jared, google R graph library 4d or http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=90 bests milton On Sun, Dec 6, 2009 at 7:31 AM, Jared Nance wrote: > Hello list, > Thanks in advance for whatever help you can give me. I have a data set > that I want to visualize that has fo

[R] Label individual points in lattice by and ID

I am using a plot to inspect data points, and I would like to identify each point with respect to an ID. At issue is that I am producing a faceted plot with many IDs (96) and the key is far to large to accurately identify points by color. 1) Can you direct me on labeling or printing data points by

Re: [R] rgl keyboard shortcut for translation on Mac?

Thanks Duncan... I see how to do it. It appears that all 3 of my mouse buttons are actually set to the trackpad button on my MacBook Pro. So if I assign any button via pan3d I get the same behavior, and can't get back to the default. I'll play with it! Bryan On 12/6/09 9:14 PM, "Duncan Murdoc

Re: [R] rgl keyboard shortcut for translation on Mac?

Bryan Hanson wrote: Hello Everyone. I¹m on a Mac, using rgl. Thanks to all the developers and maintainers on it! If I drag the mouse, I can rotate the view, if I hold down the option key or the ctrl key and drag, I get scaling (though at times I seem to alter the perspective?). Adding the shi

Re: [R] help with repeated values

Here you have an example: > colleges <- factor( c("a", "a", "b", "c", "c" )) > colleges [1] a a b c c Levels: a b c > table( colleges ) colleges a b c 2 1 2 > colleges.count <- table( colleges )[colleges] > colleges.count colleges a a b c c 2 2 1 2 2 > which( colleges.count > 1 ) a a c c 1 2 4 5

Re: [R] bootstrap help

Hi there, This seems to 'simple' for that. ( I mean the codes are too good) Can I ask for a bit explaination? or some simplier (maybe few more line codes) to archieve that goal? And the codes you gave here, first line give a table with x and y , random normal values, but the second line, what d

Re: [R] help

Tena koe Kristen all.equal(theOneVariable, theSum) HTH Peter Alspach PS It is advisable to use a more informative subject line. > -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Kristen Mark > Sent: Monday, 7 December 2009

Re: [R] help with repeated values

college=c("Columbia","Columbia","Harvard","Harvard","NYU") unique(college) table(college) unique(college)[table(college)>1] HTH Daniel - cuncta stricte discussurus - -Ursprüngliche Nachricht- Von: r-help-boun...@r-project.org [mailto:r-h

[R] Bootstrap pvalues

Hi, How do I generate bootstrap p-values for goodness of fit test. Thank you in advance for your help, Juliet [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] estimate inverse gaussian in R

noclue_ wrote: I have a one-variable data set in R. The plot of histogram of my numerical variable suggests an inverse gaussian distribution. How can I obtain best estimation for the two parameters of inverse gaussian based on my data? Thanks. See the examples for function pelp in packag

[R] help with repeated values

I have a huge college dataset and I want to be able to find which colleges appear more than once in the dataset. How would I go about doing this? The variable college is a factor. If I just do table(college) it is impossible to find where there are repetitions. Any help would be greatly apprecia

[R] help

How would I verify that one variable is the summation of two other variables? Thank you! [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] bootstrap help

Hi， I have a dataframe with x and y in pairs, ( 49 pairs i.e. 49 rows ) x y 76 80 138 143 67 67 29 50 381 464 23 48 37 63 120 115 how do i generate 50 Bootstrap samples and corresponding estimates? i have tried sample( data, size=50,replace=T) does not seem to be right Thanks.

[R] Function for MFE and MSFE

Hello, Is there a function that calculates MFE(Mean Forecast Error) and MSFE(Mean Square Forecast Error) in R? Cant find anything! Thanks Tobias -- View this message in context: http://n4.nabble.com/Function-for-MFE-and-MSFE-tp950368p950368.html Sent from the R help mailing list archive at Nab

[R] yeroon.net/ggplot2 web application v0.11

A new version of the ggplot2 web application is available at http://www.yeroon.net/ggplot2. New features include 1D geom’s (histogram, density, freqpoly), syntax mode (by clicking the tiny arrow at the bottom), and some additional facet options. Furthermore some minor improvements and fixes, most n

Re: [R] optim with constraints

Hi, Prof Nash Thanks for your comment! I modified my code to be (added an extra parametr): optim(c(1.14,0.25,0.06), weibull.like, mydata=mydata, method="L-BFGS- B", hessian = TRUE, lower = c(0, 0, 0), upper = c(Inf, Inf, 1)) But I had the following error: Error in optim(c(1.14, 0.25, 0.06), we

[R] rgl keyboard shortcut for translation on Mac?

Hello Everyone. I¹m on a Mac, using rgl. Thanks to all the developers and maintainers on it! If I drag the mouse, I can rotate the view, if I hold down the option key or the ctrl key and drag, I get scaling (though at times I seem to alter the perspective?). Adding the shift key doesn't seem to

Re: [R] estimate inverse gaussian in R

David Scott wrote: Jorge Ivan Velez wrote: Hi Tim, Take a look at http://bm2.genes.nig.ac.jp/RGM2/R_current/library/fBasics/man/dist-nigFit.html HTH, Jorge' i don't think this is what the OP wanted. The nig is the normal inverse Gaussian distribution, which is a special case of the general

[R] igraph plot - vertex colors

Hi everyone, I have successfully created and analyzed my network data. I am new to R, and Network Analysis too, but I want to color my vertex based on some of the centrality measures calculated. Can someone point me in the right direction? Again, I am new to R, but given how powerful R appe

Re: [R] bootstrap help

data<-data.frame(x=rnorm(49), y=rnorm(49)) t(sapply(1:50,function(i){colMeans(data[sample(nrow(data),size=nrow(data),replace=TRUE),])})) r-help@r-project.org 2009/12/7 casperyc : > > Hi there, > > i think that's not what i was aiming for... > > i was aked to > > Generate 50 Bootstrap samples and

[R] confint for glm (general linear model)

Hi, I have a glm gives summary as follows, Estimate Std. Errorz valuePr(>|z|) (Intercept) -2.03693352 1.449574526 -1.405194 0.159963578 A0.01093048 0.006446256 1.695633 0.089955471 N0.41060119 0.224860819 1.826024 0.

Re: [R] bootstrap help

Hi there, i think that's not what i was aiming for... i was aked to Generate 50 Bootstrap samples and corresponding estimates if i do data[sample(nrow(data),size=50,replace=TRUE),] it will give me a table of 50 rows ( 50 sets of x and y) then how do i estimate the mean? ( mean was esitmated b

Re: [R] bootstrap help

casperyc hotmail.co.uk> writes: > I have 49 pairs in my data.frame 'data' > > x y > 7680 > ... ... > > how do I get a bootstrap sample of size n=50? > data[sample(nrow(data),size=50,replace=TRUE),] __ R-help@r-project.org mailing list http

Re: [R] SAS "datalines" or "cards" statement equivalent in R?

Gary, if i were you, i would use scan(). here is a piece of code. # DO DATA INPUT IN R CONSOLE WITH SCAN() # #--# # COMPARABLE SAS CODE: # #--# # data t

Re: [R] estimate inverse gaussian in R

Jorge Ivan Velez wrote: Hi Tim, Take a look at http://bm2.genes.nig.ac.jp/RGM2/R_current/library/fBasics/man/dist-nigFit.html HTH, Jorge' i don't think this is what the OP wanted. The nig is the normal inverse Gaussian distribution, which is a special case of the generalized inverse Gaussi

[R] bootstrap help

Hi, I have 49 pairs in my data.frame 'data' x y 76 80 138 143 67 67 29 50 381 464 23 48 37 63 120 115 ... ... how do I get a bootstrap sample of size n=50? i have tried sample(data,size=50,replace=TRUE) and sample(data,replace=TRUE) both didnt see

Re: [R] Dividing a pixel image into factors - (cut.im(), cut.default())

On 4/12/2009, at 10:10 PM, Gough Lauren wrote: Hi, I have a numeric pixel image which I would like to divide into factors for analysis in Spatstat. I have found that I can use cut.im() function to divide the range of pixel values into a series of equal length intervals (e.g. if my pixels

Re: [R] Use of apply rather than a loop

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of William Dunlap > Sent: Friday, December 04, 2009 3:29 PM > To: Dennis Fisher; r-h...@stat.math.ethz.ch > Subject: Re: [R] Use of apply rather than a loop > > You could try using

Re: [R] Passing arguments to gpar

Hi Sebastien Bihorel wrote: Dear R-users, I would like to know how to pass arguments to gpar() without hard-coding them. I tried to store my arguments in a list and passed this list to gpar(), but it did find the way to do it properly. Any help would be appreciated. a<- list(fontisze=8,co

Re: [R] .libPaths(new) stopped working in 2.10

On 7/12/2009, at 8:32 AM, Alexy Khrabrov wrote: I used to have the following in my .Rprofile: if (length(.libPaths())==1) .libPaths(paste(Sys.getenv("HOME"),"/Library/R/",paste(R.version $major,as.integer(R.version$minor),sep='.'),"/library",sep='')) -- and it added my user-defined library

Re: [R] .libPaths(new) stopped working in 2.10

On Sun, Dec 6, 2009 at 3:04 PM, Duncan Murdoch wrote: > Does that directory exist?  .libPaths("foobar") silently does nothing, > because (on my system) "foobar" is not a directory. (This is documented > behaviour, though it would perhaps be friendlier if it gave a warning when > it dropped a reque

Re: [R] ifelse and wildcard matching

On 06/12/2009 3:05 PM, Jennifer Brea wrote: Hi, I'm trying to create a new variable using the following code: A$B2<-ifelse(A$B %in% "EDUCATION","EDUCATION",NA) The problem is, I want the new variable to be labeled "EDUCATION" if the original *contains* "EDUCATION", not of it is an exact match

[R] ifelse and wildcard matching

Hi, I'm trying to create a new variable using the following code: A$B2<-ifelse(A$B %in% "EDUCATION","EDUCATION",NA) The problem is, I want the new variable to be labeled "EDUCATION" if the original *contains* "EDUCATION", not of it is an exact match. I've done a lot of searching on wild card

Re: [R] .libPaths(new) stopped working in 2.10

On 06/12/2009 2:53 PM, Alexy Khrabrov wrote: I have the .libPaths() defined in my .Rprofile, and it stopped having any effect. Even if I try to do it in the REPL, as: x <- paste(Sys.getenv("HOME"),"/Library/R/",paste(R.version$major,as.integer(R.version$minor),sep='.'),"/library",sep='') x

Re: [R] .libPaths(new) stopped working in 2.10

I have the .libPaths() defined in my .Rprofile, and it stopped having any effect. Even if I try to do it in the REPL, as: > x <- > paste(Sys.getenv("HOME"),"/Library/R/",paste(R.version$major,as.integer(R.version$minor),sep='.'),"/library",sep='') > x [1] "/Users/alexyk/Library/R/2.10/library"

Re: [R] regression tests still fail for R version 2.10.0

On Sat, Dec 05, 2009 at 01:39:03PM -0500, David Winsemius wrote: > The question lacks so many specifics that would be needed to address > problem that I suspect persons with sufficient knowledge to address it > (and I am probably not among that group) simply passed it over and > waited for yo

[R] .libPaths(new) stopped working in 2.10

I used to have the following in my .Rprofile: if (length(.libPaths())==1) .libPaths(paste(Sys.getenv("HOME"),"/Library/R/",paste(R.version$major,as.integer(R.version$minor),sep='.'),"/library",sep='')) -- and it added my user-defined library directory. Then I installed packages there, so during

Re: [R] R + Hull-White model using nonlinear least squares

On Dec 6, 2009, at 1:08 PM, king123 wrote: Hi guys I have data that contains the variances vt of the yields of 1, 2, 3, 4, 5,10, 20 year bonds. Assuming the Hull-White model for the yield of a t-year zero-coupon bond, I have to estimate the σ of the Hull-White model using nonlinear leas

[R] Evaluation error

Hello, I want to write a function calculating posterior probability. Can anyone help have a look at the following program? Here function f is what I eventually need. It is related to function f1, f2 and f3. The argument x is a vector. I get evaluation errors here. Thank you very much! > f1

Re: [R] estimate inverse gaussian in R

Hi Tim, Take a look at http://bm2.genes.nig.ac.jp/RGM2/R_current/library/fBasics/man/dist-nigFit.html HTH, Jorge' On Sun, Dec 6, 2009 at 12:48 PM, noclue_ <> wrote: > > I have a one-variable data set in R. > > The plot of histogram of my numerical variable suggests an inverse > gaussian distrib

Re: [R] estimate inverse gaussian in R

On Dec 6, 2009, at 12:48 PM, noclue_ wrote: I have a one-variable data set in R. The plot of histogram of my numerical variable suggests an inverse gaussian distribution. How can I obtain best estimation for the two parameters of inverse gaussian based on my data? http://www.lmgtfy.com/?q=

Re: [R] R ANOVA question

On Dec 6, 2009, at 1:14 PM, Xin Ge wrote: Hi All, I'm fitting an intercept ANOVA model using R. I'm using following contrast statement, so my first level is considered as base level. options(contrasts = c("contr.treatment", "contr.treatment")) My data has two factors are "Month" and "WeekD

Re: [R] Error in eval(expr, envir, enclos) : object 'N' not found

hansoffate wrote: Oh yes, that's true, but I don't understand why the same exact model works fine the first time I run LSODA, then running the model again with different PARMS causes it to error with this error. I've been trying to debug or find other relevant posts with this same error, but so

[R] R + Hull-White model using nonlinear least squares

Hi guys I have data that contains the variances vt of the yields of 1, 2, 3, 4, 5,10, 20 year bonds. Assuming the Hull-White model for the yield of a t-year zero-coupon bond, I have to estimate the σ of the Hull-White model using nonlinear least squares and give a 95% confidence interval for each

Re: [R] help

On Dec 6, 2009, at 7:19 AM, Kristen Mark wrote: Hello there, How would I go about creating a function that will read a vector containing names of variables and then will output tables of proportions for all the categorical variables and summary statistics appropriate for all quantitative

[R] estimate inverse gaussian in R

I have a one-variable data set in R. The plot of histogram of my numerical variable suggests an inverse gaussian distribution. How can I obtain best estimation for the two parameters of inverse gaussian based on my data? Thanks. -- View this message in context: http://n4.nabble.com/estimate

[R] conout

In order to get user input when running a batch file with Rterm in Windows XP and Windows Server 2003, I have been using the R script below with R-2.8.1 and earlier. This opens a command line window with the statement "Please enter the ID:" and feeds the users response into the variable 'id'. cono

[R] help

Hello there, How would I go about creating a function that will read a vector containing names of variables and then will output tables of proportions for all the categorical variables and summary statistics appropriate for all quantitative variables for all variables read by the function. Tha

[R] R ANOVA question

Hi All, I'm fitting an intercept ANOVA model using R. I'm using following contrast statement, so my first level is considered as base level. options(contrasts = c("contr.treatment", "contr.treatment")) My data has two factors are "Month" and "WeekDay". Due to default alphabetic ordering, the Mont

Re: [R] Error in eval(expr, envir, enclos) : object 'N' not found

Oh yes, that's true, but I don't understand why the same exact model works fine the first time I run LSODA, then running the model again with different PARMS causes it to error with this error. I've been trying to debug or find other relevant posts with this same error, but so far I have had no

Re: [R] R has stopped plotting

On Dec 6, 2009, at 11:47 AM, Walther, Alexander wrote: @ Ista: no, R still struggles with the output. Thanks anyway. :) Is there a way to merge multiple of these PNGs into a single PDF document afterwards? I would like to refrain from bridging this gap with additional software like f.i. Adobe

[R] citing R - really a bibtex/jabref question

Although this is really a bibtex/Jabref question, I am hoping this might still be the best place for help. The bibtex entry in Jabref is like this @MANUAL{R2009, title = {R: A Language and Environment for Statistical Computing}, author = {{R Development Core Team}}, organization = {R Founda

Re: [R] barchart() {lattice} help.

Please ignore my last question. I found a way to handle that. One last thing: I'm defining my own y scales. In the process the bar starts from below the "y=0" line (or below the y-axis). Is there a way to get rid of it. Here is the code and data I'm using R Code (Data read in object "dta") dta$

Re: [R] R has stopped plotting

@ Ista: no, R still struggles with the output. Thanks anyway. :) Is there a way to merge multiple of these PNGs into a single PDF document afterwards? I would like to refrain from bridging this gap with additional software like f.i. Adobe Acrobat. Cheers Alex Ista Zahn schrieb: > You need

Re: [R] barchart() {lattice} help.

Thank you Uwe, Dennis, and Gary for your help. I have one more question: I'm using pre-defined y-scales and trying to create grid lines. As "Female" category has low sum value, its y-axis range from 0-150 whereas "Male" ranges from 0-300. Is it possible to make them on same scale. Here is the pre

Re: [R] barchart() {lattice} help.

Thanks Uwe, I got your suggestions part too. 2009/12/6 Uwe Ligges > > > Peng Cai wrote: > >> Hi, >> >> I'm plotting grouped barplot using following code and data. I need help >> with >> re-ordering the labels. >> >> 1. On x-axis the factor "AGE" is grouped in order "0--4", "15--18", >> "5--14";

Re: [R] R has stopped plotting

You need dev.off() I think. -Ista On Sun, Dec 6, 2009 at 10:53 AM, Walther, Alexander wrote: > Dear list, > > as described beforehand, i wanted to merge three single plots into one > JPEG-file. Since I use PNG now, code goes as follows: > > > png(...) > > layout(...) > > par(...) > plot(...) > >

[R] R has stopped plotting

Dear list, as described beforehand, i wanted to merge three single plots into one JPEG-file. Since I use PNG now, code goes as follows: png(...) layout(...) par(...) plot(...) par(...) plot(...) par(...) plot(...) which works out pretty well, but now a new problem occurred: R has stopped t

Re: [R] Announcing a new R news site: R-bloggers.com

Tal, For completeness: On 5 December 2009 at 22:32, Tal Galili wrote: | 1) Planet R is limited (for years) to 26 feeds only, and I don't remember | seeing it evolve to include (or allow inclusion) of new R blogs that came | around. That was already addressed by Romain. It would appear that you

Re: [R] optim with constraints

Read the help page for `optim'. You already seem to be aware that "L-BFGS-B" should be used. How much more harder is to read the help page further to figure out how to supply the constraints? There are also various other algorithms, but "L-BFGS-B" should be able to do the job. Ravi. __

Re: [R] Announcing a new R news site: R-bloggers.com

Hi Romain. 1) Thank you for the link to how to add myself to PlanetR. 2) What I meant was that you don't need an online aggregator, since you know how to use it with your own local software. Cheers, Tal Contact Details:---

Re: [R] barchart() {lattice} help.

Peng Cai wrote: Hi, I'm plotting grouped barplot using following code and data. I need help with re-ordering the labels. 1. On x-axis the factor "AGE" is grouped in order "0--4", "15--18", "5--14"; whereas I would like to have it in "0--4", "5--14", "15--18". 2. If I need to re-order "RACE"

Re: [R] Error in eval(expr, envir, enclos) : object 'N' not found

hansoffate wrote: I'm running an LSODA to generate some graphs, but I need to stop at a certain point and use those values to generate another LSODA output. This is working fine, but when I try to run the second LSODA, I get the "Error in eval(expr, envir, enclos) : object 'N' not found". Any

Re: [R] MLE in R

Please read ?optim and about its arguments lower, upperBounds on the variables for the "L-BFGS-B" method. Uwe Ligges Liang Wang wrote: Hi, dear R users I am a newbie in R and I need to use the method of meximum likelihood to fit a Weibull distribution to my survival data. I use "optim"

Re: [R] Incorrect number of dimensions

Sarah wrote: I am trying to fit a linear model with seasonal effects but keep getting 'Incorrect number of dimensions'. I have no idea what this means or how to fix it. I am following instructions received from an instructor and it worked for him, so I assume it has something to do with my da

Re: [R] R packages and assess to data in packages

Sharpie wrote: Tariq Perwez wrote: Hi Everyone, I have two very basic questions and would appreciate your help. 1. I would like to see/access the data that comes with a given R package by using a function like read.table(). For example, I just installed car package from CRAN. I know that som

Re: [R] SAS "datalines" or "cards" statement equivalent in R?

Thanks Dennis and Jorge! On Sat, Dec 5, 2009 at 8:24 PM, Jorge Ivan Velez wrote: > Dear Gary, > > Here is a suggestion using read.table() and textConnection(): > > toread <- "id sex age inc r1 r2 r3 > 1 F 35 17 7 2 2 > 17 M 50 14 5 5 3 > 33 F 45 6 7 2 7 > 49 M 24 14 7 5 7 > 65 F

Re: [R] optim with constraints

As per ?optim Usage: optim(par, fn, gr = NULL, ..., method = c("Nelder-Mead", "BFGS", "CG", "L-BFGS-B", "SANN"), lower = -Inf, upper = Inf, control = list(), hessian = FALSE) Note that the optimx() function on R-forge http://r-forge.r-project.org/R/?grou

Re: [R] 4D plotting

3rd and 4th dimension can easily be accommodated in a 2d plot by using different color, shape of points etc. Therefore you can start with ggplot package. You might also want to look into here : http://had.co.nz/ggplot2/geom_point.html Best, Jared Nance wrote: > > Hello list, > Thanks in advanc

[R] 4D plotting

Hello list, Thanks in advance for whatever help you can give me. I have a data set that I want to visualize that has four important dimensions (x,y,z,t). What I would like to be able to do is plot the three spatial coordinates using, say, scatterplot3d (or similar), and have the color of the p

[R] MLE in R

Hi, dear R users I am a newbie in R and I need to use the method of meximum likelihood to fit a Weibull distribution to my survival data. I use "optim" as follows: optim(c(1.14,0.25),weibull.like,mydata=mydata,method="L-BFGS-B",hessian = TRUE) My question is: how do I setup the constraints tha

[R] optim with constraints

Hi, dear R users I am a newbie in R and I wantto use the method of meximum likelihood to fit a Weibull distribution to my survival data. I use "optim" as follows: optim(c(1, 0.25),weibull.like,mydata=mydata,method="L-BFGS-B",hessian = TRUE) My question is: how do I setup the constraints so that

[R] hi there!

hi there! how u doin?...hope fine. look, i've a difficulty in simulation especially for the bayesian analysis using R. Using the famous Markov Chain Monte Carlo (MCMC) method, Methropolis hasting algorithm i want to simuate different distrbution like: 1/ Weibull-exponential 2/Gamma-Gamma 3/Poisson-

Re: [R] Announcing a new R news site: R-bloggers.com

On 12/05/2009 11:45 PM, Tal Galili wrote: Hi Romain, First - I'd like to say that you where one of the bloggers I wanted to e-mail and didn't get to it. I still remember our short session in useR 2009 where you showed me your plans for an R editor - I am still waiting to see it once it comes out.

Re: [R] plot data from tapply

I left off the statement to load the reshape package. If you don't have it, install it from CRAN: x=c(1,2,3,1) y=c(1,2,3,1) ss=c(55,NA,55,88) ss_byxy_test=tapply( ss, list( x, y), mean, na.rm=TRUE) ss_byxy_test # use the 'reshape' package library(reshape) # now 'melt' the data to get it into a f

Re: [R] plot data from tapply

Here is one way of doing it: x=c(1,2,3,1) y=c(1,2,3,1) ss=c(55,NA,55,88) ss_byxy_test=tapply( ss, list( x, y), mean, na.rm=TRUE) # use the 'reshape' package ss_byxy_test # now 'melt' the data to get it into a format for plotting (ss_melt <- melt(ss_byxy_test)) # create the plot area so you can ad