Re: [R] glm poisson function

ng 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 commented, minimal, self-contained, reproducible code. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e

Re: [R] comparing proportions

n/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.

Re: [R] comparing proportions

rform fairly well, as the Edgeworth series converges rapidly to normal for binomial distributions with p within 0.15 to 0.85 and 10+ replicates, as I stated before. I'd be interested in seeing the results of these 3 fits myself just for curiosity. At 01:21 PM 2/10/2011, array chip wro

Re: [R] comparing proportions

You need to change models 2 and 3 to use ~ group + subject. You left subject out as a fixed factor. At 05:17 PM 2/10/2011, array chip wrote: Robert, thank you! I tried all 3 models you suggested. Since each subject only has one line of data in my dataset, would including Subject as a factor

[R] expected behavior when parsing lines with special characters

285.2 6 32414251XP_327605.1 7 30683862NP_850592.1 8 15218215NP_175645.1 9 30699476NP_178116.2 10 115466360 NP_001056779.1 I would appreciate any feedback. Thanks, -Robert > sessionInfo() R version 2.12.1 (2010-12-16) Platform: x86_64-pc-mingw32/x64 (64-bit) locale: [1] LC_COLLA

Re: [R] [R-sig-ME] Fwd: Re: ANOVA and Pseudoreplication in R

org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. __ 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 an

Re: [R] problem with glm(family=binomial) when some levels have only 0 proportion values

___ 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 commented, minimal, self-contained, reproducible code. =

[R] defining a method for multiple classes

) object@x) setMethod("x", signature("myClass2"), function(object) object@x) I was wondering if it is possible to do the equivalent of: setMethod("x", signature(c("myClass1","myClass2")), function(object) object@x) Or am I always required

[R] Splitting Datasets

ing to learn fast. Thanks in advance, Robert [[alternative HTML version deleted]] __ 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.h

[R] saving contour() plot info

{ I think this message got rejected at the 1st attempt - trying again} R 2.15.1 , windows XP I have a very non-stationary bivariate time-series - say {xt,yt} t=1 ... lots. I want to do a bivariate density contour-plot of the whole series and then step through the series 1 second at a time

Re: [R] saving contour() plot info

Many thanks for these ideas ... I'll try them, and report back Cheers Bob Kinley -Original Message- From: Duncan Murdoch [mailto:murdoch.dun...@gmail.com] Sent: 03 July 2012 15:54 To: Robert Douglas Kinley Cc: r-help@r-project.org Subject: Re: [R] saving contour() plot

Re: [R] help with filled.contour() -

Take a look at the code for filled.contour(). You'll find a line beginning .Internal(filledcontour( You can adapt this line and the lines around it to achieve what you want. Good luck Bob Kinley -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-pro

Re: [R] saving contour() plot info

-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Robert Douglas Kinley Sent: 03 July 2012 16:17 To: Duncan Murdoch Cc: r-help@r-project.org Subject: Re: [R] saving contour() plot info Many thanks for these ideas ... I'll try them, and report back Cheers Bob K

Re: [R] differences between survival models between STATA and R

On 7/9/2012 9:17 AM, Javier Palacios Fenech wrote: Please. After, Terry's response I guess I was expecting to hear how your comparison between R and STATA went when you used the R function, survreg() for your analysis. We still don't know what your data look like. The posting guide asks for

Re: [R] Using MCMC sampling to estimate p values with a mixed model

, self-contained, reproducible code. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Bea

Re: [R] analysing a three level reponse

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 commented, minimal, self-contained, reproducible code. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS

Re: [R] coxph() - unexpected result using Crawley's seedlings data (The R Book)

uide.html and provide commented, minimal, self-contained, reproducible code. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Dri

[R] qpcR propagate problem

I am using the propagate function of the qpcR package to estimate the standard error for an expression. This expression is simple enough that I am able to calculate the first-order propagation of error, which is what the documentation on propagate states that it does. However, the results are not

Re: [R] AKIMA color contouring

__ 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 commented, minimal, self-contained, reproducible code.

[R] survival: fitting equation to survival curve?

Given a survfit object, is it possible to fit an equation to the resulting survival curve? What about with a coxph or survreg object? TIA, Rob [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mail

[R] survival: fitting equation to survival curve?

Given a survfit object, is it possible to fit an equation to the resulting survival curve? Is this possible? What about with a coxph or survreg object? TIA, Rob [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

Re: [R] survival: fitting equation to survival curve?

ities. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Lancaster, Robert (Orbitz) Sent: Wednesday, October 26, 2011 1:40 PM To: r-help@r-project.org Subject: [R] survival: fitting equation to survival curve? Given a survfit object, is it possib

Re: [R] How to plot survival data from multiple trials (simulations)?

I'm assuming you want to plot the Kaplan-Meier curves. First, load your data into a data.frame, e.g. surv_data It looks like ST is your survival time and ind is your event indicator. The following code should get you started: # create your survival object: mySurv = Surv(surv_data$ST,surv_data

Re: [R] How to plot survival data from multiple trials (simulations)?

Oh sorry, I misread. You'll only be able to get the confidence interval out of the object returned by survfit. From: Zheng Lu [mailto:luy...@hotmail.com] Sent: Saturday, October 29, 2011 1:04 PM To: Lancaster, Robert (Orbitz); r-help@r-project.org Subject: RE: [R] How to plot survival data

Re: [R] Kaplan Meier - not for dates

I think it really depends on what your event of interest is. If your event is that the patient got better and "left treatment" then I think this could work. You would have to mark as censored any patient still in treatment or any patient that stopped treatment w/o getting better (e.g. in the c

[R] [ reached getOption("max.print") -- omitted 8178 rows ]]

Hi, I have a weird thing I don¹t understand. Here¹s what I did: I read some data: data=read.table("fileName²) then I printed the data to the screen: data But it didn¹t finish: lot¹s of data was written out, but not all of it... Then it interrupted and said: [ reached getOption("max.print") -- omitt

Re: [R] Web based R-help not a list

On 11/9/11 7:22 PM, "David Winsemius" wrote: > > On Nov 9, 2011, at 5:38 PM, Cem Girit wrote: > >> Hello, >> >>Is there a web version of this R-Help user group >> (such as >> the ones under Google Groups) such that > > There is Nabble. It's not going to make you any friends

[R] subreddit for R related stuff

Thought there might be some Redditors lurking on these mailing lists. I created a sub-reddit for R (and by extension Bioconductor) discussions, links, etc. http://www.reddit.com/r/Rsoftware/ This will be the first and only shameless plug. -Robert Robert M. Flight, Ph.D. University of

Re: [R] Solving 100th order equation

Is this a homework problem? Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia

Re: [R] How to make R running faster

d minimize the number of lines of R. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-467-294

[R] Is there a good package for multiple imputation of missing values in R?

ut how to get a new data frame with the imputed values replaced (I don't have Herrell's book). Any pointers would be appreciated. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Form

Re: [R] Is there a good package for multiple imputation of missing values in R?

At 03:02 AM 6/30/2008, Robert A. LaBudde wrote: I'm looking for a package that has a start-of-the-art method of imputation of missing values in a data frame with both continuous and factor columns. I've found transcan() in 'Hmisc', which appears to be possibly suite

Re: [R] shapiro wilk normality test

ke you to the asymptotic region. Use sample() to select the data at random from within your data set to avoid bias in using the test. E.g., shapiro.test(sample(mydata, 1000, replace=TRUE)) ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e

Re: [R] shapiro wilk normality test

tribution. If you wish to test how the distribution is nonnormal, within some family of nonnormal distributions, you will have to specify such a null hypothesis and test for deviation from it. E.g., testing for coefficient of skewness = 0.

Re: [R] peak finding

short function might be useful in picking off peaks (and symmetrically modified, valleys): #findpeaks() #Copyright 2007 by Robert A LaBudde, all rights reserved #find peaks in (x,y) curve findpeaks<- function(x, y) { nx<- length(x) ny<- length(y) if (nx != ny) { print (&#x

Re: [R] generate random numbers subject to constraints

e random variates for X1, X2, based upon whatever unspecified distribution you wish. 2. Solve the two equations for X3 and X4. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.

Re: [R] generate random numbers subject to constraints

the interior >of a general (k-r)-dimensional simplex embedded in >k dimensions, with (k+1) given vertices? > The method of "rejection": 1. Generate numbers randomly in the hypercube. 2. Test to see if the point falls within the prescribed area. 3. Accept the point if it does. 4.

[R] Problem installing and using package "tseries"

;tseries' > There doesn't seem to be anything unusual. The package is installed like any other, and the required file is at C:\Program Files\R\R-2.6.2\library\tseries\libs\tseries.dll. It is 36,864 bytes and looks superficially okay. Any ideas as to what's wrong here? Thanks.

Re: [R] Problem installing and using package "tseries"

he problem and recompiled tseries which will appear on >CRAN master within 12 hours. > >Best wishes, >Uwe Ligges > > >Robert A. LaBudde wrote: >>I have R 2.6.2, and have tried downloading and installing the >>package "tseries". I get the

[R] Multidimensional contingency tables

7 black black yes 6 8 black blackno97 The xtabs() function doesn't seem appropriate, as it has no means of using 'count'. This must be a common problem. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail:

Re: [R] Multidimensional contingency tables

Now that is simple and elegant. Thanks! PS. Is there a course available for learning how to read R help information? :) At 10:52 PM 4/21/2008, Gabor Grothendieck wrote: > xtabs(count ~., prob1) > >On Mon, Apr 21, 2008 at 10:46 PM, Robert A. LaBudde <[EMAIL PROTECTED]> wrote: &

Re: [R] bootstrap for confidence intervals of the mean

https://stat.ethz.ch/mailman/listinfo/r-help >PLEASE do read the posting guide http://www.R-project.org/posting-guide.html >and provide commented, minimal, self-contained, reproducible code. Robert A. LaBudde,

Re: [R] quasi-random sequences

ex.ac.uk/atto > >__ >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 commented, minimal, self-contained, reproducible code. ===

Re: [R] Simulation study in R

to know how to sample X in the first place. Is its distribution uniform, or something else? ======== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Ti

Re: [R] Averaging 'blocks' of data

PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Fo

Re: [R] t-test between percentages

. __ 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 commented, minimal, self-contained, reproducible code. Robert

[R] R function which finds confidence interval for binomial variance

#x27;? Thanks. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-467-2947 "Vere

Re: [R] R function which finds confidence interval for binomial variance

ieve reasonable coverage. So I'm still looking for a reliable method for all p and for reasonable n. The proportion-based method is the best I've found, so far. ======== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTE

[R] Confidence interval for binomial variance

t. #09.26.08 02.50 binomVarCI.r #copyright 2008 by Robert A LaBudde, all rights reserved #CI for binomial sample variance #created: 09.26.08 by r.a. labudde #changes: require('binGroup') binomVarCI<- function (n, x, conf=0.95) { p<- x/n #proportion if (p<0.25 | p>0.75

Re: [R] FW: logistic regression

not take into account the stepwise >>>> procedure above that was used to torture the data until they >>>> confessed. >>>> >>>> Frank >>> Please book this as a fortune. >>> >>> Dieter >> Seconded! >> T

Re: [R] Located Latent Class Analysis (Uebersax)

ting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. ======== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.

Re: [R] Logistic regression problem

Tura, M.D,MPH,Ph.D National Institute of Cardiology Brazil __ 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 commente

[R] termplot

I have used glm and stepAIC to choose a best model. I can use termplot to assess the contribution of each explanatory variable in the glm. However the final model after running stepAIC includes interaction terms, and when I do termplot I get "Error in `[.data.frame`(mf, , i) : undefined columns s

[R] glm StepAIC with all interactions and update to remove a term vs. glm specifying all but a few terms and stepAIC

Problem: I am sorting through model selection process for first time and want to make sure that I have used glm, stepAIC, and update correctly. Something is strange because I get a different result between: 1) a glm of 12 predictor variables followed by a stepAIC where all interactions are consi

[R] Unscaled factanal()

most instances, but occasionally I'd like to perform the analysis directly upon the covariance matrix itself. Thanks. ======== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.

[R] How to do multi-factor stratified sampling in R

ction marginal proportions as well for some pairs. This must be a common problem in data mining, but I don't seem to be able to locate the proper library or function for doing this in R. Thanks for any help. ======== Robert A. LaBudde, PhD

Re: [R] How to do multi-factor stratified sampling in R

At 03:54 PM 3/8/2008, David Winsemius wrote: >"Robert A. LaBudde" <[EMAIL PROTECTED]> wrote in >news:[EMAIL PROTECTED]: > > > Given a set of data with a number of variables plus a response, I'd > > like to obtain a randomized subset of the rows such

Re: [R] stepAIC and polynomial terms

ge to see what the dependence is. Frequently curve is flat up to a certain age, and then linear thereafter. This gives rise to a pseudo-quadratic relationship. You should be able to fit it better with the split plus a linear term. 5. Think about how age should affect your response before trying model

Re: [R] removing a specific number of digist from a character string

test <- c("060907_17_3_5_1_1_2909.tif", "060907_17_3_5_2_1_2910.tif", "060907_17_3_5_3_1_2911.tif") sub('[[:digit:]][[:digit:]][[:digit:]][[:digit:]]\.tif', '', test) or test <- c("060907_17_3_5_1_1_2909.tif", "060907_17_3_5_2_1_2910.tif", "060907_17_3_5_3_1_2911.tif") sub('[[:digit:]]{4}\.ti

Re: [R] data frame

I believe you're looking for: dim(a) dim(a)[1] # Number of observations, in your example, 12 dim(a)[2] # Number of variables per observation, in your example, 9 --Jeff. On Sep 17, 2007, at 12:05 PM, Alfredo Alessandrini wrote: > Hi everybody, > > If I've a data frame like this: > > datafram

Re: [R] removing a specific number of digist from a character string

For the sake of absolute correctness: > sub('[[:digit:]]{4}\.tif', '', test) should be sub('[[:digit:]]{4}\\.tif', '', test) -- Jeff. On Sep 17, 2007, at 11:59 AM, Jeffrey Robert Spies wrote: > test <- c("060907_17_3_5_1_1_2909.tif&quo

Re: [R] Need help on "date"

And one using regular expressions: x <- "2005-09-01" pattern <- '([[:digit:]]{4})-([[:digit:]]{2})-([[:digit:]]{2})' y <- sub(pattern, '\\1', x) m <- sub(pattern, '\\2', x) d <- sub(pattern, '\\3', x) -- Jeff. On Sep 18, 2007, at 5:00 AM, Arun Kumar Saha wrote: > Dear all, > > I have a variable

Re: [R] Need help on "date"

, this site is quite useful: http://www.cs.utah.edu/dept/old/texinfo/regex/regex_toc.html. Make sense? Jeff. On Sep 18, 2007, at 10:46 AM, Arun Kumar Saha wrote: > Dear Jeffrey, > > Your syntax looks very extraordinary to me. I would be very happy > if you can explain this notation. &g

Re: [R] A matrix multiplication

How about this: a <- matrix(cbind(rep(2, 500), rep(3, 500)), 500, 2) b <- matrix(cbind(rep(5, 500), rep(6, 500), rep(7, 500)), 500, 3) matrix(apply(a, c(2), "*", b), nrow=500, ncol=6) We apply the multiplier (quoted as specified in the apply help) with argument b to every column of a as specif

Re: [R] simple matching with R

I hope this will change by and by. The more examples you see and play with, the more you'll understand. > So I would be very pleased if you could help me once again. > > Greetings > > Birgit Cheers, Jeff. > > Am 28.09.2007 um 18:25 schrieb Jeffrey Robert Spies: &

Re: [R] simple matching with R

Not sure how you want to handle the NAs, but you could try the following: #start MalVar29_37 <- read.table(textConnection("V1 V2 V3 V4 V5 V6 V7 V8 V9 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 NA NA NA NA NA NA NA NA NA 0 1 0 0 0 1 0 0 0"), header=TRUE)

Re: [R] convert time series to data.frame

Hi Edna, Can you send a small subset of the data as an example and the function call you used to read the data in originally? It might be helpful in understanding why you're losing the "time element". Jeff. On Oct 1, 2007, at 12:27 AM, Edna Bell wrote: > Dear R gurus > > I would like to ta

Re: [R] mean of subset of rows

You were on the right track with the for loop, but often you can do the same thing looplessly (I know, it's not really a word) in R: If your data is like this: data<-data.frame(ID=rep(letters[1:4], 5), size=runif(20)) then apply either tapply(data$size, data$ID, mean) or aggregate(data$size

Re: [R] standard errors of regression coefficients

stderr_int <- summary(lm(y ~ x))$coefficients[1,2] stderr_slope <- summary(lm(y ~ x))$coefficients[2,2] Jeff. On Oct 3, 2007, at 3:01 AM, Alexander Moreno wrote: > Hi, > > If I have two vectors x and y and I do lm(y~x) and now I want to > define > variables that are the standard errors of the

[R] [R} Getting 'tilting' confidence intervals in R

[Resend with proper subject line] I am trying to compute bootstrap confidence intervals for a sample using R 2.5.1 for Windows. I can get "Normal", "Basic", "Percentile", "BCa" and "ABC" from boot.ci() and boot() in the Davison & Hinkley "boot" package. But I can't figure out how to use tilt.b

[R] Getting 'tilting' confidence intervals in R

n the tilt.boot() call. This must be simple, but I'm not finding out the answer from the documentation. Also, use of any other packages to accomplish the tilting interval would also be useful. Thanks. Robert A. LaBudde, PhD,

Re: [R] Inferences for ratios of non-normal means

t;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 commented, minimal, self-contained, reproducible code.

Re: [R] applying math/stat functions to rows in data frame

rest? Or, must these summary functions be programmed >separately to work on a row? Try using t() to transpose the matrix, and then apply the column function of interest. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [

[R] Formatting multi-way ANOVA output for spectra analysis

I've studied R a little bit, although I haven't used it in some time (except via RCommander).  I'm working on my dissertation project and have spectrometer data that I need to evaluate.  I need to find a way to simplify the output from multi-way ANOVA so I can reduce the areas of the spectrum t

[R] Segmentation Fault on Unix box with nloptr, works on Windows

I have an R script that is giving me a Segmentation Fault depending on the size of the dataset. It is only happening on our Unix installation of R Server. I am able to run it against a Windows server with the exact same data and script successfully. The Segmentation Fault occurs when I call n

[R] help with line graphs - rather lengthy to explain need

I am trying to figure out the best way to organize and plot data generated by a Excel spreadsheet (one driving a sample turntable and collecting optical spectra). The output of the equipment and software is an excel spreadsheet with sample numbers in the first row, and in the first column ther

Re: [R] Dealing with NaN's in data frames

Hi Jon, Here's one way. > x <- c(1,2,3,4,NaN) > y <- c(1,2,NaN,4,5) > > myDF <- data.frame(x,y) > myDF x y 1 1 1 2 2 2 3 3 NaN 4 4 4 5 NaN 5 > > myDF[ is.na(myDF) ] <- NA > myDF x y 1 1 1 2 2 2 3 3 NA 4 4 4 5 NA 5 Cheers, Bob =

Re: [R] teaching R

Hello all, I am going to be running a small statistics workshop using R sometime in November. I am restricted to R because of the specific libraries I will be using - a good thing in my book - however the attendees are unfamiliar with R. I plan on giving as little R information as possible - ju

Re: [R] SPSS Statistics-R Integration Plug-In

Hi Michael, I've used the R plug-in and really like it. You can read my instructions on how to install and use it by going to Amazon.com, searching for the book, "R for SAS and SPSS Users" and then "search inside the book" for the section, "Running R from SPSS". I've only got about 3 pages on it (

Re: [R] How good is R at making publication quality tables?

Hi Paul, Sorry I didn't get to that subject in the first edition of R for SAS and SPSS Users. Several of the options people have mentioned will be in the second edition, although that's about a year off. I did get them added to R for Stata Users, due out in early April. Cheers, Bob >-Ori

[R] Windowing issue with diagram package & R 9.1

,1]<-1 ;M[4,2]<-2;M[3,4]<-3;M[1,3]<-4 Col <- M Col[] <- "black" Col[4,2] <- "darkred" pp<-plotmat(M,pos=c(1,2,1),curve=0.2,name=names,lwd=1,box.lwd=2,cex.txt=0.8, arr.type="triangle",box.size=0.1,box.type="hexa",box.prop=0.25, arr.co

[R] Reading data entered within an R program

Dear R-helpers, I know of two ways to reading data within an R program, using textConnection and stdin (demo program below). I've Googled about and looked in several books for comparisons of the two approaches but haven't found anything. Are there any particular advantages or disadvantages to thes

Re: [R] Reading data entered within an R program

Since stdin seemed simpler I figured textConnection must have some advantage. Thanks! Bob From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com] Sent: Saturday, July 11, 2009 6:00 PM To: Muenchen, Robert A (Bob) Cc: R-help@r-project.org Subject: Re: [R] Reading data entered within an

Re: [R] Rcmdr GUI goes into loop via alt+backspace under Windows OS

It just produces the "bell" sound on my 32-bit windows XP machine runing R 2.9.1 Is this really a "standard" combination? I know a lot of programs that use ctrl-Z, but I've never come across this shortcut key combination for undo. - Original Message - From: "tradenet" To: Sent:

[R] different L2 regularization behavior between lrm, glmnet, and penalized?

packages has insight into why the results differ, I'd appreciate hearing about it. I'm new to R. Thanks. (Version info below. Same results on Windows and Solaris 8, except that I haven't gotten glmnet to compile on the latter.) Robert V (Bob) Sasseen [1]sass...@ai.sri.co

[R] Frequencies, proportions & cumulative proportions

Dear R-Helpers, I've looked high and low for a function that provides frequencies, proportions and cumulative proportions side-by-side. Below is the table I need. Is there a function that already does it? Thanks, Bob > # Generate some test scores > myValues <- c(70:95) > Score <- ( sample( myVal

Re: [R] Frequencies, proportions & cumulative proportions

Ted, I know how to do that. It's just such a standard display in SAS, SPSS and Stata that I figured someone had done it and I had just overlooked it. Thanks! Bob I don't think there is a ready-made one, but it is very little effort to make your own: mkMyTable <- function(X){ Table <- data.f

Re: [R] Frequencies, proportions & cumulative proportions

0.043 | |---| ... -Original Message- From: David Scott [mailto:d.sc...@auckland.ac.nz] Sent: Friday, October 16, 2009 8:42 PM To: Muenchen, Robert A (Bob) Cc: ted.hard...@manchester.ac.uk; r-help@r-project.org Subject: Re: [R] Frequencies, proportions & cumulative proportions Muenchen, Robert A (Bob)

[R] R Packages Crack the 3,000 Mark!

Hi All, I don't know if this has been reported before, but according to Henrique Dallazuanna's program (below) the number of R packages has exceeded the 3,000 mark. The count today is 3,175. I ran this just a couple of months ago & the number was still in the high 2,000s, so it must be fairly r

Re: [R] R Packages Crack the 3,000 Mark!

: Muenchen, Robert A (Bob) Cc: r-help@r-project.org Subject: Re: [R] R Packages Crack the 3,000 Mark! Hello On 11/24/09, Muenchen, Robert A (Bob) wrote: > I don't know if this has been reported before, but according to Henrique > Dallazuanna's program (below) the number of R packages

Re: [R] R Packages Crack the 3,000 Mark!

I thought that the unique function would eliminate duplicate package names. Is there a better way to count the number of packages? Thanks, Bob -Original Message- From: Gabor Grothendieck [mailto:ggrothendi...@gmail.com] Sent: Wednesday, November 25, 2009 10:40 AM To: Muenchen, Robert A

Re: [R] different L2 regularization behavior between lrm, glmnet, and penalized?

The author of the penalized package, j.j.goeman at lumc.nl, kindly replied to my message. He also responded to another question I asked him. -- The differences have to do with different scaling defaults. lrm by default standardizes the covariates to unit sd before applying pen

Re: [R] renaming column names

First, data is an R function so it is better to avoid it as a name for another object. For example, use dat instead. Try this: data(iris) dat=iris colnames(iris) [1] "Sepal.Length" "Sepal.Width" "Petal.Length" "Petal.Width" "Species" colnames(dat) [1] "Sepal.Length" "Sepal.Width" "Petal.Le

[R] Updated comparison table for SAS-SPSS Add-ons and R Functions

Hi All, I have substantially expanded the table that compares SAS and SPSS add-on modules to somewhat equivalent R packages. This new version is at: http://r4stats.com/add-on-modules and I would very much appreciate any feedback you might have on it. The site http://r4stats.com is the replacemen

Re: [R] Updated comparison table for SAS-SPSS Add-ons and R Functions

>From: b.rowling...@googlemail.com [mailto:b.rowling...@googlemail.com] On >Behalf Of Barry Rowlingson >Sent: Wednesday, January 13, 2010 7:03 PM >To: Muenchen, Robert A (Bob) >Cc: r-help@r-project.org >Subject: Re: [R] Updated comparison table for SAS-SPSS Add-ons and R Func

Re: [R] Updated comparison table for SAS-SPSS Add-ons and R Functions

Hi Liviu, Thanks for those suggestions. I've made the changes and added you to the list of contributors. Cheers, Bob > -Original Message- > From: Liviu Andronic [mailto:landronim...@gmail.com] > Sent: Thursday, January 14, 2010 7:06 AM > To: Muenchen, Robert

Re: [R] Would Like Some Advice

Besides monitization, Windows has a few other things that infuriate me... (1) VERY hard to control updates, (2) "sneaker" updates - things installed that people don't want (like trying to force Windows computer owners to update - and sometimes wrecking the computer when it does), (3) bad update

[R] R Online Workshops October 7-11

Learn R and/or data mangement at home October 7 through 11 http://r4stats.com/2013/09/11/learn-r-andor-data-management-from-home-october-7-11/ == Bob Muenchen (pronounced Min'-chen) Accredited Professional Statistician(tm) Manager, Researc

[R] How Rcmdr or na.exclude blocks TukeyHSD

Dear R-Helpers, I was calling the TukeyHSD function and not getting confidence intervals or p-values. It turns out this was caused by missing data and the fact that I had previously turned on R Commander (Rcmdr). John Fox knew that Rcmdr sets na.action to na.exclude, which causes the problem. I

[R] Popularity of R, SAS, SPSS, Stata...

Hi All, I've been fiddling around with various ways to estimate the popularity of R, SAS, SPSS, Stata, JMP, Minitab, Statistica, Systat, BMDP, S-PLUS, R-PLUS and Revolution R. It's not an easy task. You can see what I've come up with so far at http://r4stats.com/popularity . I'm sure people will h

<    3   4   5   6   7   8   9   >