[R] package license questions

I am looking for some advice on licenses. Here is my situation: Over the last couple years, I have developed a rather large number of fire department analysis functions. I am in the process of trying to publish some packages to make these functions available to the public. I am trying to

Re: [R] caculate the frequencies of the Amino Acids

I sent a private reply to the poster yesterday with the solution to his problem, bit evidently he hasn't seen it, so ... # Function to tabulate an amino acid sequence # using strsplit and table: AAtab2 - function(x) { # Creates a frequency table of the characters in a string AA -

Re: [R] Run R in server web

Hi Jindo, I don't know how to do what you ask, but I recommend you to have a look at this: http://www.r-statistics.com/2009/12/a-web-application-of-rs-ggplot2/ And try to follow the links to the websites behind it, maybe somewhere in it there will be more information. Best, Tal

[R] questions on generic functions

Hi all, I'm trying to build a package, but run into problems with generic functions. I have two classes aa and bb, and defined the function summary.bb. Now I want to make summary.bb the default method for both classes aa and bb. That is, if I have two objects object.aa and object.bb from

Re: [R] Help with function fitdistr in MASS

Please read the footer of this message. ?fitdistr says x: A numeric vector. and setting the class does not make it a numeric vector (it is just a label). And fitdistr early on does if (missing(x) || length(x) == 0L || mode(x) != numeric) stop('x' must be a non-empty

Re: [R] xyplot: problems with column names legend

Hi, Using backticks might work to some extent, library(lattice) `my variable` = 1:10 y=rnorm(10) xyplot(`my variable` ~ y) but if your data is in a data.frame the names should have been converted, make.names('my variable') [1] my.variable HTH, baptiste 2010/1/3 Jay josip.2...@gmail.com:

Re: [R] Help with function fitdistr in MASS

On Sat, 2010-01-02 at 23:20 -0800, Saji Ren wrote: Hi, R users: I want to fit my data into a normal distribution by using the command fitdistr in MASS. I changed my data class from ts to numeric by class(mydata)=numeric but after using fitdistr, I got the result below

Re: [R] Rscript: how to suppress all output

Thanks! I'm using Mac 10.6.2. This did it: ~/Documents cat test.r #! /usr/bin/Rscript --vanilla --default-packages=utils suppressMessages(require(Hmisc)) suppressMessages(require(gsubfn)) a - 1:100 print(a) ~/Documents Rscript test.r /dev/null ~/Documents -J 2010/1/3 Prof Brian Ripley

Re: [R] [Off-topic] problem with Tinn-R editor

Ask directly to Jose-Claudio Faria (jcfa...@users.sourceforge.net). He is the author and maintainer of Tinn-R. Best, Philippe Grosjean ..°})) ) ) ) ) ) ( ( ( ( (Prof. Philippe Grosjean ) ) ) ) ) ( ( ( ( (Numerical Ecology of Aquatic

Re: [R] Anova in 'car': SSPE apparently deficient rank

Colleen F. Moore wrote: I have design with two repeated-measures factor, and no grouping factor. I can analyze the dataset successfully in other software, including my legacy DOS version BMDP, and R's 'aov' function. I would like to use 'Anova' in 'car' in order to obtain the sphericity

Re: [R] Please help me!!!! Error in `[.data.frame`(x, , retained, drop = FALSE) : undefined columns selected

thanks ，now i reduce the number of the descriptors, It is ok. 在2010-01-03，Max Kuhn [via R] ml-node+997387-1292262...@n4.nabble.com 写道： -原始邮件- 发件人:Max Kuhn [via R] ml-node+997387-1292262...@n4.nabble.com 发送时间:2010年1月3日 星期日

[R] RHTMLForms Package

Can anyone tell me from where to download the package RHTMLForms ? CRAN seems not holding this anymore. I am using Windows Vista OS. Thanks -- View this message in context: http://n4.nabble.com/RHTMLForms-Package-tp997672p997672.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] RHTMLForms Package

On Sun, 3 Jan 2010, Ron_M wrote: Can anyone tell me from where to download the package RHTMLForms ? CRAN seems not holding this anymore. I am using Windows Vista OS. It never did: it is an Omegahat package, see http://www.omegahat.org/RHTMLForms/. (I don't understand why a Google search on

Re: [R] xyplot: problems with column names legend

Thanks, the backtickes got the code working. However, now I cant get it to draw the legend/key. For example, look at this figure: http://osiris.sunderland.ac.uk/~cs0her/Statistics/xyplot5.png My graph is similar, but instead of 1,2,...,8 as the names of the series I want it to say Data one (a

Re: [R] package license questions

On 02/01/2010 11:19 PM, Markus Weisner wrote: I am looking for some advice on licenses. Here is my situation: Over the last couple years, I have developed a rather large number of fire department analysis functions. I am in the process of trying to publish some packages to make these

[R] Scores in Dimensions in Correspondence Analysis

Hi R, I was experimenting with the CA package for correspondence analysis. data(author) ca(author) plot(ca(author)) How can I get the values used by the plot for the two axis for biplot? In other words, how will I be able to obtain the scores in the dimensions (which are obtained

Re: [R] Scores in Dimensions in Correspondence Analysis

Does the $rowcoord and $colcoord help? Thanks and Regards, Shubha Karanth -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Shubha Vishwanath Karanth Sent: Sunday, January 03, 2010 6:25 PM To: r-help@r-project.org Subject: [R]

[R] F77_CALL, F77_NAME definition

I give up. Maybe it is my search (Windows) but I cannot seem to find the definition of the F77_CALL or F77_NAME macros. Either there are too many matches or the search just doesn't find it. For example where is the source for: F77_CALL(dpotri) ? Thank you. Kevin

Re: [R] F77_CALL, F77_NAME definition

On 03/01/2010 9:11 AM, rkevinbur...@charter.net wrote: I give up. Maybe it is my search (Windows) but I cannot seem to find the definition of the F77_CALL or F77_NAME macros. Either there are too many matches or the search just doesn't find it. For example where is the source for:

Re: [R] Scores in Dimensions in Correspondence Analysis

Got that... Thanks and Regards, Shubha Karanth | Amba Research Ph +91 80 3980 8283 | Mob +91 94 4886 4510 Bangalore * Colombo * London * New York * San José * Singapore * www.ambaresearch.com   -Original Message- From: r-help-boun...@r-project.org

Re: [R] F77_CALL, F77_NAME definition

On Jan 3, 2010, at 9:11 AM, rkevinbur...@charter.net wrote: I give up. Maybe it is my search (Windows) but I cannot seem to find the definition of the F77_CALL or F77_NAME macros. Either there are too many matches or the search just doesn't find it. For example where is the source for:

Re: [R] F77_CALL, F77_NAME definition

On Sun, Jan 3, 2010 at 2:11 PM, rkevinbur...@charter.net wrote: I give up. Maybe it is my search (Windows) but I cannot seem to find the definition of the F77_CALL or F77_NAME macros. Either there are too many matches or the search just doesn't find it. For example where is the source for:

Re: [R] questions on generic functions

Wayne (Yanwei) Zhang wrote: Hi all, I'm trying to build a package, but run into problems with generic functions. I have two classes aa and bb, and defined the function summary.bb. Now I want to make summary.bb the default method for both classes aa and bb. If you are using S3 approach,

Re: [R] F77_CALL, F77_NAME definition

On 03/01/2010 9:43 AM, Barry Rowlingson wrote: On Sun, Jan 3, 2010 at 2:11 PM, rkevinbur...@charter.net wrote: I give up. Maybe it is my search (Windows) but I cannot seem to find the definition of the F77_CALL or F77_NAME macros. Either there are too many matches or the search just doesn't

Re: [R] F77_CALL, F77_NAME definition

On Sun, Jan 3, 2010 at 3:20 PM, Duncan Murdoch murd...@stats.uwo.ca wrote: I think it's better to use a reasonable text editor here; I used Textpad.  I don't think there's anything too special about it, but it does have Search | Find in files, and I can list the file pattern (obviously *.h for

Re: [R] F77_CALL, F77_NAME definition

On 03/01/2010 10:32 AM, Barry Rowlingson wrote: On Sun, Jan 3, 2010 at 3:20 PM, Duncan Murdoch murd...@stats.uwo.ca wrote: I think it's better to use a reasonable text editor here; I used Textpad. I don't think there's anything too special about it, but it does have Search | Find in files,

Re: [R] Anova in 'car': SSPE apparently deficient rank

Dear Peter and Colleen, I think that Peter realizes this, but what Anova() does in this case is equivalent to the MANOVA anova(multmodel, M = ~ Afac*Bfac, X = ~Afac + Bfac, idata=poke.idata) Error in anova.mlm(multmodel, M = ~Afac * Bfac, X = ~Afac + Bfac, idata = poke.idata) : residuals

[R] update packages from local

Hi all, I have an old package installed, say abc. Now I made some changes to the source, and built a new version of the source code abc_0.1.1.tar.gz. How can I update the old package to this newer version from the local tar.gz file? I was running the following, but it did not work.

Re: [R] update packages from local

On Jan 3, 2010, at 11:08 AM, Wayne (Yanwei) Zhang wrote: Hi all, I have an old package installed, say abc. Now I made some changes to the source, and built a new version of the source code abc_0.1.1.tar.gz. How can I update the old package to this newer version from the local tar.gz file?

Re: [R] [R-sig-eco] How to create axes in arbitrary positions?

Dear Zongshan Li, I'm forwarding this back to r-help ... it's generally best to keep the conversation going in public, so that (a) other people can chime in with ideas and suggestions, (b) your question doesn't get lost if the original respondent doesn't have time to deal with it, (c) the

[R] R2HTML Report number format, or Better Way?

Here I am again with question I'll feel foolish for asking, when I see the answer. I'm trying to produce a report and here's where I get stuck: How do I get R2HTML to produce the same number format? Particularly remove the decimal places for Par and Sal. Are there better methods to produce

Re: [R] Ordering variables in a parallel coordinates plot

On Sat, 2 Jan 2010, Tal Galili wrote: Hello all, I am searching for a way in R to re-order variables before presenting them in a parallel coordinates plot. So far I didn't find anything within a R related context on how to do this. I did find some texts talking about how it should be done in

[R] problems about Quantile Regression for Longitudinal Data

Hello, I uesd Dr. Koenker's code on the following website for quantile regression of longitudinal data: http://www.econ.uiuc.edu/~roger/research/panel/rq.fit.panel.R I can run the code when there are three explanatory variables in X matrix (ncol(X)=4, because there is a vector 1 at the first

Re: [R] update packages from local

Thanks. But I after I restarted R, it is still the old version of the package that is loaded. Am I using the update.packages() correctly? Regards, Wayne (Yanwei) Zhang __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] update packages from local

1. Please quote the old threat, at least I do not keep old messages. 2. Try install.packages() from a freshly started R (that has not loaded the package in advance). 3. Since I am doing things in the shell for package development anyway, I tend to call package installation from the shell as

Re: [R] update packages from local

On Jan 3, 2010, at 12:42 PM, Wayne (Yanwei) Zhang wrote: Thanks. But I after I restarted R, it is still the old version of the package that is loaded. Am I using the update.packages() correctly? Apparently not. I'm not an expert at this so it's perfectly possible that someone smarter

[R] ADF test : how to deal with the missing values ?

[[alternative(swapped) HTML version deleted]] This is MIME Epilogue __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

[R] ADF test : how to deal with the missing values ?

hi everyone, and happy new year. i want to perform an Augmented Dickey Fuller test. the basic adf.test function doesn't work on my dataset as data are missing. has a package a function enabling missing values ? does any one know a trick to get through ? thanks, édouard.

[R] bio3d package not installing

Is it just me, or is the bio3d package no longer available? Is there another way to do a reasonable nucleic acid sequence alignment in R? Thanks, Jonathan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] update packages from local

Sorry. The original message is below. Uwe, do you mean if i just want to update the old package abc, i have to run R INSTALL abc_0.1.1.tar.gz in the shell? I know that will overwrite that package, but can I do something with the update.packages() function within R? From what I read in the

Re: [R] ADF test : how to deal with the missing values ?

On Jan 3, 2010, at 1:27 PM, Edouard Tallent wrote: hi everyone, and happy new year. i want to perform an Augmented Dickey Fuller test. the basic adf.test function ... from package? Yeah, I know I could look it up but it is considered good posting form to mention the package name.

Re: [R] Ordering variables in a parallel coordinates plot

Hi Charles, Thanks for answering - you are right about the posting guide (sorry). In any case, my question wasn't on how to reorder the columns (that is simple), but on how to choose what order to put them in. Thanks, Tal Contact

Re: [R] Anova in 'car': SSPE apparently deficient rank

John Fox wrote: Dear Peter and Colleen, I think that Peter realizes this, but what Anova() does in this case is equivalent to the MANOVA anova(multmodel, M = ~ Afac*Bfac, X = ~Afac + Bfac, idata=poke.idata) Error in anova.mlm(multmodel, M = ~Afac * Bfac, X = ~Afac + Bfac, idata =

[R] plot question

Hi, I am new to R so forgive me if the following query is somewhat simple.  I have a small tab-separated file with n records of glucose values, one record per day, seven measurements per day. It looks like this: date    sober    no    vm    nm    va    na    vs 20091229    NA    6.8    NA   

Re: [R] Ordering variables in a parallel coordinates plot

Hi Charles, You're solution is great (and is actually what my professor suggested me to do today). In the meantime I searched even more and found this article: http://davis.wpi.edu/~xmdv/docs/tr0313_osf.pdf That gives a good description of the problem and of his attempts at solving it. I

Re: [R] update packages from local

On Sun, 3 Jan 2010, Wayne (Yanwei) Zhang wrote: Sorry. The original message is below. Uwe, do you mean if i just want to update the old package abc, i have to run R INSTALL abc_0.1.1.tar.gz in the shell? I know that will overwrite that package, but can I do something with the update.packages()

[R] Questions regarding sem using hetcor() function from polycor and diagrams

Hello R Users, While I have attempted to dig into the R help files and I have not identified the answer to these questions, I apologize in advance if my questions were answered in the past. I also recognize that one of my questions unfortunately verges on statistical rather than code

Re: [R] bio3d package not installing

Is this link of use: http://mccammon.ucsd.edu/~bgrant/bio3d/download/download.html -Peter Ehlers Jonathan wrote: Is it just me, or is the bio3d package no longer available? Is there another way to do a reasonable nucleic acid sequence alignment in R? Thanks, Jonathan [[alternative

[R] Artificial Neural Networks

Hi! I am studying to use some R libraries which are applied for working with artificial neural neworks (amore, nnet). Can you recommend some useful, reliable and easy to get example data to use in R for creating and testing a neural network? And what library will you advise?

Re: [R] Questions regarding sem using hetcor() function from polycor and diagrams

Dear Stella, -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Stella Copeland Sent: January-03-10 3:40 PM To: r-help@r-project.org Subject: [R] Questions regarding sem using hetcor() function from polycor and diagrams Hello

Re: [R] plot question

I gather that each row defines one series, right? In that case try this: Lines - datesobernovmnmvanavs 20091229NA6.8NA2.711.72.76.2 20091230NA6.8NA2.711.72.76.2 20091231NA6.8NA2.711.72.7

Re: [R] Ordering variables in a parallel coordinates plot

On Sun, 3 Jan 2010, Tal Galili wrote: Hi Charles, Thanks for answering - you are right about the posting guide (sorry). In any case, my question wasn't on how to reorder the columns (that is simple), but on how to choose what order to put them in. OK, and I see you included a reference, so

Re: [R] Ordering variables in a parallel coordinates plot

Hi, Tal, Tal Galili wrote: In any case, my question wasn't on how to reorder the columns (that is simple), but on how to choose what order to put them in. So is your question answered? Or what exactly is the problem using the example below? library(lattice) parallel(~iris[1:4] |

[R] Thin bars in R hist !!??

Hi, I am trying to plot a histogram with my dataset that has 68 elements, 67 of which are zero and the last one is 18. It can be reproduced as follows: x-array(0, dim=(68)) x[1] = 18 I am trying to plot its histogram using: hist(x, breaks=10, xlim=c(0, 100), axes=F) axis(2, at=seq(0, 70, 5))

Re: [R] R2HTML Report number format, or Better Way?

On 01/04/2010 03:36 AM, L.A. wrote: Here I am again with question I'll feel foolish for asking, when I see the answer. I'm trying to produce a report and here's where I get stuck: How do I get R2HTML to produce the same number format? Particularly remove the decimal places for Par and Sal.

Re: [R] R2HTML Report number format, or Better Way?

L.A. wrote: Here I am again with question I'll feel foolish for asking, when I see the answer. I'm trying to produce a report and here's where I get stuck: srtype-cbind(Par,Sal,Median,COD,PRD,LowerCI,UpperCI) srtype Chances are better to get a reply when you supply the data as

Re: [R] Thin bars in R hist !!??

On 01/04/2010 08:25 AM, Whalojazz wrote: Hi, I am trying to plot a histogram with my dataset that has 68 elements, 67 of which are zero and the last one is 18. It can be reproduced as follows: x-array(0, dim=(68)) x[1] = 18 I am trying to plot its histogram using: hist(x, breaks=10, xlim=c(0,

Re: [R] Thin bars in R hist !!??

Just seting breaks=2 should fix your 'issue'. But what could such a histogram possibly tell anyone? Unless this is a pathological case in a more elaborate analysis, I can't see any sense in what you're doing. -Peter Ehlers Whalojazz wrote: Hi, I am trying to plot a histogram with my dataset

Re: [R] Thin bars in R hist !!??

Whalojazz wrote: Hi, I am trying to plot a histogram with my dataset that has 68 elements, 67 of which are zero and the last one is 18. It can be reproduced as follows: x-array(0, dim=(68)) x[1] = 18 I am trying to plot its histogram using: hist(x, breaks=10, xlim=c(0, 100), axes=F) axis(2,

Re: [R] Artificial Neural Networks

Hi Alex, I am not a NN user, but I would suggestion you give more information about what are you planning to do. May be it helps the helper. Happy New Year. miltinho On Sun, Jan 3, 2010 at 3:54 PM, Alex Olafson alex.olaf...@yahoo.com wrote: Hi! I am studying to use some R libraries which are

Re: [R] Thin bars in R hist !!??

Peter Ehlers wrote: Just seting breaks=2 should fix your 'issue'. But what could such a histogram possibly tell anyone? Unless this is a pathological case in a more elaborate analysis, I can't see any sense in what you're doing. Yes it did! However, this requires that I explicitly set

Re: [R] Thin bars in R hist !!??

Peter Dalgaard wrote: Were you perhaps looking for hist(x, breaks=seq(0, 100,10))? (Notice that if breaks is a single number, it will generate breaks based on range(x), not xlim.) Great! That's *exactly* what I needed. Thanks a lot! -- View this message in context:

[R] Emacs vs Eclipse vs Rcmdr

Hi everyone, I could have posted this on R-devel or the GUI list, however don't feel it's that serious. Hence, decided R-help would be the most appropriate. I'm not so much interested in which is the best user interface for R. Rather which is the best ***platform*** for developing ***new*** user

Re: [R] plot question

Hi: My understanding is that you want a separate plot over time of each of your repeated measures variables. If that's the case, then perhaps something along the following lines will work. Given the manufactured data set Gabor provided (g from lines), let's use the reshape package and then plot

[R] calculations on columns with partially matching names

Is there a command for partial matching of character strings? Specifically, I'd like to be able to calculate the mean of the values in any columns in a data frame or matrix that have identity in part of their column names. For example, columns labeled mpw06a and mpw06b match on the first five

[R] CRAN (and crantastic) updates this week

CRAN (and crantastic) updates this week New packages * dcemriS4 (0.20.1) Brandon Whitcher http://crantastic.org/packages/dcemriS4 A collection of routines and documentation that allows one to perform a quantitative analysis of dynamic contrast-enhanced or

Re: [R] R2HTML Report number format, or Better Way?

Thanks for the replys. I've played with both suggestions, but must not understand as I could not get them to work. In the first, do I have to save my results as a *.dat file? Where do I find the HTML file? For the second, As a beginner, is this what you mean by supply the data as code?

Re: [R] calculations on columns with partially matching names

On Jan 3, 2010, at 6:09 PM, Jim Bouldin wrote: Is there a command for partial matching of character strings? Specifically, I'd like to be able to calculate the mean of the values in any columns in a data frame or matrix that have identity in part of their column names. For example,

Re: [R] R2HTML Report number format, or Better Way?

On 01/04/2010 10:49 AM, L.A. wrote: Thanks for the replys. I've played with both suggestions, but must not understand as I could not get them to work. In the first, do I have to save my results as a *.dat file? No, I simply showed how I obtained the data frame that you gave as an

Re: [R] Emacs vs Eclipse vs Rcmdr

Hello On 1/3/10, Charlotte Maia mai...@gmail.com wrote: Rcmdr Pros and Cons: - I haven't used it for a long time, so can't really comment. - However, I was surprised by how many reverse dependencies it has. So I will assume it has some potential. Rcmdr is probably not the best editor

[R] Assistance with boot() Package

Good Evening R Community: I believe I understand the basics of using the boot() bootstrap resampling function in the boot() package. I have not had any trouble creating a boot.object to which I apply the boot.ci() function to calculate one or all of the available confidence intervals. What I

Re: [R] Questions bout SVM

2010/1/2 Nancy Adam nancyada...@hotmail.com: Hi Steve, Thanks a lot for your reply. 1)I’m still confused which equation (1- sqrt(mean(mymodel$MSE)) OR 2- mean(sqrt(mymodel$MSE)) )is equivalent to sqrt(mean(error**2))? So, as I mentioned before, mymodel$MSE is a vector that's as long as the

Re: [R] bio3d package not installing

On Jan 3, 2010, at 1:41 PM, Jonathan wrote: Is it just me, or is the bio3d package no longer available? Is there another way to do a reasonable nucleic acid sequence alignment in R? Sounds like a BioC question to me. Thanks, Jonathan David Winsemius, MD Heritage Laboratories West

Re: [R] bio3d package not installing

On Jan 3, 2010, at 2:11 PM, David Winsemius wrote: On Jan 3, 2010, at 1:41 PM, Jonathan wrote: Is it just me, or is the bio3d package no longer available? Is there another way to do a reasonable nucleic acid sequence alignment in R? Sounds like a BioC question to me. My (faulty)

Re: [R] Artificial Neural Networks

1) there are plenty of data for neural net testing in R. you might check datasets package on CRAN. 2) which neural net are you talking about, BP, RBF, LVQ, or something else. the world of neural nets is pretty much like a zoo. without knowing which animal you are talking about, nobody can help

[R] function in aggregate applied to specific columns only

I want to use aggregate with the mean function on specific columns gender - factor(c(m, m, f, f, m)) student - c(0001, 0002, 0003, 0003, 0001) score - c(50, 60, 70, 65, 60) basicSub - data.frame(student, gender, score) basicSubMean - aggregate(basicSub, by=list(basicSub$student), FUN=mean,

Re: [R] function in aggregate applied to specific columns only

On Jan 3, 2010, at 10:46 PM, david hilton shanabrook wrote: I want to use aggregate with the mean function on specific columns gender - factor(c(m, m, f, f, m)) student - c(0001, 0002, 0003, 0003, 0001) score - c(50, 60, 70, 65, 60) basicSub - data.frame(student, gender, score) basicSubMean -

Re: [R] function in aggregate applied to specific columns only

Hi: Perhaps the plyr package would be useful. It contains functions colwise(), numcolwise() and catcolwise() that will perform the same operation on the stated type of object. In this case, numcolwise() is appropriate: str(basicSub) 'data.frame': 5 obs. of 3 variables: $ student: num 1 2 3

Re: [R] function in aggregate applied to specific columns only

You want this? basicSubMean - aggregate(basicSub[c(score)], by=list(basicSub$student), FUN=mean, na.rm=TRUE) basicSubMean Group.1 score 1 1 55.0 2 2 60.0 3 3 67.5 bests milton On Sun, Jan 3, 2010 at 10:46 PM, david hilton shanabrook dhsha...@acad.umass.edu wrote: I

Re: [R] function in aggregate applied to specific columns only

Here are 6 ways: 1. aggregate aggregate(basicSub[score], basicSub[student], mean) student score 1 1 55.0 2 2 60.0 3 3 67.5 2. tapply with(basicSub, tapply(score, student, mean)) 123 55.0 60.0 67.5 3. summaryBy in doBy package library(doBy) summaryBy(. ~

Re: [R] function in aggregate applied to specific columns only

Just for the fun of it, here are two more: by and ave. with(basicSub, by(score, student, mean)) student: 1 [1] 55 student: 2 [1] 60 student: 3 [1] 67.5 Not my favorite

Re: [R] Likelihood Ratio Tests

Is there any package available in R to do the following hypothesis tests? Testing the means of two Poissons (equivalent to the difference of two proportions) Testing the equality of two proportions from binomials Testing the equality of proprtions of two negative binomials (both conditional and

[R] Extract vector elements until cumsum = x

Hi All, I have a vector n, and for each n[i] I want to extract n[i], n[i+1], n[i+2]..., until the cumulative sum of n[i] and subsequent elements exceeds a CheckValue, whereupon I move to the next index and repeat. I am trying to find a Vectorized approach, and have seen similar posts where

[R] Interpreting coefficient in selection and outcome Heckman models in sampleSelection

Hi there Within sampleSelection, I'm trying to calculate the marginal effects for variables that are present in both the selection and outcome models. For example, age might have a positive effect on probability of selection, but then a negative effect on the outcome variable. i.e.

[R] rggobi: Installation problems on OS X Snow Leopard

Hello: I have successfully installed both R and R64 for Apple OS X from CRAN plus GTK+2.14 and GGobi 2.1.8 from: http://www.ggobi.org/downloads I launch either R or R64 and try to install the rrgobi package and get (with R64) the following: R version 2.10.1 (2009-12-14) Copyright (C) 2009 The

[R] how to draw abline correctly?

Hello, I am frastruated with this graph, just cannot get what I need. Thank you for any suggestions or help. I really appreciate it. I wrote the following code, but there are 3 problems 1, the red line is added on the graph but without any marker on the y-axis. I want to display the number

Re: [R] R2HTML Report number format, or Better Way?

L.A. wrote: For the second, As a beginner, is this what you mean by supply the data as code? The actual data is 17000 rows and 31 columns. In your case, it is irrelevant if the data are 17000x31 or 2x3, so an example with few columns is sufficient. And if you supply the data in code,

Re: [R] xyplot: several plots in one creates y-scale problem

Hi As you did not provide an example as per posting guide here is one dataset cobbled from zoo help. It is fairly rough but is optioned up to give you an idea. The legend takes is cues from trellis.par.set() and you can set this for the plot by using par settings It makes use of strip.left

Re: [R] how to draw abline correctly?

Hi there, May be part of your problem is solved. size=c(80, 40, 20, 16, 10, 8, 4, 2) estimate=c(0.1353, 0.1337, 0.1237, 0.1164, 0.1058, 0.0983, 0.0759 , 0.0607) plot(size, estimate, type=o, lwd=1.5, col=blue, ylim=c(0, 0.2), xlim=c(0,80), pch=19, axes=F) abline(h=0.1361, col=red, lwd=1.5)