On Wed, 4 Aug 2004, Duncan Mackay wrote:
> Would someone like to comment on the relation between the "buffer bytes"
> and "lines" settings in the console settings? Are they interdependent?
Both are limits, and you need enough of each.
> If so, how? And why allow both to be tweaked? If not, what
On Tue, 3 Aug 2004, Thomas Lumley wrote:
> I can't remember the historical reason why utils isn't loaded in binary
> form in 1.9.1. In r-devel they are all loaded from binary databases when
> used.
It's size. The R code for utils is relatively small (relative to stats
and graphics, or lattice
Hoi Tiago,
--On dinsdag 3 augustus 2004 18:50 +0100 Tiago R Magalhaes
<[EMAIL PROTECTED]> wrote:
Once again I'm sorry for these basic questions and since predictably
I'll have more of those if there's a basic-questions-list I would
love to know more about it
Recently we discussed, on this list, t
Thanks to all who helped.
I used your ideas and code samples to write the following (for the
benefit of people who will search this list later):
rbind.case <- function(..., name="case", values) {
dots <- list(...);
if (missing(values)) values <- 1:length(dots);
if (length(values)!=length
I am curious as to why I quite often receive responses to questions on the help before
receiving the actual question.
For example Mr Graves has responded to Mr Lumley's response to Mr Wegelin yet I ave
not received Mr. Lumley's post while I have received the other two.
Cheers,
Jim
[
You may also wish to consider "VarCorr". This function exists in
both nlme and lme4. Someone else mentioned "getVarCov". However, that
is not in the current lme4.
hope this helps. spencer graves
Thomas Lumley wrote:
On Tue, 3 Aug 2004, Jacob Wegelin wrote:
The print method
Hi all,
Would someone like to comment on the relation between the "buffer bytes"
and "lines" settings in the console settings? Are they interdependent?
If so, how? And why allow both to be tweaked? If not, what is the effect
of changing "lines", but leaving "buffer bytes" the same?
I wanted to be
I believe getVarCov() will give you what you need.
Harold
-Original Message-
From: [EMAIL PROTECTED] on behalf of Thomas Lumley
Sent: Tue 8/3/2004 8:36 PM
To: Jacob Wegelin
Cc: [EMAIL PROTECTED]
Subject: Re: [R] lme fitted correlation
On Tue, 3 Aug 2004, Jacob Wegelin wrote:
>
> The print method for lme *prints out* the fitted correlation matrix for
> the random effects. Is there any way to get these values as an object in
> R? I have examined the components of the lme object (called "junk" in the
> example below) and the comp
Dear Prof Ripley
Thanks for your reply and clarification. However:
1. Regarding model.tables() returning "Design is unbalanced". Setting contrasts to
Helmert does indeed make the design balanced, but model.tables() still returns "Design
is unbalanced":
> options()$contrasts
unorder
Thanks for the tip.
-F.
On Aug 3, 2004, at 1:21 AM, Prof Brian Ripley wrote:
On Tue, 3 Aug 2004, Francisco J. Bido wrote:
I have R 1.9.1 on MacOSX. The automated package download appears
broken in this version.
Your version? It is not broken in the version R-core distributes, but
you
may be usin
The print method for lme *prints out* the fitted correlation matrix for
the random effects. Is there any way to get these values as an object in
R? I have examined the components of the lme object (called "junk" in the
example below) and the components of summary(junk) without finding these
numbe
On Tue, 3 Aug 2004, Rolf Turner wrote:
>
> Thomas Lumley wrote:
>
> > Because, as you noted, the base packages are stored in binary form. This
> > already speeds things up, and will have even more impact in 2.0.0 with
> > "lazy loading" of functions.
>
> This CANNOT be the (complete) explan
Thomas Lumley wrote:
> > However attr(arima,"source") was still NULL (and of course comments
> > were non-existant). Apparently the functions in stats get
> > ***loaded*** from a binary file ``all.rda'' (in the
> > .../library/stats/R directory) rather than getting read in from ascii
> > files.
Thomas Lumley <[EMAIL PROTECTED]> writes:
> Personally, I find this the easiest way to look at the source. You can
> browse the source on https://svn.r-project.org/R/ if you are using a
> binary distribution, rather than downloading the whole thing.
Except that you probably don't want to, i.e. gr
On Tue, 3 Aug 2004, Thomas Lumley wrote:
> On Tue, 3 Aug 2004, Jack Tanner wrote:
>
> > > write.table(NULL)
> > Error in which(unlist(lapply(x, function(x) is.character(x) ||
> > is.factor(x : argument to "which" is not logical
> >
> > Is this correct behavior? It seems harsh to abort an ent
> > On Mon, 2 Aug 2004 17:50:44 -0500, "Bickel, David"
> > <[EMAIL PROTECTED]> wrote:
> >
> > >After reading the help page on set.seed, I am unsure about
> > how to solve the following problem. I need to call function f
> > a thousand times. The list of values returned by f, should be
> > as random
On Tue, 3 Aug 2004, Jack Tanner wrote:
> > write.table(NULL)
> Error in which(unlist(lapply(x, function(x) is.character(x) ||
> is.factor(x : argument to "which" is not logical
>
> Is this correct behavior? It seems harsh to abort an entire run just
> because one of the tables you generated h
On Tue, 3 Aug 2004, Rolf Turner wrote:
>
> However, after substantial experimentation, I am still at a loss to
> get at the original source code for arima() (replete with comments)
> save by accessing the original, ``unmade'' source for R. Not too
> hard, but not exactly ``immediate''.
Personally
> write.table(NULL)
Error in which(unlist(lapply(x, function(x) is.character(x) ||
is.factor(x : argument to "which" is not logical
Is this correct behavior? It seems harsh to abort an entire run just
because one of the tables you generated happened to be NULL.
-JT
_
Recently Brian Ripley had occasion to mock my inability to see a
comment in the code for arima(), in the stats package. After
considerable dredging around in the r-news archives I found reference
to keep.source() and keep.source.pkgs(), which I conjectured just
***might*** possibly be the ``obvio
It is in MASS4, and even in the R scripts for MASS3 in MASS/scripts3 in
your R installation.
And that is explained in the R complements on the book's web site.
On Tue, 3 Aug 2004, Sven C. Koehler wrote:
> Dear list!
>
> I am interested in learning about MLE and I wonder whether it is possible
Graphically, you can use coplot() or lattice to see how the relationship
between height and weight changes with age.
With a linear model, you can do something like:
lm(height ~ weight * age, ...)
and work out the interpretation of the coefficients.
HTH,
Andy
> From: Peter Tait
>
> Hi,
>
> I
As you know, R is not S-Plus, though they for the most part share a common syntax
(S) and a lot of (nearly) identical functions with (nearly) identical behavior.
Nevertheless, there's lots of stuff in S-Plus that's not in R and vice-versa. So
you shouldn't expect to find all S-Plus functions in R.
Hi,
I want to look at the relationship between 2 continuous/quantitative variables while
adjusting for a third continuous variable. For example: relationship between height
and weight adjusted for age.
As I understand it using a glm would work if I had a categorical variable (I could
make age i
> From: Tiago R Magalhaes
>
> Hi
> I have two basic questions, and here they go, but I was wondering as
> well where can I ask these basic questions without bothering you
> people
>
>
> I've used Splus and now I'm using R and there's some functions that I
> can't simply find
>
> one:
> sort
Hi
I have two basic questions, and here they go, but I was wondering as
well where can I ask these basic questions without bothering you
people
I've used Splus and now I'm using R and there's some functions that I
can't simply find
one:
sort.col
that allows data.frames to be sort by a given
Dear list!
I am interested in learning about MLE and I wonder whether it is possible
to use the examples for maximum likelihood estimation given in 8.7 in
MASSv3 with R? AFAIU R does not have a direct replacement for S-PLUS's
ms() which the examples use for the fitting, but optim() may be of help
My guess is that
http://cran.r-project.org/src/contrib/Descriptions/rlecuyer.html might be
relevant, too.
Cheers,
Andy
> From: Duncan Murdoch
>
> On Mon, 2 Aug 2004 17:50:44 -0500, "Bickel, David"
> <[EMAIL PROTECTED]> wrote:
>
> >After reading the help page on set.seed, I am unsure about
> ho
Doran, Harold air.org> writes:
:
: I am working with a longitudinal data set in the long format. This data
: set has three observations per grade level per year. Here are the first
: 10 rows of the data frame:
:
: >tenn.dat[1:10,]
:
: year schid type grade gain se new cohort
:
: 6 2001 100
Have you considered "cbind" and "rbind"? If your data.frame has
factors, they could present problems with "rbind". Try
'sapply(tenn.dat, class)'. If you have only class character or numeric,
use "cbind" to match the columns of tenn.dat both in name and class,
then use "rbind" to combine it w
I am working with a longitudinal data set in the long format. This data
set has three observations per grade level per year. Here are the first
10 rows of the data frame:
>tenn.dat[1:10,]
year schid type grade gain se new cohort
6 2001 155 4 33.1 3.5 4 3
7 2002 1000
On Tue, 3 Aug 2004, Eve Mathieu wrote:
> Dear R-help group,
>
> I have to maximize a likelihood with 40 parameters and I want to compare
> the MLE given by "nlminb" (Splus2000, on Windows) with those given by
> "optim" (R, on Unix).
>
> 1) On Splus,
> The algorithm "nlminb" seems to converge (the
On Tuesday 03 August 2004 06:39, 8rino-Luca Pantani wrote:
> Hi all,
> first of all thanks for the answer to my previous question on lattice
> some time ago.
> In particular to Patrick Connolly for advices on netiquette (I hope
> this time I'm doing right)
> and to Deepayan Sarkar fro the help
hello.
I have a query about the Matrix package for R. I wrote some code a while
ago using Matrix version 1.6.2 with an early version of R, to do some least
squares for regression:
xn
[,1] [,2] [,3] [,4]
[1,]1 0.7429352 0.5519528 0.4100652
[2,]1 0.7443713 0.5540886 0.4
Duncan already gave the answer. Here's a bit more detail:
rbindGroup <- function(...) {
theList <- list(...)
nGroup <- length(theList)
size <- sapply(theList, nrow)
result <- rbind(...)
result <- cbind(result, rep(1:nGroup, each=size)
result
}
You probably want to add mor
Dear R-help group,
I have to maximize a likelihood with 40 parameters and I want to compare
the MLE given by "nlminb" (Splus2000, on Windows) with those given by
"optim" (R, on Unix).
1) On Splus,
The algorithm "nlminb" seems to converge (the parameters stabilize) , it
stops after several itera
Hi all,
first of all thanks for the answer to my previous question on lattice some
time ago.
In particular to Patrick Connolly for advices on netiquette (I hope this
time I'm doing right)
and to Deepayan Sarkar fro the help on lattice.
Now, my nowaday problem.
Please consider the following
m
Hello,
I've a problem aligning tickmarks to an image. I've created a correlation matrix for
84 datasets. I'm visualizing the matrix as an image with colour coding according to
the correlation coefficient.
The 84 datasets are distributed over three factors, but the desgin is unbalanced, so
that
On Tue, 3 Aug 2004 15:51:43 +1000, "Briggs, Meredith M"
<[EMAIL PROTECTED]> wrote:
>
>Hello
>
>I have to show that R has more functionality than SAS for MonteCarlo type exercises
>involving matrix algebra. As I has little experience with SAS is anyone able to help.
>I have a model written in R b
> "Sundar" == Sundar Dorai-Raj <[EMAIL PROTECTED]>
> on Mon, 02 Aug 2004 12:53:43 -0500 writes:
Sundar> F Duan wrote:
>> I am using par(mfrow=c()) to plot multi-figures in the
>> same window. And I like to put a common title (and xlab,
>> ylab) for all of plots. I have
Hoi Jinsong,
--On dinsdag 3 augustus 2004 1:42 -0700 Jinsong Zhao <[EMAIL PROTECTED]>
wrote:
For instance, I hope to remove the V3~6 column, for all the value in
those colume is zero.
V3 V4 V5 V6 V7 V8 V9V10
1 0 0 0 0 0.000 0.000 0.000 0.000
2 0 0 0 0 0.000 0.000 0
On Tue, 3 Aug 2004, Jinsong Zhao wrote:
> I hope to remove a whole column from a data frame or matrix (> 2000
> columns). All value in the column are same. The first thing is to
> select those columns.
>
> For instance, I hope to remove the V3~6 column, for all the value in
> those colume is zer
What about:
my.df<-data.frame(x=rnorm(20),y=rep(0,20),z=rep(1,20))
my.df[,!apply(my.df,2,sd)==0,drop=F]
This uses the fact that a column with eaqual values has a standard deviation
of 0.
Regards
Wayne
-Original Message-
From: Jinsong Zhao [mailto:[EMAIL PROTECTED]
Sent: 03 Augu
Dear all,
I hope to remove a whole column from a data frame or matrix (> 2000
columns). All value in the column are same. The first thing is to
select those columns.
For instance, I hope to remove the V3~6 column, for all the value in
those colume is zero.
V3 V4 V5 V6 V7 V8 V9
Hi,
see ? quantile to obtain deciles of variable X1
see ? cut to divide the range of 'x' into intervals
and codes the values in 'x' according to which
interval they fall.
se ? table to use the cross-classifying factors to
build a contingency table of the counts at each
combination of factor leve
46 matches
Mail list logo