[R] Est. Component Size with AIC/BIC under Gamma Distribution

Dear all, I am trying to model number of samples from a given series. The series are modelled according Gamma function. In order to estimate the # samples, I use BIC/AIC with MLE (computed from dgamma function). Here is the code I have. __BEGIN__ mlogl <- function( x_func, theta_func, samp) {

Re: [R] Where to download BRugs

install.packages("BRugs") still works on Windows -- the binaries have been moved to the CRAN extras collections. If you want a direct URL http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.7/BRugs_0.4-1.zip On Thu, 22 May 2008, qqw wrote: Hi all, I tried to follow an online tutorial

Re: [R] 1501-511 Compilation failed for file ch2inv.f - R on AIX 5.2.

Look like you need to use a Fortran compiler set to accept Fortran 77 (which the R sources are). We don't know what 'f95' is, and this is not really an R issue so please ask your sysadmins for help. On Thu, 22 May 2008, Senthilkumaran wrote: Hi, While executing the make after successful c

Re: [R] Pros and Cons of R

Completely agreed! BR, Shubha Shubha Karanth | Amba Research Ph +91 80 3980 8031 | Mob +91 94 4886 4510 Bangalore * Colombo * London * New York * San José * Singapore * www.ambaresearch.com -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Ruben Roa Ureta

Re: [R] Pros and Cons of R :GUI

R has very good GUI packages that I have used - R Commander RCmdr and Rattle rattle . Since I work on multiple packages with constraints of time, I almost always use the GUI rather go through the intricacies of command line . The log of these GUIs shows the relevant R command that was used, so you

[R] How to import package into R script

Hi all, What's the construct to import a built-in package (e.g. ARIMA) into the R script. I can't seem to find them in the documentation. Regards, Edward __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Where to download BRugs

Hi, I just installed it from http://cran.wustl.edu/bin/windows/contrib/2.7/BRugs_0.4-1.zip Change 2.7 to 2.6 and you get the older version. HTH, Aldi > > Thanks Aldi, do you know anywhere we could download previous package of > BRugs? Do you have a copy which you could send to me? (I only foun

Re: [R] Where to download BRugs

On 22/05/2008 8:35 PM, qqw wrote: Thanks Aldi, do you know anywhere we could download previous package of BRugs? Do you have a copy which you could send to me? (I only found one old version on Splus site which is not usable under R2.6 :() Thanks! You can get older versions from CRAN, in http:/

Re: [R] About Passing Arguments to Function

Dear Mark and Haris, Thanks so much to you both. I really appreciate it. - Edward On Fri, May 23, 2008 at 10:54 AM, <[EMAIL PROTECTED]> wrote: > to add to what Haris said: I think you are trying to find the mle's of the > gamma, given some data, and either what you're doing is not the right

Re: [R] Where to download BRugs

Thanks Aldi, do you know anywhere we could download previous package of BRugs? Do you have a copy which you could send to me? (I only found one old version on Splus site which is not usable under R2.6 :() Thanks! -- View this message in context: http://www.nabble.com/Where-to-download-BRugs-tp1

Re: [R] About Passing Arguments to Function

On May 22, 2008, at 9:14 PM, Edward Wijaya wrote: Hi, Below I have a function mlogl_k, later it's called with "nlm" . __BEGIN__ vsamples<- c(14.7, 18.8, 14, 15.9, 9.7, 12.8) mlogl_k <- function( k_func, x_func, theta_func, samp) { tot_mll <- 0 for (comp in 1:k_func) { curr_ml

[R] how to compute "explained variance" in metaMDS(vegan)?

Dear all, I used metaMDS (vegan) to generate two NMDS axis and now I would like to know the amount of variance explained for each axis. Is there a way to do that? kind regards, miltinho [[alternative HTML version deleted]] __ R-help@r-project

Re: [R] help with table

Is this what you want: > x <- by(frame, list(frame$p.a, frame$p.b), function(.row){ + c(f=sum(.row$r)/ nrow(.row), p.a=.row$p.a[1], p.b=.row$p.b[1],freq=nrow(.row)) + }) > do.call(rbind, x) f p.a p.b freq [1,] 0.50 10 22 [2,] 0.50 20 32 [3,] 0.25 50 54 On Thu, M

[R] About Passing Arguments to Function

Hi, Below I have a function mlogl_k, later it's called with "nlm" . __BEGIN__ vsamples<- c(14.7, 18.8, 14, 15.9, 9.7, 12.8) mlogl_k <- function( k_func, x_func, theta_func, samp) { tot_mll <- 0 for (comp in 1:k_func) { curr_mll <- (- sum(dgamma(samp, shape = x_func, scale=theta_

[R] Using unicode characters with paste or sprintf

Hi, I'm having difficulties with Unicode characters with R 2.7.0 on Windows XP SP2. I'd like to use a delta character, which I'd also like to mix with other formatted output using paste() or sprintf(). On my system, I can successfully do this with a static string: print("\u0394Q = 2.2 L/s")

Re: [R] mixed model results from SAS and R

Hi, To continue the test of the differences between mixed model SAS (9.1.3) and lme (R 7.0), I removed the intercept of the random effect in the R case. In that case lme is showing that beta coefficients are almost all of them identical (-0.08424262) and somehow all of them are spread along the

Re: [R] 3D Area Diagram

Thank you all for your help. I will check the packages you suggested asap and maybe I find a good way to get a such a plot. The picture I posted was only an example taken from the microsoft-website, just to make clear what kind of diagram I am looking for. I know that one part is hidden by the ot

Re: [R] Where to download BRugs

The BRugs package is maintained by Uwe Ligges. So it is possible he forgot to place it in the new version of R/ repositories. Aldi Charles Annis, P.E. wrote: Could you mean RBugs? Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3265 http://www.StatisticalEngineering.

Re: [R] Where to download BRugs

Could you mean RBugs? 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 Behalf Of qqw Sent: Thursday, May 22, 2008 3:10 PM To: r-help@r-project.org

[R] plot 7 * 3 matrix on DIN A4 pdf

Hello, I want to plot 21 scatter plots in a 7*3 matrix on a single A4 page (using mfcol). Below is an example, which -by and large- looks very close to the desired result. But I do not get it to fill out the whole page. If I specify the a different page size, the plots still do not fill out the wh

Re: [R] 3D Area Diagram

AlGates algates.de> writes: > > > Hello, > > maybe someone can help me. I am looking for a possibility to plot a 3D area > diagram like in Excel: > http://www.microsoft.com/germany/mac/excel/images/chartbefore.jpg Watch > this! > This would be theoretically possible, although not that eas

Re: [R] 3D Area Diagram

On May 22, 2008, at 3:04 PM, AlGates wrote: Hello, maybe someone can help me. I am looking for a possibility to plot a 3D area diagram like in Excel: http://www.microsoft.com/germany/mac/excel/images/chartbefore.jpg Watch this! Would be nice if someone had any idea about that. I don't

Re: [R] 3D Area Diagram

On 22/05/2008 3:04 PM, AlGates wrote: Hello, maybe someone can help me. I am looking for a possibility to plot a 3D area diagram like in Excel: http://www.microsoft.com/germany/mac/excel/images/chartbefore.jpg Watch this! Would be nice if someone had any idea about that. I'd say the world w

Re: [R] AMOVA results from ade4 different than in the reference publication

Hi, Ben is correct. The alternative hypothesis is not good for this test (within pop). I will modify it and submit a new version of the package on CRAN (probably next week). BTW, for question about ade4 please send your email to adelist (http://listes.univ-lyon1.fr/wws/info/adelist) Cheers

[R] 1501-511 Compilation failed for file ch2inv.f - R on AIX 5.2.

Hi, While executing the make after successful configuration (./configure --with-readline=no --without-iconv), the following error occurs: mbf2n11s (Regatta) /saswork/R/R-2.7.0$ make Target "R" is up to date. Target "R" is up to date. Target "R" is up to date. Target "R" is up to date. Target "Ma

[R] Where to download BRugs

Hi all, I tried to follow an online tutorial to run openBUgs but the package BRugs has been removed from R repository. Could someone provide a link for where to download BRugs? Thanks a lot! -- View this message in context: http://www.nabble.com/Where-to-download-BRugs-tp17411410p17411410.htm

[R] 3D Area Diagram

Hello, maybe someone can help me. I am looking for a possibility to plot a 3D area diagram like in Excel: http://www.microsoft.com/germany/mac/excel/images/chartbefore.jpg Watch this! Would be nice if someone had any idea about that. Thank you Alex -- View this message in context: http://www

[R] help with table

Dear All, I would be most grateful for help on this problem. I have a binary response variable. It is classified by various factors. I would like to have a table with; column 1 frequency response is 1 / (frequency the response is 1 + frequency the response is 0) #given the classifying factors

Re: [R] Alternatives to rJava and JRI

Although this is an update of the same code base, RServe is nice and has newer objects that one may use. For example, the exception handling has been improved and now there is a external server ability. Before you give up on JRI look into the newer code, it might be more useful. http://www.rforge.

Re: [R] Alternatives to rJava and JRI

There's SJava, http://www.omegahat.org/RSJava. On Thu, May 22, 2008 at 12:50 PM, Munir, Danish <[EMAIL PROTECTED]> wrote: > > Has anyone come across any alternatives to rJava and JRI? Are they any > good? Better perhaps? > > Please give your reasons. > > Thanks, > Danish > - - - - - - - - - - - -

Re: [R] Alternatives to rJava and JRI

I do not know if this is "better" as we are just starting to look into it (and "better" is in the eye of the beholder) but look at: http://tolstoy.newcastle.edu.au/R/e3/help/07/11/4253.html in the archives. The link to the software is http://www.ebi.ac.uk/ microarray-srv/frontendapp/ HTH,

[R] Alternatives to rJava and JRI

Has anyone come across any alternatives to rJava and JRI? Are they any good? Better perhaps? Please give your reasons. Thanks, Danish - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - This message is intended only for the personal and confidential use of t

Re: [R] length functions

Never mind! In fact I just discovered that the underscore and the "<-" are used equivalently in SPLUS. Applejus wrote: > > Hello, > > Does anyone know what the functions p_length and s_length do in SPLUS and > what is their equivalent in R? > > Thanks! > -- View this message in context:

Re: [R] Extracting slots from ROCR prediction objects

Hi Stacey, ROCR uses S4 classes. The elements are accessed using "@" instead of "$". You can find an example on slide 12 of the following slide deck: http://rocr.bioinf.mpi-sb.mpg.de/ROCR_Talk_Tobias_Sing.ppt Also have a look at the R code that appears when you type "demo(ROCR)" in R which contai

Re: [R] Pros and Cons of R

On Thu, May 22, 2008 at 02:07:01PM -0400, R P Herrold wrote: > On Thu, 22 May 2008, Monica Pisica wrote: [...] > >When a new R version is in place you > >cannot up-grade your old R version, you have to do a new > >installation and re-load all the packages you used to have > >and delete / un-in

Re: [R] setting specific colours for different groups of points on a scatter plot

On 5/22/08, Sarah Kinnings <[EMAIL PROTECTED]> wrote: > Hi, > > I have a table with three columns of data, let's call these a, b and c. I > have produced a scatter plot of a against b, but I would like to colour the > points according to their corresponding c values. c takes the value of > either

Re: [R] Pros and Cons of R

Hi, I have no idea about SAS but i know that when i've started using R i've struggled quite a bit until i started to get some insight and results. I cannot say it is easier now, but certainly i like it much better, i would not be ale to do the work i am doing without it and i would hate

Re: [R] Extracting slots from ROCR prediction objects

Hi Stacey, Try this: library(ROCR) data(ROCR.simple) pred <- prediction(ROCR.simple$predictions,ROCR.simple$labels) perf <- performance(pred,"sens","spec") cuts <- unlist(slot(perf,"alpha.values")) cuts HTH, Jorge On Thu, May 22, 2008 at 2:08 PM, Stacey Burrows <[EMAIL PROTECTED]> wrote: > Hi

Re: [R] Plotting a Quadratic...

Hi Ryan, See the example in ?predict.lm . I think it would be useful in this case. HTH, Jorge On Wed, Dec 31, 1969 at 8:16 PM, Ryan Harrigan <[EMAIL PROTECTED]> wrote: > I have an equation describing the best-fit model for a set of points (just > 2 > axes) that is in the form: > > y=b+mx+px^2

[R] Extracting slots from ROCR prediction objects

Hi, I have an object from the prediction function from the ROCR package and I would like to extract one of the slots from the object, for example the cutoffs slot. However the usual techniques ($, [["name"]]) of subsetting don't work. How can I assess the lists in the slots? Here is an example

[R] Pros and Cons of R

On Thu, 22 May 2008, Monica Pisica wrote: I don't consider packages as up-grades, i consider them as enhancements . When a new R version is in place you cannot up-grade your old R version, you have to do a new installation and re-load all the packages you used to have and delete / un-inst

Re: [R] Pros and Cons of R

There is a new book coming out this year on spatial data analysis in R. http://www.springer.com/statistics/statistical+theory+and+methods/book/978-0-387-78170-9 I am not sure there are any other book length treatments of spatial data analysis that utilize any of the other major statistical pa

Re: [R] Pros and Cons of R

On Thu, May 22, 2008 at 12:00 PM, Monica Pisica <[EMAIL PROTECTED]> 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

[R] Executing a command at interrupt

Colleagues, I am executing a large batch file in R (2.7.0) in OS X (Leopard) or Linux (Fedora). Usually, R completes normally and executes wrap-up commands. However, occasionally, R is "killed" externally, e.g., "kill " from the OS. When that happens, the wrap-up is not executed. Is

[R] Plotting a Quadratic...

I have an equation describing the best-fit model for a set of points (just 2 axes) that is in the form: y=b+mx+px^2 Where b is the intercept, m is the slope describing a linear term, and p is a slope of the quadratic term. I would like to plot this equation on a curve (I know the equation is y=(

[R] setting specific colours for different groups of points on a scatter plot

Hi, I have a table with three columns of data, let's call these a, b and c. I have produced a scatter plot of a against b, but I would like to colour the points according to their corresponding c values. c takes the value of either 0,1,2,3 or 4 and so I would like to colour all of the poi

Re: [R] Pros and Cons of R

Hi, IMHO, "- Everything you want to do is a command line, minimal GUI", this should be considered a "Pro". The commandline is much more powerful than a GUI. Cheers, John Monica Pisica wrote: > > Charles, > > I don't consider packages as up-grades, i consider them as enhancements . >

Re: [R] Pros and Cons of R

I think it is a fairly safe bet that at any given time there are numerous 64-bit R sessions using more than 4GB of memory -- I periodically have one. In my opinion 64-bit implementation is definitely in the "pro" column. It is as transparent going from 32-bit to 64-bit as moving between 32-bit pl

Re: [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 > h

Re: [R] AMOVA results from ade4 different than in the refer ence publication

Kim Milferstedt uiuc.edu> writes: > > Hello, > > I am trying to run some AMOVA analyses with the amova function in the > package ade4. > > When running the example dataset provided in ade4, I noticed a > difference between the published results from the same data (Excoffier > et al. 1992) a

Re: [R] Pros and Cons of R

Hi, Monica Pisica wrote: > - There is no perfect “beginner” book. How about - Crawley, Michael (2007). The R book, Wiley & Sons. - Maindonald, John & John Braun (2007): Data Analysis and Graphics Using R (2nd edition), Cambridge University Press. As a political scientist (with programming exp

Re: [R] Pros and Cons of R

Thanks for the heads up - i will add that to the "pros". Monica> From: [EMAIL PROTECTED]> To: [EMAIL PROTECTED]; r-help@r-project.org> Date: Thu, 22 May 2008 09:32:06 -0700> Subject: RE: [R] Pros and Cons of R> > On the 64-bit part, I tested 2.7.0 on a Dual Core Lenovo ThinkPad and was able t

Re: [R] Pros and Cons of R

Charles, I don't consider packages as up-grades, i consider them as enhancements . When a new R version is in place you cannot up-grade your old R version, you have to do a new installation and re-load all the packages you used to have and delete / un-install the old version of course

[R] AMOVA results from ade4 different than in the reference publication

Hello, I am trying to run some AMOVA analyses with the amova function in the package ade4. When running the example dataset provided in ade4, I noticed a difference between the published results from the same data (Excoffier et al. 1992) and what ade4 calculates. Below are the data for "wi

[R] new mailing list: R-sig-genetics

Dear list, A new mailing list *R-sig-genetics* has just been created (thanks to Martin Maechler). This list is dedicated to discussions about the *analysis of genetic markers* using the R software. The covered topics, along with various information about the list, are available from the R-si

Re: [R] Pros and Cons of R

This is about the 8761st time that this old warhorse has been discussed on this list. I suggest you search the archives (e.g. via CRAN, Gmane, Google,...) for these prior discussions. -- Bert Gunter Genentech P.S. Incidentally, the answer to your question is "It depends..." -- on what your appli

[R] Packfor library-problem downloading

Hello, I'm a new user to R, and I was trying to install "packfor" in my R console and I always run into an error. I tried to search, in in the forum, someone that already might have made the same question, but all I could find was a general information how to install libraries in R. Pass through

Re: [R] Binning groups

On Thu, 22 May 2008, francogrex wrote: f > f > Hi, this is probably quite simple but I can't seem to do it correctly. I have f > a data frame of counts of infections in different ages; something like: f > count=c(1, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 7, 8, 8, 9, 9, f > 10, 11, 15, 17,

Re: [R] disaggregate frequency table into flat file

Marc, it's the second "expansion" type transformation I was after, although your expand.dft looks quite complicated? here's what I finaly came up with - the bold lines correspond to what expand.dft does? > orig<-matrix(c(40,5,30,25), c(2,2)) > orig [,1] [,2] [1,] 40 30 [2,]5 25 >

Re: [R] Pros and Cons of R

On the 64-bit part, I tested 2.7.0 on a Dual Core 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

Re: [R] Plotting 3 Time Series

Hi ngupta: i didn't get what you mean with "If B than xts..If T yts..." i suppose that you want to plot the prices of type A, then the prices of type B and so on... # well, to extract prices of tipe A, you can use A <- which(table[,1] == "A") # and then A.prices <- table[,2][A] hope that helps,

Re: [R] Pros and Cons of R

What do you mean, "there are no up-grades?" There are 1,401 ancillary packages - all free. That sounds like upgrades to me. Of course there is no *perfect* beginner's book, but Peter Dalgaard's Introductory Statistics with R (Paperback), Springer, 3d printing edition (January 9, 2004) is pretty

Re: [R] Pros and Cons of R

Hello, I enjoy using R, and find the excellence of this newsgroup removes most of the altitude of the "learning curve". Any learning curve, IMHO, is more than offset by the benefits and capabilities of the language and its awesome plotting tools. A grateful advocate, John Kevin E. Thorpe wro

Re: [R] Pros and Cons of R

Wensui Liu wrote: agree, i think R is more like a standard program language than SAS. however, SAS programmer might not feel intuitive to pick up R. That says more about SAS programmers than it does about R. __ R-help@r-project.org mailing list htt

Re: [R] Pros and Cons of R

agree, i think R is more like a standard program language than SAS. however, SAS programmer might not feel intuitive to pick up R. On Thu, May 22, 2008 at 12:14 PM, Kevin E. Thorpe <[EMAIL PROTECTED]> wrote: > Monica Pisica wrote: >> >> Cons: >> >> - R has a very steep learning curve. > > I don't

Re: [R] Pros and Cons of R

Monica Pisica wrote: > > Cons: > > - R has a very steep learning curve. I don't think the learning curve is any steeper than SAS programming, it is just a different kind of curve. -- Kevin E. Thorpe Biostatistician/Trialist, Knowledge Translation Program Assistant Professor, Department of Pub

Re: [R] Pros and Cons of R

On Thu, May 22, 2008 at 5:00 PM, Monica Pisica <[EMAIL PROTECTED]> wrote: > - It seems that it is hard to install correctly under Linux. Actually, it is quite easy to install R under Linux, at least in some distributions. For instance, on Fedora: yum install R R-devel and that is it. Paul

[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] R static is dynamically linked!!

On Tue 20 May 2008 at 11:52:32 AM -0400, George Georgalis wrote: > >But this makes me realize, since our OS is generally >backwards compatible, I may be able to build on the old >systems and use on the newer ones as well. Will give >that a shot. > that doesn't work, I think because gcc changed fro

[R] length functions

Hello, Does anyone know what the functions p_length and s_length do in SPLUS and what is their equivalent in R? Thanks! -- View this message in context: http://www.nabble.com/length-functions-tp17406498p17406498.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] How to use classwt parameter option in RandomForest

The weights given should correspond to the ordering of levels(y) where y contains the class labels. If in doubt, you can also give the classwt as a named vector (e.g., classwt=c(B=3, A=2, C=1)). Search in the R-help archive to see other options and why you probably shouldn't use classwt. Andy F

[R] Plotting 3 Time Series

Hi, I am fairly new to R so this may be easy for some. I have one column that is the Type="A, B' T" and I have another column that lists Price='894, 895, 896' I wanted to plot three different series based upon the Type...If B than xts..If T yts... How exactly do I do this? Thanks in advance --

Re: [R] automation of R? running an R script at a certain time eachnight?

I do this all the time. See http://finzi.psych.upenn.edu/R/Rhelp02a/archive/105099.html (Vista is a bit harder than XP, but similar.) As far as determining when to proceed with the script, you could have your db updater create a semaphore file your batch or R script could wait for. HTH, David L.

Re: [R] disaggregate frequency table into flat file

Is this what you want? > xtabs(Freq ~ Var1 + Var2, data = orig) Var2 Var1 A B A 40 30 B 5 25 See ?xtabs Or is this what you want? expand.dft <- function(x, na.strings = "NA", as.is = FALSE, dec = ".") { DF <- sapply(1:nrow(x), function(i) x[rep(i, each = x$Freq[i]), ],

Re: [R] disaggregate frequency table into flat file

On May 22, 2008, at 8:56 AM, maiya wrote: sorry, my mistake! the data frame should read: orig<-as.data.frame.table(orig) orig Var1 Var2 Freq 1AA 40 2BA5 3AB 30 4BB 25 but basicaly i would simply like a sample of the original matrix ( which is a fre

Re: [R] disaggregate frequency table into flat file

sorry, my mistake! the data frame should read: orig<-as.data.frame.table(orig) orig Var1 Var2 Freq 1AA 40 2BA5 3AB 30 4BB 25 but basicaly i would simply like a sample of the original matrix ( which is a frequency table/contingency table/crosstabulation)

Re: [R] Stripchart and Boxplots side-by-side

Dieter Wirz: > with(InsectSprays, boxplot(count ~ spray, boxwex = 0.3)) > with(InsectSprays, stripchart(count ~ spray, col = "red", vertical = > TRUE, add = TRUE)) > > But the dots from the stripchart are plotted over the Boxes. > Is there any possibility to have Stripchart and Boxplots side-by-s

Re: [R] disaggregate frequency table into flat file

Not exactly clear what you are asking for. Your data.frame.table does not seem related to the original 'orig'. What exactly are you expecting as output? On Wed, May 21, 2008 at 10:16 PM, maiya <[EMAIL PROTECTED]> wrote: > > i appologise for the trivialness of this post - but i've been searching

Re: [R] secondary axis in a plot

If you need to do it: plot(PLOT1, bty='c') par(new=TRUE) plot(PLOT2, axes=FALSE, ylab='', xlab='', bty='c') axis(4) mtext("secondary y-axis label",4) On Thu, May 22, 2008 at 7:16 AM, MarySR <[EMAIL PROTECTED]> wrote: > > Hello, > I want to do a plot of two series, but adding a secondary axis in

Re: [R] secondary axis in a plot

MarySR wrote: Hello, I want to do a plot of two series, but adding a secondary axis in the rigth side to take into account the different scales of the two series. Anyone can tell me haw can I do that in R? Hi Maria, Have a look at twoord.plot in the plotrix package. Jim _

Re: [R] secondary axis in a plot

MarySR: > I want to do a plot of two series, but adding a secondary axis in the > rigth side to take into account the different scales of the two series. > Anyone can tell me haw can I do that in R? I would *strongly* advice you to read this article on dual-scaled axes before proceeding: http://

[R] Stripchart and Boxplots side-by-side

Dear all - With the following code I get Boxplots and Stripcharts in one Plot: with(InsectSprays, boxplot(count ~ spray, boxwex = 0.3)) with(InsectSprays, stripchart(count ~ spray, col = "red", vertical = TRUE, add = TRUE)) But the dots from the stripchart are plotted over the Boxes. Is there any

[R] Stripchart and Boxplots side-by-side

Dear all - With the following code I get Boxplots and Stripcharts in one Plot: with(InsectSprays, boxplot(count ~ spray, boxwex = 0.3)) with(InsectSprays, stripchart(count ~ spray, col = "red", vertical = TRUE, add = TRUE)) But the dots from the stripchart are plotted over the Boxes. Is there any

[R] deconvoluting hclust objects

I have a hclust object that looks like this: > test77 Call: hclust(d = input) Cluster method : complete Number of objects: 11 > test77$height [1] 0.000 0.000 0.000 0.000 0.000 0.000 0.900 [8] 0.9473684 1.7894737 8.5771948 > test77$merge [,1] [,2] [1,] -

[R] secondary axis in a plot

Hello, I want to do a plot of two series, but adding a secondary axis in the rigth side to take into account the different scales of the two series. Anyone can tell me haw can I do that in R? Thank you very much Maria -- View this message in context: http://www.nabble.com/secondary-axis-in-a-plo

[R] mixed model resuts from SAS and R

Hi, I was wondering if there is a way to figure out why in SAS random beta coefficients are 0 vs. in R the beta-s are non zero. The variables of the data are nidl, time, and sub (for subject). Time and nidl are continuous variables. I am applying random coefficients model. Any input is greatly app

Re: [R] Binning groups

If you just want to bin the current ages: > count=c(1, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 7, 8, 8, 9, 9, + 10, 11, 15, 17, 17, 17, 17, 19, 19, 19, 19, 20, 20, 20, 21, 21, + 21, 22, 22, 22, 22, 23, 23, 23, 23, 23, 24, 24, 25, 27, 31, 33 + ) > age=c(3, 8, 9, 7, 56, 58, 10, 13, 53, 55, 11, 12,

Re: [R] How to account for autoregressive terms?

Please do not send HTML (as the posting guide asks) -- it has made your code unreadable. Also, without seeing an example we don't know if the 'completely different' is in the output or corresponds to very different fitted models. There is no 'armaFit' in R: it seems you may be using package

Re: [R] Binning groups

Is this what you want: > count=c(1, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 7, 8, 8, 9, 9, + 10, 11, 15, 17, 17, 17, 17, 19, 19, 19, 19, 20, 20, 20, 21, 21, + 21, 22, 22, 22, 22, 23, 23, 23, 23, 23, 24, 24, 25, 27, 31, 33 + ) > age=c(3, 8, 9, 7, 56, 58, 10, 13, 53, 55, 11, 12, 14, 51, 54, + 15,

[R] How to account for autoregressive terms?

Hi, how to estimate a the following model in R: y(t)=beta0+beta1*x1(t)+beta2*x2(t)+...+beta5*x5(t)+beta6*y(t-1)+beta7*y(t-2)+beta8*y(t-3) 1) using "lm" : dates <- as.Date(data.df[,1]) selection<-which(dates>=as.Date("1986-1-1") & dates<=as.Date("2007-12-31")) dep <- ts(data.df[selection,c("dep

Re: [R] help (boot)

On 5/21/2008 4:52 PM, Sharai Gomez-Raygoza wrote: Hi!!! I am trying to run a bootstrap but I have troubles. My code is as follows theta <- function(x){mean(x,)} boot1<-boot(x,theta,R=1000) Error in statistic(data, original, ...) : unused argument(s) (1:487) Any help is welcome Thanks When s

Re: [R] Cumulative average

does e.g. x <- c(8,5,4,2) cumsum(x) / 1:length(x) do what you need? Albert Am Mittwoch, den 21.05.2008, 20:36 -0500 schrieb Jacques Wagnor: > Dear List, > > Does there exist a function that calculates a cumulative average? > Neither running() from library(gregmisc) nor running.mean() from > libra

[R] Binning groups

Hi, this is probably quite simple but I can't seem to do it correctly. I have a data frame of counts of infections in different ages; something like: count=c(1, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 7, 8, 8, 9, 9, 10, 11, 15, 17, 17, 17, 17, 19, 19, 19, 19, 20, 20, 20, 21, 21, 21, 22, 22, 22

[R] disaggregate frequency table into flat file

i appologise for the trivialness of this post - but i've been searching the forum wothout luck - probably simply because it's late and my brain is starting to go.. i have a frequency table as a matrix: orig<-matrix(c(40,5,30,25), c(2,2)) orig [,1] [,2] [1,] 40 30 [2,]5 25 i basic

[R] Cumulative average

Dear List, Does there exist a function that calculates a cumulative average? Neither running() from library(gregmisc) nor running.mean() from library(igraph) seems to be able to give a cumulative average. Any help or pointers would be greatly appreciated. Jacques ___

[R] help (boot)

Hi!!! I am trying to run a bootstrap but I have troubles. My code is as follows theta <- function(x){mean(x,)} boot1<-boot(x,theta,R=1000) Error in statistic(data, original, ...) : unused argument(s) (1:487) Any help is welcome Thanks [[alternative HTML version deleted]]