[R] summary of linear fixed effects model is different than the HSAUR book

2007-07-22 Thread Christopher W. Ryan
Running R 2.5.1 and a newly downloaded lme4 package on WinXP I'm trying to work my way through Everitt and Hothorn's Handbook of Statistical Analyses Using R, c 2006. (No, it's not homework.) Chapter 10 discusses linear mixed effects models for longitudinal data. I've called my long data frame

Re: [R] summary of linear fixed effects model is different than the HSAUR book

2007-07-22 Thread Peter Dalgaard
Christopher W. Ryan wrote: But on page 169, summary() is shown to produce additional columns in the fixed effects section, namely degrees of freedom and the P-value (with significance stars). How can I produce that output? Am I doing something wrong? Has lme4 changed? The latter. To

[R] summary statistics for groups

2007-07-16 Thread Kirsten Beyer
I have a list of ICD9(disease) codes and associated exposure measures. I am interested in a script that would calculate the mean exposure for each ICD9 code and the combined mean for all other codes excluding that code, then attaching those measures, and a count of how many numbers were included

[R] summary and min max

2007-04-23 Thread Sebastian P. Luque
Hi, I came across a case where there's a discrepancy between minimum and maximum values reported by 'summary' and the 'min' and 'max' functions: --cut here---start- R str(tt) num [1:1397] 1952 1970 1976 1967 1946 ... R summary(tt) Min. 1st Qu

Re: [R] summary and min max

2007-04-23 Thread Ben Bolker
Sebastian P. Luque spluque at gmail.com writes: Hi, I came across a case where there's a discrepancy between minimum and maximum values reported by 'summary' and the 'min' and 'max' functions: By default summary only lists 3 significant digits ... see ?summary Ben Bolker (is

Re: [R] summary and min max

2007-04-23 Thread Mike Prager
Sebastian P. Luque [EMAIL PROTECTED] wrote: I came across a case where there's a discrepancy between minimum and maximum values reported by 'summary' and the 'min' and 'max' functions: summary() rounds by default. Thus its reporting oddball values is considered a feature, not a bug. -- Mike

Re: [R] summary and min max

2007-04-23 Thread Robert Duval
Has anyone created an alternative summary method where the rounding is made only for digits to right of the decimal point? I personally don't like the way summarize works on this particular issue, but I'm not sure how to modify it generically... (of course one can always set digits=something_big

Re: [R] summary and min max

2007-04-23 Thread Bert Gunter
PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Robert Duval Sent: Monday, April 23, 2007 4:16 PM To: r-help@stat.math.ethz.ch Subject: Re: [R] summary and min max Has anyone created an alternative summary method where the rounding is made only for digits to right of the decimal point? I personally

Re: [R] summary polr

2007-02-20 Thread Michael Dewey
At 15:21 19/02/2007, Paolo Accadia wrote: Hi all, I have a problem to estimate Std. Error and t-value by “polr” in library Mass. They result from the summary of a polr object. I can obtain them working in the R environment with the following statements: temp - polr(formula =

Re: [R] summary polr

2007-02-20 Thread Michael Dewey
At 14:41 20/02/2007, you wrote: Please do not just reply to me, 1 - I might not know 2 - it breaks the threading Hi here there is an example extracted from polr help in MASS: The function could be: temp - function(form, dat) { house.plr - polr(formula = form,

Re: [R] summary polr

2007-02-20 Thread Prof Brian Ripley
That's not it (the function is 'coef' not 'coeff', and R can tell functions and lists apart). If you read the help page for polr you will see you could have used Hess=TRUE. It works then. THAT is why we needed an example, to see how you used the function. On Tue, 20 Feb 2007, Michael

Re: [R] summary polr

2007-02-20 Thread Paolo Accadia
Thank you very much Paolo Prof Brian Ripley [EMAIL PROTECTED] 20/02/07 4:00 PM That's not it (the function is 'coef' not 'coeff', and R can tell functions and lists apart). If you read the help page for polr you will see you could have used Hess=TRUE. It works then. THAT is why we needed

[R] summary polr

2007-02-19 Thread Paolo Accadia
Hi all, I have a problem to estimate Std. Error and t-value by “polr” in library Mass. They result from the summary of a polr object. I can obtain them working in the R environment with the following statements: temp - polr(formula = formula1, data = data1) coeff - summary(temp),

Re: [R] summary of the effects after logistic regression model

2007-01-26 Thread Frank E Harrell Jr
andrea evangelista wrote: Dear all, my aim is to estimate the efficacy over time of a treatment for headache prevention. Data consist of long sequences of repeated binary outcomes (1 if the subject has at least 1 episode of headache , 0 otherwise) on subjects randomized to placebo or

[R] summary of the effects after logistic regression model

2007-01-25 Thread andrea evangelista
Dear all, my aim is to estimate the efficacy over time of a treatment for headache prevention. Data consist of long sequences of repeated binary outcomes (1 if the subject has at least 1 episode of headache , 0 otherwise) on subjects randomized to placebo or treatment. I have fit a logistic

[R] Summary Tables for t.test?

2006-12-19 Thread Robert Lischke
Hi there, kind of a newbie-question, but anyway: is there a way to aggregate a nice table with input and results of a t-Test? I'm looking for something like a summary of Mean, Median, SD (seperated for both distributions) and p-value, df, etc. for the t-Test. I'm asking because i need a dozen

Re: [R] Summary shows wrong maximum

2006-12-07 Thread Mike Prager
Bert Gunter Nonclinical Statistics 7-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Mike Prager Sent: Wednesday, December 06, 2006 11:46 AM To: r-help@stat.math.ethz.ch Subject: Re: [R] Summary shows wrong maximum I don't know about

Re: [R] Summary shows wrong maximum

2006-12-06 Thread Prof Brian Ripley
'Unfortunately' you give no credentials for your ex cathedra pronouncement. E.g. http://en.wikipedia.org/wiki/Significant_digits says The situation regarding trailing zero digits that fall to the left of the decimal place in a number with no digits provided that fall to the right of the

Re: [R] Summary shows wrong maximum

2006-12-06 Thread rolf
Brian Ripley wrote: 'Unfortunately' you give no credentials for your ex cathedra pronouncement. E.g. http://en.wikipedia.org/wiki/Significant_digits says The situation regarding trailing zero digits that fall to the left of the decimal place in a number with no digits provided that

Re: [R] Summary shows wrong maximum

2006-12-06 Thread Bert Gunter
Folks: Is So this is at best a matter of opinion, and credentials do matter for opinions. -- Brian Ripley an R fortunes candidate? -- Bert Gunter On Tue, 5 Dec 2006, Oliver Czoske wrote: On Mon, 4 Dec 2006, Uwe Ligges wrote: Sebastian Spaeth wrote: Hi all, I have a list with a

Re: [R] Summary shows wrong maximum

2006-12-06 Thread Bert Gunter
... -- Bert Bert Gunter Nonclinical Statistics 7-7374 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Mike Prager Sent: Wednesday, December 06, 2006 11:46 AM To: r-help@stat.math.ethz.ch Subject: Re: [R] Summary shows wrong maximum I don't know about candidacy

Re: [R] Summary shows wrong maximum

2006-12-05 Thread Oliver Czoske
On Mon, 4 Dec 2006, Uwe Ligges wrote: Sebastian Spaeth wrote: Hi all, I have a list with a numerical column cum_hardreuses. By coincidence I discovered this: max(libs[,cum_hardreuses]) [1] 1793 summary(libs[,cum_hardreuses]) Min. 1st Qu. MedianMean 3rd Qu.Max.

[R] Summary shows wrong maximum

2006-12-04 Thread Sebastian Spaeth
Hi all, I have a list with a numerical column cum_hardreuses. By coincidence I discovered this: max(libs[,cum_hardreuses]) [1] 1793 summary(libs[,cum_hardreuses]) Min. 1st Qu. MedianMean 3rd Qu.Max. 1 2 4 36 141790 (note the max value of 1790)

Re: [R] Summary shows wrong maximum

2006-12-04 Thread Uwe Ligges
Sebastian Spaeth wrote: Hi all, I have a list with a numerical column cum_hardreuses. By coincidence I discovered this: max(libs[,cum_hardreuses]) [1] 1793 summary(libs[,cum_hardreuses]) Min. 1st Qu. MedianMean 3rd Qu.Max. 1 2 4 36 14

Re: [R] Summary shows wrong maximum

2006-12-04 Thread Gavin Simpson
On Mon, 2006-12-04 at 12:04 +0100, Sebastian Spaeth wrote: Hi all, I have a list with a numerical column cum_hardreuses. By coincidence I discovered this: max(libs[,cum_hardreuses]) [1] 1793 summary(libs[,cum_hardreuses]) Min. 1st Qu. MedianMean 3rd Qu.Max. 1

Re: [R] Summary shows wrong maximum

2006-12-04 Thread Sebastian Spaeth
Uwe Ligges wrote: max(libs[,cum_hardreuses]) [1] 1793 summary(libs[,cum_hardreuses]) Min. 1st Qu. MedianMean 3rd Qu.Max. 1 2 4 36 141790 (note the max value of 1790) Ouch this is bad! Anything I can do to remedy this? Known bug? No, it's a

[R] Summary, was Re: Confidence interval for relative risk

2006-11-22 Thread Michael Dewey
At 14:43 10/11/2006, Michael Dewey wrote: After considerable help from list members and some digging of my own I have prepared a summary of the findings which I have posted (see link below). Broadly there were four suggestions 1 - Wald-type intervals, 2 - transforming the odds ratio confidence

[R] summary linear regression

2006-11-09 Thread Gorka Merino
Good morning, I am a PhD student in Barcelona working with R. My question is about the summary of linear regressions. The output of that summary gives some statistical parameters of the regression. One of them is the R-squared. In the help menu i have read that the manner to calculate the

Re: [R] summary linear regression

2006-11-09 Thread Peter Alspach
Subject: [R] summary linear regression Good morning, I am a PhD student in Barcelona working with R. My question is about the summary of linear regressions. The output of that summary gives some statistical parameters of the regression. One of them is the R-squared. In the help menu i have

[R] SUMMARY control L to clear the Rgui screen in Windows

2006-09-28 Thread Charles Annis, P.E.
Eternal thanks to Jim Holtman and to Gabor Grothendieck who pointed me to this concise piece of code: cls - function() { require(rcom) wsh - comCreateObject(Wscript.Shell) comInvoke(wsh, SendKeys, \014) invisible(wsh) } Perfect! Thanks! Charles Annis, P.E.

Re: [R] summary(glm) for categorical variables

2006-09-12 Thread John Vokey
You want: anova(my.glm) On 12-Sep-06, at 4:00 AM, [EMAIL PROTECTED] wrote: Suppose we have a data.frame where variables are categorical and the response is categorical eg: my.df=NULL for(i in LETTERS[1:3]){my.df[[i]]=sample(letters, size=10)} my.df=data.frame(my.df)

[R] summary(glm) for categorical variables

2006-09-11 Thread Eleni Rapsomaniki
Dear list people Suppose we have a data.frame where variables are categorical and the response is categorical eg: my.df=NULL for(i in LETTERS[1:3]){my.df[[i]]=sample(letters, size=10)} my.df=data.frame(my.df) my.df$class=factor(rep(c(pos, neg), times=5)) my.glm=glm(class ~ ., data=my.df,

[R] Summary --- Local library under Windoze.

2006-09-01 Thread Rolf Turner
Prof. Brian Ripley solved the problem. He wrote: I was not aware that this works with relative paths for any version of R. Try using a full path, which always works for me. I tried it using a full path, and bingo! It worked like a charm. Under Unix the relative path

Re: [R] Summary --- Local library under Windoze.

2006-09-01 Thread Prof Brian Ripley
install.packages() in R 2.4.0 will make a full path out of a relative path, to avoid any confusion. (normalizePath is a good way to do that, BTW). On Fri, 1 Sep 2006, Rolf Turner wrote: Prof. Brian Ripley solved the problem. He wrote: I was not aware that this works with relative paths

[R] Summary and thanks: .Rprofile under Windoze.

2006-08-31 Thread Rolf Turner
Many thanks to all who responded to my cri de coeur: Charles Annis, John Bollinger, Petr Pikal, Brian Ripley, Duncan Murdoch, Dan Nordlund, and Richard Heiberger. (I don't think --- I hope --- I haven't missed anyone.) The essential solution to my problem was indeed to eliminate that @[EMAIL

[R] summary(lm ... conrasts=...)

2006-08-22 Thread Ted Harding
Hi Folks, I've encountered something I hadn't been consciously aware of previously, and I'm wondering what the explanation might be. In (on another list) using R to demonstrate the difference between different contrasts in 'lm' I set up an example where Y is sampled from three different normal

Re: [R] summary(lm ... conrasts=...)

2006-08-22 Thread Prof Brian Ripley
On Tue, 22 Aug 2006, [EMAIL PROTECTED] wrote: Hi Folks, I've encountered something I hadn't been consciously aware of previously, and I'm wondering what the explanation might be. Try contr.helmert(letters[1:3]) [,1] [,2] a -1 -1 b1 -1 c02

[R] summary statistics on an entire data frame

2006-08-10 Thread Patrick Aber
Hello, I am wondering how to compute functions like mean(), sd(), etc. on a data frame, but instead of getting a vector with the summary stat calculated individually for each column, get one number for the whole data set. I have noticed that coercing the frame to a time series (ts) gives the

Re: [R] summary statistics on an entire data frame

2006-08-10 Thread Petr Pikal
: Thu, 10 Aug 2006 11:10:35 -0400 Subject:[R] summary statistics on an entire data frame Hello, I am wondering how to compute functions like mean(), sd(), etc. on a data frame, but instead of getting a vector with the summary stat calculated individually for each column, get one

Re: [R] Summary method needed?

2006-08-03 Thread Heinz Tuechler
At 21:04 02.08.2006 +0100, Prof Brian Ripley wrote: On Wed, 2 Aug 2006, Christian Hennig wrote: Thank you Brian! I'm updating my fpc package at the moment and will add some new functions. I learned that there should be print and summary methods for the key functions. for

Re: [R] Summary method needed?

2006-08-03 Thread Martin Maechler
HeinzT == Heinz Tuechler [EMAIL PROTECTED] on Thu, 03 Aug 2006 09:39:35 +0100 writes: HeinzT At 21:04 02.08.2006 +0100, Prof Brian Ripley wrote: On Wed, 2 Aug 2006, Christian Hennig wrote: Thank you Brian! I'm updating my fpc package at the moment and will

[R] Summary method needed?

2006-08-02 Thread Christian Hennig
Hi list, I'm updating my fpc package at the moment and will add some new functions. I learned that there should be print and summary methods for the key functions. The purpose of the summary methods seems to be to reduce the possibly incredibly complex information in the function's output and

Re: [R] Summary method needed?

2006-08-02 Thread Prof Brian Ripley
On Wed, 2 Aug 2006, Christian Hennig wrote: Hi list, I'm updating my fpc package at the moment and will add some new functions. I learned that there should be print and summary methods for the key functions. for 'classes', I think. The purpose of the summary methods seems to be to

Re: [R] Summary method needed?

2006-08-02 Thread Christian Hennig
Thank you Brian! I'm updating my fpc package at the moment and will add some new functions. I learned that there should be print and summary methods for the key functions. for 'classes', I think. Yes. But in some cases the print method will make use of more or less all the output

Re: [R] Summary method needed?

2006-08-02 Thread Prof Brian Ripley
On Wed, 2 Aug 2006, Christian Hennig wrote: Thank you Brian! I'm updating my fpc package at the moment and will add some new functions. I learned that there should be print and summary methods for the key functions. for 'classes', I think. Yes. But in some cases the print

[R] SUMMARY: R crash with ''library(Matrix); as(x, dgCMatrix)' [was: Warning while subsetting...]

2006-07-08 Thread Thaden, John J
With thanks to Matrix package co-author Martin Maechler, I'm happy to report satisfactory closure of two recent threads I initiated about that package: - Warning while subsetting with Matrix - R Crash with 'library(Matrix);as(x,dgCMatrix) In the first, I reported seeing a warning

[R] Summary Statistics for data.frame

2006-07-08 Thread justin rapp
I apologize for my constant questions but I am new to R and trying to gain an appreciation for its capabilities. The following task is easy in Excel and I was hoping somebody could give me a quick explanation for how it can be acheived in R so I can avoid having to switch between the two

Re: [R] Summary Statistics for data.frame

2006-07-08 Thread Duncan Murdoch
On 7/8/2006 3:44 PM, justin rapp wrote: I apologize for my constant questions but I am new to R and trying to gain an appreciation for its capabilities. The following task is easy in Excel and I was hoping somebody could give me a quick explanation for how it can be acheived in R so I can

Re: [R] Summary Statistics for data.frame

2006-07-08 Thread justin rapp
When I attempt to use the mysummary function, I obtain the following error: Error in var(x) : missing observations in cov/cor When I use: by(data.logistic,data.logistic$Ydrafted,summary) I receive no errors. I cut and pasted your mysummary function directly into my r console. Should I have

Re: [R] Summary Statistics for data.frame

2006-07-08 Thread Duncan Murdoch
On 7/8/2006 4:55 PM, justin rapp wrote: When I attempt to use the mysummary function, I obtain the following error: Error in var(x) : missing observations in cov/cor var() gives that error if it sees NA values. You can get it to remove them by using var(x, na.rm = TRUE) instead of var(x).

[R] SUMMARY: making contour plots using (x,y,z) data

2006-07-01 Thread Ajay Narottam Shah
Folks, A few days ago, I had asked a question on this mailing list about making a contour plot where a function z(x,y) is evaluated on a grid of (x,y) points, and the data structure at hand is a simple table of (x,y,z) points. As usual, R has wonderful resources (and subtle complexity) in doing

[R] summary(lm(x~y)) difference between R-2.2.1 and R-2.3.0

2006-04-27 Thread Robin Hankin
Hi [macOSX 10.4.6; R-2.3.0] I have encountered a difference in behaviour between R-2.2.1 and R-2.3.0 when performing a linear model. Transcript follows for R-2.3.0 (R-2.2.1 worked as expected). How to make R-2.3.0 perform as R-2.2.1 did? dput(x) c(29.13, 29.88, 30.09, 29.99, 29.74,

Re: [R] summary(lm(x~y)) difference between R-2.2.1 and R-2.3.0

2006-04-27 Thread roger bos
I don't know if this is causing the error, by in your traceback I saw sort() was used and sort now removes all the attributes when it sorts. I used to use sort() to sort dates in character format and now it turns them into integers and is breaking all my code. The problem with upgrading is your

Re: [R] summary(lm(x~y)) difference between R-2.2.1 and R-2.3.0

2006-04-27 Thread Robin Hankin
Hi Roger good point. I can reproduce the error much more simply: x - 1:10 y - 10:1 summary(lm(x~y)) Call: lm(formula = x ~ y) Residuals: *** caught bus error *** address 0x18, cause 'invalid alignment' Traceback: 1: sort(x, partial = unique(c(lo, hi))) 2: quantile.default(resid) 3:

Re: [R] summary(lm(x~y)) difference between R-2.2.1 and R-2.3.0

2006-04-27 Thread roger bos
Robin, I still don't know what the problem is, but your example works on my R-2.3.0pat (created on 4/24/2006). You may want to try the patched version if you are not already running it. Thanks, Roger x - 1:10 y - 10:1 summary(lm(x~y)) Call: lm(formula = x ~ y) Residuals:

Re: [R] summary(lm(x~y)) difference between R-2.2.1 and R-2.3.0

2006-04-27 Thread Martin Maechler
roger == roger bos [EMAIL PROTECTED] on Thu, 27 Apr 2006 11:06:36 -0400 writes: roger Robin, I still don't know what the problem is, but roger your example works on my R-2.3.0pat (created on roger 4/24/2006). You may want to try the patched version roger if you are not

[R] summary stats

2006-04-16 Thread Neil Hepburn
I have a data set that has student test scores along with several categorical variables. I would like to generate a set of summary stats (mean, variance, n) for the data grouped by school authority and by exam topic. I have tried the by() function but that seems to only be able to handle one

Re: [R] summary stats

2006-04-16 Thread Christos Hatzis
the results: cbind(agg.mean, agg.var, agg.n) -Christos -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Neil Hepburn Sent: Sunday, April 16, 2006 2:56 PM To: r-help@stat.math.ethz.ch Subject: [R] summary stats I have a data set that has student test scores along

Re: [R] summary stats

2006-04-16 Thread John Fox
: Sunday, April 16, 2006 1:56 PM To: r-help@stat.math.ethz.ch Subject: [R] summary stats I have a data set that has student test scores along with several categorical variables. I would like to generate a set of summary stats (mean, variance, n) for the data grouped by school authority

[R] summary(lm) * lack of a parameter estimate

2006-04-11 Thread Luis Ridao Cruz
R-help, I'm using a lm model for some data. The code is below: fitData - lm(formula = log(f1) ~ year + mon + pair + rek, weights = w, data = data) I want to get the fitted values of log(f1) for every year so I do: tapply(fitted(fitData), list(data96$year), mean) but there is no parameter

[R] summary tables with tapply

2006-03-11 Thread Stefanie von Felten, IPWIfU
Hello everyone, #I have a long table with a factor Plotfac: Plotfac[1:20] [1] 1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 5 5 6 6 #there is the species diversity of each plot (Plotfac): DataSort$div[1:20] [1] 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 2 2 2 2 #and the

Re: [R] summary tables with tapply

2006-03-11 Thread Uwe Ligges
Stefanie von Felten, IPWIfU wrote: Hello everyone, #I have a long table with a factor Plotfac: Plotfac[1:20] [1] 1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 5 5 6 6 #there is the species diversity of each plot (Plotfac): DataSort$div[1:20] [1] 4 4 4

Re: [R] summary tables with tapply

2006-03-11 Thread Gabor Grothendieck
I am not sure I understand the question but is the situation that you have two vectors: x and y such that for each level of y x is constant so that for each level of y you want to find that value of x? In that case: x - c(A, A, A, B, B) y - c(1,1,2,3,3) unique(data.frame(x,y)) or tapply(x, y,

[R] Summary()

2006-03-02 Thread Dan Chan
Hi, I have a dataframe(FireDanger) that contains weather data (T, RH, WS, etc) measured at 27 stations. I want to do a quantile summary for each variable at each station. quantile(FireDanger$RH, probs=seq(0,1,0.1), na.rm=T) summary(FireDanger) are close, but not quite. Is there a single

Re: [R] Summary()

2006-03-02 Thread Liaw, Andy
Probably something like this: x - as.data.frame(matrix(rnorm(27 * 5), 27, 5)) sapply(x, quantile, seq(0, 1, .1)) V1 V2 V3 V4 V5 0% -3.18125816 -2.3857797 -1.0564583 -2.5085757 -2.142863890 10% -1.63832807 -1.3854193 -0.7863538 -2.0054480

[R] SUMMARY: aggregate vs. tapply

2006-02-12 Thread Joseph LeBouton
Hi all; Thanks for the responses to my query of how to make tapply into a table instead of an n-dimensional array. Summary of responses follows: Peter Dalgaard: as.data.frame(with(tmp,as.table(tapply(C,list(A=A,B=B),sum Phil Spector wrote: z = tapply(y,list(var1,var2,var3,var4),sum)

[R] Summary functions to dataframe

2006-01-03 Thread Mike Bock
I have written a few different summary functions. I want to calculate the statistics by groups and I am having trouble getting the output as a dataframe. I have attached one example with a small dataset that calculates summary stats and percentiles, I have others that calculate upper confidence

Re: [R] Summary functions to dataframe

2006-01-03 Thread Gabor Grothendieck
Try this: Pstats - function(x) c(Max = max(x), Min = min(x), AMean = mean(x), AStdev = sd(x), Samples = length(x), quantile(x, 1:9/10, na.rm = TRUE)) res - with(areas, by(AdRes, N_Type, Pstats)) do.call(rbind, res) Also, check out summaryBy in the doBy package at

Re: [R] summary[[r.squared]] gives strange results

2005-12-07 Thread Sundar Dorai-Raj
Felix Flory wrote: I am simulating an ANOVA model and get a strange behavior from the summary function. To be more specific: please run the following code and see for yourself: the summary()[[r.squared]] values of two identical models are quite different!! ## 3 x 3 ANOVA of two factors x

[R] summary nls output

2005-09-25 Thread Chun-Ying Lee
Dear R user: I bulid a package, and in the package I use the function nls to solve some questions. If I have two sets of data, and I want to summary these two data's nls output, I write the command in the package source code like: { .. summary(fm1) summary(fm2) } then i compiler the

Re: [R] summary nls output

2005-09-25 Thread Prof Brian Ripley
On Sun, 25 Sep 2005, Chun-Ying Lee wrote: Dear R user: I bulid a package, and in the package I use the function nls to solve some questions. If I have two sets of data, and I want to summary these two data's nls output, I write the command in the package source code like: { ..

Re: [R] summary nls output

2005-09-25 Thread Andrew Robinson
This is because R only permits one object to be output. If you weould like two objects, you might try something like { .. list(fm1 = summary(fm1), fm2 = summary(fm2)) } Good luck, Andrew On Sun, Sep 25, 2005 at 06:49:26PM +0800, Chun-Ying Lee wrote: Dear R user: I bulid a

Re: [R] summary of problem with mCall function.

2005-09-10 Thread Douglas Bates
On 9/7/05, Bill Shipley [EMAIL PROTECTED] wrote: Last week I posted a question concerning the mCall function, which is used to create self-starting functions and is described in the book by Pinheiro, J.C. and Bates, D.M. (Mixed-effects models in S and S-PLUS). On page 345 one finds the

[R] Summary: Unexpected result of read.dbf

2005-08-19 Thread Susumu Tanimura
Hi there, This is summary and patch for a bug in read.dbf, demonstrating in Message-Id: [EMAIL PROTECTED]. After consulting Rjpwiki, a cyber-community of R user in Japan, the cause was found, and the patch of solution was proposed. Overflowing occurs when we use read.dbf for reading a dbf

Re: [R] Summary: Unexpected result of read.dbf

2005-08-19 Thread Prof Brian Ripley
It really isn't clear that this is correct. The reason is correct: read.dbf treats numeric files with no decimals as integers, and that _is_ as stated on the help page. So it is definitely not a `bug', and reading the help would have shown the reason for the original question. [I in general

[R] summary: cbind a list of matrices

2005-07-16 Thread Mauro Gasparini
Thanks to all kind people who answered. Attached is a useful reply, among many others. -- Mauro Gasparini Professore Straordinario di Statistica Dipartimento di Matematica, Politecnico di Torino Corso Duca degli Abruzzi 24 I-10129 Torino, Italy tel: +39 011 564 7546 fax: +39 011 564 7599

[R] summary statistics for lists of matrices or dataframes

2005-05-10 Thread David Kane
Is there a simple way to calculate summary statistics for all the matrices or dataframes in a list? For example: z - list(matrix(c(2,2,2,2), ncol = 2), matrix(c(4,4,4,4), ncol = 2)) z [[1]] [,1] [,2] [1,]22 [2,]22 [[2]] [,1] [,2] [1,]44 [2,]44 I

RE: [R] summary statistics for lists of matrices or dataframes

2005-05-10 Thread McGehee, Robert
AM To: r-help@stat.math.ethz.ch Subject: [R] summary statistics for lists of matrices or dataframes Is there a simple way to calculate summary statistics for all the matrices or dataframes in a list? For example: z - list(matrix(c(2,2,2,2), ncol = 2), matrix(c(4,4,4,4), ncol = 2)) z [[1

Re: [R] summary statistics for lists of matrices or dataframes

2005-05-10 Thread Dimitris Rizopoulos
, May 10, 2005 4:03 PM Subject: [R] summary statistics for lists of matrices or dataframes Is there a simple way to calculate summary statistics for all the matrices or dataframes in a list? For example: z - list(matrix(c(2,2,2,2), ncol = 2), matrix(c(4,4,4,4), ncol = 2)) z [[1]] [,1] [,2] [1

Re: [R] summary statistics for lists of matrices or dataframes

2005-05-10 Thread Patrick Burns
You could use 'do.call' with 'bind.array' (from S Poetry) or 'abind' to convert your list of matrices into a three-dimensional array. Patrick Burns Burns Statistics [EMAIL PROTECTED] +44 (0)20 8525 0696 http://www.burns-stat.com (home of S Poetry and A Guide for the Unwilling S User) David Kane

Re: [R] summary statistics for lists of matrices or dataframes

2005-05-10 Thread Patrick Burns
No, you can't. Because 'bind.array' doesn't take an arbitrary number of arguments. Robert's solution does what I had in mind. Patrick Burns wrote: You could use 'do.call' with 'bind.array' (from S Poetry) or 'abind' to convert your list of matrices into a three-dimensional array. Patrick Burns

RE: [R] summary statistics for lists of matrices or dataframes

2005-05-10 Thread Huntsinger, Reid
] [mailto:[EMAIL PROTECTED] On Behalf Of David Kane Sent: Tuesday, May 10, 2005 10:04 AM To: r-help@stat.math.ethz.ch Subject: [R] summary statistics for lists of matrices or dataframes Is there a simple way to calculate summary statistics for all the matrices or dataframes in a list? For example: z

[R] Summary: My question about factor levels versus factor labels.

2005-05-09 Thread Ajay Narottam Shah
Yesterday, I had asked for help on the list. Brian Ripley and Bruno Falissard had most kindly responded to me. Here is the solution. factorlabels - c(School, College, Beyond) # 1 2 3 education.man - c(1,2,1,2,1,2,1,2) # PROBLEM: Level 3 doesn't

[R] summary(as.factor(x) - force to not sort the result according factor levels

2005-05-02 Thread Christoph Lehmann
Hi The result of a summary(as.factor(x)) (see example below) call is sorted according to the factor level. How can I get the result not sorted but in the original order of the levels in x? test - c(120402, 120402, 120402, 1323, 1323,200393, 200393, 200393, 200393, 200393)

Re: [R] summary(as.factor(x) - force to not sort the result according factor levels

2005-05-02 Thread Tony Plate
Christoph Lehmann wrote: Hi The result of a summary(as.factor(x)) (see example below) call is sorted according to the factor level. How can I get the result not sorted but in the original order of the levels in x? by creating the factor with the levels in the order you want: test - c(120402,

Re: [R] summary(as.factor(x) - force to not sort the result according factor levels

2005-05-02 Thread Sundar Dorai-Raj
Christoph Lehmann wrote on 5/2/2005 3:29 PM: Hi The result of a summary(as.factor(x)) (see example below) call is sorted according to the factor level. How can I get the result not sorted but in the original order of the levels in x? test - c(120402, 120402, 120402, 1323, 1323,200393,

[R] Summary: GLMMs: Negative Binomial family in R

2005-04-13 Thread nflynn
Here is a summary of responses to my original email (see my query at the bottom). Thank you to Achim Zeileis , Anders Nielsen, Pierre Kleiber and Dave Fournier who all helped out with advice. I hope that their responses will help some of you too. *

[R] summary method in URCA package doesn't work

2005-02-25 Thread Chalasani, Prasad
I can't figure out how to get the summary method in the URCA package to work. E.g. when I use the following code fragment in the help for the ca.jo function, it always tries to use the summary method from the base package, not the urca package. How do I force it use the summary method of the

Re: [R] summary method in URCA package doesn't work

2005-02-25 Thread Achim Zeileis
Please contact (or at least Cc) the maintainer of contributed packages when reporting problems with contributed packages. (I've added Bernhard to Cc: now.) I can't figure out how to get the summary method in the URCA package to work. E.g. when I use the following code fragment in the help for

Re: [R] summary(aov(...)) into a string?

2005-02-15 Thread Prof Brian Ripley
It doesn't print anything: the summary.aov (or summary.aovlist) print method does. ?summary.aov tells you the structure of the objects they return. On Tue, 15 Feb 2005, RenE J.V. Bertin wrote: I'd like to annotate a plot with the output of summary(aov(model)), ideally just with the significant

Re: [R] summary(aov(...)) into a string?

2005-02-15 Thread Gabor Grothendieck
RenE J.V. Bertin rjvbertin at hotmail.com writes: : : Hello, : : I'd like to annotate a plot with the output of summary(aov(model)), ideally just with the significant : effects. I don't find a means to redirect what that command prints into a string. : Is this possible, and if so, how? :

Re: [R] summary(aov(...)) into a string?

2005-02-15 Thread Peter Dalgaard
Prof Brian Ripley [EMAIL PROTECTED] writes: It doesn't print anything: the summary.aov (or summary.aovlist) print method does. ?summary.aov tells you the structure of the objects they return. Yes. but wouldn't capture.output() and par(family=mono) be closer to the mark? On Tue, 15

[R] summary output for inverse Gaussian GLM

2004-08-10 Thread Edward Dick
I'm getting inconsistent output about the link function from the summary() command when fitting an inverse Gaussian GLM: summary(glm(ig.formula, family=inverse.gaussian(link = log), + data=mydata, start=start.vals))$call glm(formula = ig.formula, family = inverse.gaussian(link = log), data =

Re: [R] summary output for inverse Gaussian GLM

2004-08-10 Thread Thomas Lumley
On Tue, 10 Aug 2004, Edward Dick wrote: I'm getting inconsistent output about the link function from the summary() command when fitting an inverse Gaussian GLM: Yes, it's a bug in the name, not in the result, fortunately. The inverse.gaussian family object always has the name of the link set

Re: [R] summary() doesn't work with Date class objects

2004-07-14 Thread Uwe Ligges
Scott Waichler wrote: The handy function summary() doesn't work correctly with Date class objects: R.version.string [1] R version 1.9.1, 2004-06-21 b - as.Date(c(2002-12-26, 2002-12-27, 2002-12-28, 2002-12-29, 2002-12-30)) b [1] 2002-12-26 2002-12-27 2002-12-28 2002-12-29 2002-12-30 summary(b)

[R] summary() doesn't work with Date class objects

2004-07-13 Thread Scott Waichler
The handy function summary() doesn't work correctly with Date class objects: R.version.string [1] R version 1.9.1, 2004-06-21 b - as.Date(c(2002-12-26, 2002-12-27, 2002-12-28, 2002-12-29, 2002-12-30)) b [1] 2002-12-26 2002-12-27 2002-12-28 2002-12-29 2002-12-30 summary(b) Min.

RE: [R] text editor for R - summary on R-WinEdt

2004-07-12 Thread Benjamin . STABLER
] Subject: Re: [R] text editor for R - summary on R-WinEdt Dear all, let me try to summarize this thread's R-WinEdt related messages and give a few comments: Murray Jorgensen wrote: I tried R-WinEdt a few years ago, but as I remember it interfered with my usual use of WinEdt which

Re: [R] text editor for R - summary on R-WinEdt

2004-07-08 Thread Uwe Ligges
Dear all, let me try to summarize this thread's R-WinEdt related messages and give a few comments: Murray Jorgensen wrote: I tried R-WinEdt a few years ago, but as I remember it interfered with my usual use of WinEdt which is as a front end to MiKTeX. Is there a way to use WinEdt both ways?

[R] SUMMARY: elementary sapply question

2004-06-22 Thread Ajay Shah
I am grateful to Andy Liaw, Douglas Grove, Brian Ripley, Tony Plate, Dirk Eddelbuettel and Sundar Dorai-Raj all of whom got together and drilled sense into my skull. I would like to take some effort into explaining what the question was, that I was grappling with, and the (nice) R way of solving

[R] summary table newbie question

2004-05-12 Thread Jeff D. Hamann
I've got a newbie question and I got a little lost in the table helps. I've got a data.frame I would like to summarize as a (and pardon for the lack of correct vernacular) data collection matrix. My data looks like, stand siteindex age acres pct.acres 1232 116 45 8477.3105

  1   2   >