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
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
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
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
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
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
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
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
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 =
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,
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
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
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),
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
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
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
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
'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
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
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
...
-- 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
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.
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)
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
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
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
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
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
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
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.
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)
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,
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
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
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
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
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
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
: 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
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
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
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
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
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
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
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
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
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
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
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).
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
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,
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
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:
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:
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
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
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
: 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-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
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
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
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,
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
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
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)
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
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
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
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
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:
{
..
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
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
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
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
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
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
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
, 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
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
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
]
[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
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
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)
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,
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,
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.
*
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
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
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
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?
:
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
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 =
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
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)
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.
]
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
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?
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
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 - 100 of 157 matches
Mail list logo