Iain Pardoe said the following on 2005-05-04 20:08:
I'd like to model the relationship between m responses Y1, ..., Ym and a
single set of predictor variables X1, ..., Xr. Each response is assumed
to follow its own regression model, and the error terms in each model
can be correlated. My
Paul Smith wrote:
Dear All
I would like to know whether it is possible with R to get the
mathematical expression of the density of a sum of two independent
continuous random variables.
No, that corresponds to a convolution of the two densities, and R can't
do any symbolic integration.
You could
Have you considered package distr? It will do something similar to
what you request, I think; it may or may not be adequate for your
purposes.
spencer graves
Duncan Murdoch wrote:
Paul Smith wrote:
Dear All
I would like to know whether it is possible with R to get the
Angus Repper [EMAIL PROTECTED] writes:
Hello
I am a long-time SAS user, but am new to R. I specifically am looking for
information pertaining to generating graphics for web output. I would like
to create dynamic graphics (in the form of generalized reports) for my web
site that is
Dear,
I am new user of 'R' (package e1071) and the first
thing I did was reading the manuals. Although it was
not clear for me how it works.
My problem is the following:
I have a data set of customers data of insurances with
variables Class, Group Distance Claim. The file
called 'xx.dat'. Using
Dear list,
I do not know if this is a bug.
Let's suppose mat is a matrix derived from this code:
x-rnorm(10)
y-rnorm(10)
names(x)-LETTERS[1:10]
names(y)-LETTERS[1:10]
mat-cbind(x,y)
Now in R 2.0.1 I have:
names(dist(mat))
A B C D E F G H I J
In R 2.1.0 I have:
names(dist(mat))
NULL
I have tried
Elio Mineo [EMAIL PROTECTED] writes:
In R 2.1.0 I have:
names(dist(mat))
NULL
I have tried this on Windows (precompiled binary distribution) and on
Linux Mandrake 10.0 (R compiled from source code).
Is this a bug or not?
Not.
From the NEWS file:
o The S3 methods for getting
Dear All
I would like to get the histogram for the following model with
discrete and continuous random variables:
* with probability 1/3, a random number is drawn from the continuous
uniform distribution (min=0, max=1);
* with probability 2/3, a random number is drawn from a different
continuous
On 5/5/05, Spencer Graves [EMAIL PROTECTED] wrote:
Have you considered package distr? It will do something similar
to
what you request, I think; it may or may not be adequate for your
purposes.
Thanks, Spencer and Duncan. Maybe, the best choice is to use Maple or
MuPAD for that.
See http://finzi.psych.upenn.edu/R/Rhelp02a/archive/48428.html
Andy
From: Paul Smith
Dear All
I would like to get the histogram for the following model with
discrete and continuous random variables:
* with probability 1/3, a random number is drawn from the continuous
uniform
I
heard that 'R' does not do a very good job at handling large datasets, is
this true?
importing huge datasets in a data.frame with e.g. a subsequent step of
conversion of some columns into factors may lead into memory troubles
(probably due to memory overhead when building out factors). But
Thanks, Peter.
So, in R 2.1.0 there is no way to get names from a dist object, isn't
there?
Peter Dalgaard wrote:
Elio Mineo [EMAIL PROTECTED] writes:
In R 2.1.0 I have:
names(dist(mat))
NULL
I have tried this on Windows (precompiled binary distribution) and on
Linux Mandrake 10.0
You gave no detail for anyone to provide useful help. What exact commands
did you use, and what, if any, error/warning did you get? Please (re-)read
the posting guide. You may also want to do
vignette(svmdoc)
Andy
From: Dominique Knockaert
Dear,
I am new user of 'R' (package e1071)
From: Elio Mineo
Thanks, Peter.
So, in R 2.1.0 there is no way to get names from a dist
object, isn't
there?
If you are referring to the row names of the matrix from which distances are
computed, see the Value section of ?dist, which has:
The object has the following attributes (besides
On Thu, 5 May 2005, Christoph Lehmann wrote:
I heard that 'R' does not do a very good job at handling large
datasets, is this true?
importing huge datasets in a data.frame with e.g. a subsequent step of
conversion of some columns into factors may lead into memory troubles
(probably due to
Hi there. I have a data frame, X, with n+m columns. I want to regress each of
the first n columns on the last m. This is what I am trying to do:
for ( i in 1:n )
assign( paste (reg,1:14,sep=)[i] , lm( X[,i] ~ X[,i+1] + ... +
X[,i+m], data= X ) )
It happens that some of the
On 04-May-05 Roger Bivand wrote:
On Wed, 4 May 2005, Sean Davis wrote:
see ?aggregate.
Or maybe tapply, or its close relative, by:
by(df, list(df$station, df$date), function(x)
+ x$row[which.max(x$chlorophyll)])
: Ancona
: 21/06/01
[1] NA
Hi
I am applying the wilcoxon paired signed rank test to many sub-groups of
a data frame using by() and wilcox.test, and sometimes the test uses
continuity correction and sometimes it doesn't. What I want to know is:
1) what is continuity correction? The docs say:
correct: a logical
Hi,
I have a data.frame which looks like this:
id est0 est1est2est3est4est5est6est7
1 1 2 3 1 7 9 3 4
2 4 1 1 7 6 5 1 2
[...]
I would like to reorder it to obtain
Hi,
I have been returning to the same problem a number of times without success
and now look for help with the following:
How do I fit a distribution function with the same number of parameters as
there are quantiles and values such that I get an exact solution as opposed
to a minimum least
Thanks to Andy Liaw to suggest me the right solution (at least for me):
attributes(dist(mat))$Labels
or
attr(dist(mat),Labels)
Bye,
Elio
--
--
Prof. Angelo M. Mineo
Dipartimento di Scienze Statistiche e Matematiche S.
See ?reshape
Uwe Ligges
Achaz von Hardenberg wrote:
Hi,
I have a data.frame which looks like this:
id est0 est1 est2 est3 est4 est5 est6est7
1 1 2 3 1 7 9 3 4
2 4 1 1 7 6 5 1 2
[...]
I would like to reorder it to obtain the following:
id est VALUE
1 0 1
Dimitri Joe wrote:
Hi there. I have a data frame, X, with n+m columns. I want to regress each of
the first n columns on the last m. This is what I am trying to do:
for ( i in 1:n )
assign( paste (reg,1:14,sep=)[i] , lm( X[,i] ~ X[,i+1] + ... +
X[,i+m], data= X ) )
It happens that some of
Use reshape(), as:
dat
id est0 est1 est2 est3 est4 est5 est6 est7
1 112317934
2 241176512
dat2 - reshape(dat, varying=list(names(dat)[-1]), times=0:7,
direction=long)
dat2
id time est0
1.0 101
2.0 204
1.1
Hi
R d -
data.frame(id=1:3,est0=c(1,10,100),est1=c(6,7,8),est2=c(-5,-5,-5))
R d
id est0 est1 est2
1 116 -5
2 2 107 -5
3 3 1008 -5
R out - as.vector(t(d[,-1]))
R
cbind(id=rep(d$id,each=nrow(d)),est=as.vector(row(as.matrix(d[,
-1]))),out)
id est out
[1,]
On Fri, 2005-05-06 at 15:26 +0200, Achaz von Hardenberg wrote:
Hi,
I have a data.frame which looks like this:
id est0 est1 est2est3est4est5est6est7
11 2 3 1 7 9 3 4
24 1 1 7 6 5
On 5/3/05, Rolf Turner [EMAIL PROTECTED] wrote:
Sean Davis wrote:
I just did some of this learning myself. Here are a couple of links
that I found useful:
http://www.stat.auckland.ac.nz/S-Workshop/Gentleman/S4Objects.pdf
http://eeyore.ucdavis.edu/stat250/OOP.html
I found the
On 5/5/05, Ted Harding [EMAIL PROTECTED] wrote:
On 05-May-05 Peter Dalgaard wrote:
[...]
Both systems are victims of the curse of the rectangular data set to
some extent. Prototypically, you record the sex of a rat along with
every single measurement on it, as if the rat could change sex
Hello,
Is there any possibility in R to see the body of the non-visible
function, for
example princomp?
If I do :
methods(princomp)
so, I get that princomp.default and princomp.formula are non-visible
functions and
body(princomp.default) doesnt show it.
In particular, I guess I have a very
Here is an approach using optim:
tmpfunc - function(param){
ml-param[1]
sl-param[2]
(qlnorm(.15,ml,sl)-10)^2 + (qlnorm(.5,ml,sl)-30)^2
}
res - optim(c(1,2), tmpfunc)
res - optim(res$par, tmpfunc)
res
hope this helps,
Greg Snow, Ph.D.
Statistical Data Center, LDS
getS3method(princomp, default).
getAnywhere(princomp.default)
The first of these is in the See Also of ?methods.
On Thu, 5 May 2005, Anna Oganyan wrote:
Hello,
Is there any possibility in R to see the body of the non-visible
function, for
example princomp?
If I do :
methods(princomp)
so, I
Anna Oganyan wrote on 5/5/2005 7:42 AM:
Hello,
Is there any possibility in R to see the body of the non-visible
function, for
example princomp?
If I do :
methods(princomp)
so, I get that princomp.default and princomp.formula are non-visible
functions and
body(princomp.default) doesnt show
I have a few problems with pairs(). I want to add the variables names to the
plot, but I don't know how. I tried different things with the labels
argument but nothing worked out. This is the code:
read.table(bay1.bay,header=T)-Rdata
(histogram part skipped)
pairs(Rdata,
On Thu, 5 May 2005, Liaw, Andy wrote:
From: Anna Oganyan
Hello,
Is there any possibility in R to see the body of the non-visible
function, for
example princomp?
If I do :
methods(princomp)
so, I get that princomp.default and princomp.formula are non-visible
functions and
body(princomp.default)
On 5/5/05, Uwe Ligges [EMAIL PROTECTED] wrote:
See ?reshape
Uwe Ligges
Achaz von Hardenberg wrote:
Hi,
I have a data.frame which looks like this:
id est0 est1 est2est3est4est5est6est7
11 2 3 1 7 9 3 4
(Lazy answer, not checking your code in detail)
Try using splom() in the lattice package, as it gives more extensive control
to do the sorts of things you seem to want.
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the statistician is to catalyze the
Why does the following not work?
by(warpbreaks, warpbreaks$tension, function(x) add1(lm(breaks ~ 1,
data=x), ~ . + wool))
Error in inherits(x, data.frame) : Object x not found
Is there another efficient way to apply add1() to subsets of a data frame?
Thanks.
--
Chuck Cleland, Ph.D.
NDRI, Inc.
71
On Wednesday 04 May 2005 10:30, Sander Oom wrote:
Dear R graphics gurus,
Another question about lattice graphics. This time I would like to plot
means and confidence intervals by group factor in a lattice graph. I can
not find any working lattice examples. Maybe a custom panel function is
Could anyone help me with the following issue.
Using the GSL library in R I define the following code:
#
library(gsl);
S-function(n)
{ r-0:n;
ans-sum(gsl_sf_choose(n,r)*(-1)^r*2^(2*r)*gamma_inc(6-2*r,2))
ans }
#
SS(10) yields 34.91868
SS(40) yields 5.340422
SS(60) yields
I was wondering if someone can tell me the best way to search through a
matrix and choose certain rows (based on certain conditions) to put into a
separate matrix.
What I have tried so far is very slow for a large dataset I'm working with.
e.g., I have this piece of code to create a new matrix
I'm sure it's possible to display the variable names by making a small
change to the code. Someone did this last week but I can't contact this
person now to ask how. I tried this
vars - names(Rdata)
pairs(Rdata, labels=vars,...)
but that doesn't work (formal argument labels matched by multiple
For the following code below, the x-axis ticks are 1,2,3,4,5,6,7 when I was
expection them to be 1,2,8,9,10,11,12. Please help me figure out where is the
mistake.
library(lattice)
testdata - as.data.frame(t(structure(c(
1,2005,9.24,6.18,634,
2,2005,8.65,6.05,96,
8,2004,6.81,6.51,16,
[EMAIL PROTECTED] wrote:
I was wondering if someone can tell me the best way to search through a
matrix and choose certain rows (based on certain conditions) to put into
a separate matrix.
What I have tried so far is very slow for a large dataset I'm working
with. e.g., I have this piece of
The summands can get pretty large: the 2^(2*r) factor will erase a lot of
precision. Is there another way to compute this? Perhaps use the recurrence
relation for the binomial coefficient to get a recurrence relation for your
function?
Reid Huntsinger
-Original Message-
From: [EMAIL
Yikes, that looks complicated. I think what you're trying to do is available
already, see:
?subset
--
Tyler Smith
PhD Candidate
Plant Science Department
Faculty of Agricultural and Environmental Sciences
McGill University
MacDonald Campus
21,111 Lakeshore
Ste. Anne de Bellevue, Quebec, Canada
If you want variable names, why do you have labels = NULL in your call??
Your highlighting criterion is not clear. If you want to choose one variable
and highlight points in all plots for which the value of that variable meets
your criterion, then determine the indices (outside of pairs) for
Ghosh, Sandeep wrote on 5/5/2005 11:10 AM:
For the following code below, the x-axis ticks are 1,2,3,4,5,6,7 when I was
expection them to be 1,2,8,9,10,11,12. Please help me figure out where is the
mistake.
library(lattice)
testdata - as.data.frame(t(structure(c(
1,2005,9.24,6.18,634,
On Thursday 05 May 2005 13:10, Ghosh, Sandeep wrote:
For the following code below, the x-axis ticks are 1,2,3,4,5,6,7 when I was
expection them to be 1,2,8,9,10,11,12. Please help me figure out where is
the mistake.
[...]
colnames(testdata) - c('month', 'year', 'mean','stdDev','miceCount')
Hi Sundar,
Greatly appreicate your response, but I get this error on running the barchart
cmd
with(testdata, print(barchart(mean ~ month | year, data = testdata,
layout = c(1, nlevels(testdata$year)),
horizontal = FALSE,
ylim = c(1, max(testdata$mean) +
Deepayan Sarkar wrote on 5/5/2005 12:03 PM:
On Thursday 05 May 2005 13:10, Ghosh, Sandeep wrote:
For the following code below, the x-axis ticks are 1,2,3,4,5,6,7 when I was
expection them to be 1,2,8,9,10,11,12. Please help me figure out where is
the mistake.
[...]
colnames(testdata) - c('month',
The following line when pasted into an R Console window causes Windows XP to flash a blue screen
and then restart.
R version 2.1.0
plot(x=c(1:10),y=sin(c(1:10)), type=l)
Windows XP SP2
I installed the precompiled version of R 2.1.0
This isn't the data I was originally trying to graph,
From: Pieter Provoost
I'm sure it's possible to display the variable names by making a small
change to the code. Someone did this last week but I can't
contact this
person now to ask how. I tried this
vars - names(Rdata)
pairs(Rdata, labels=vars,...)
but that doesn't work (formal
It worked for me..Here are the commands..
library(lattice)
testdata - as.data.frame(t(structure(c(
1,2005,9.24,6.18,634,
2,2005,8.65,6.05,96,
8,2004,6.81,6.51,16,
9,2004,9.0,7.29,8,
10,2004,8.84,6.18,524,
11,2004,8.54,6.35,579,
12,2004,9.97,6.3,614,
12,2005,8.75,5.84,32,
), .Dim=c(5,8
Deepayan,
Thank you for your response. I tried your suggestion to manually specify the
y-axis label for the factor but I am stuck. My attempt is below. I once again
appreciate any input as well as your time. lattice/R-project is great.
Thanks once again,
David Arenas
# Example dataframe
Yup, it crashed my Windows 2000 laptop! In fact, so bad, I had to remove the
battery to restart.
My guess is that this is deep in the heart of the C code in the display
libraries that 2.1.0 uses. Good luck to the R gurus in finding it. In case
it's of any use, here are particulars of my computer:
I can confirm that on my system running Windows 2000 that the
command using R 2.1.0 hangs -- the graphics window appears but
nothing happens subsequently. Under 2.0.0 it works okay.
Patrick Burns
Burns Statistics
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry
Hi all,
As far as I can tell, the only canned way to do an intersect between two
vectors of ints is the intersect(vec1, vec2) function -- is there another
function I'm missing for intersecting more than two vectors??
TIA,
Ken
__
In comparing 2 GAM models are there some guidelines to determine which model
might be a better fit by comparing the 2 UBRE scores? I assume that a
larger value is indicative of a better fit? What magnitude of difference is
a significant one? I imagine that one could also use the percent of
Here's a brute-force approach:
%i% - intersect
x1 - c(1, 3, 4, 7)
x2 - c(3, 7, 8)
x3 - c(7, 1)
x1 %i% x2 %i% x3
[1] 7
Andy
From: Ken Termiso
Hi all,
As far as I can tell, the only canned way to do an intersect
between two
vectors of ints is the intersect(vec1, vec2) function --
I had the same problem (on R2.1.0 on Win XPSP2), let me just add, in case it
is of any help:
1) that even when the plot works in 2.1.0 it takes much longer than in
R2.0.1
2) this is worse if plot(x,type=l) is used (17sec versus 1sec for the
type=l plot with a vector of length 15000)
example
I'm really sorry, I didn't notice the existing labels argument in the
script. I will have another try with the highlighting. This is the complete
code:
setwd(c:/data/marelac/thesis/r/)
read.table(bay1.bay,header=T)-Rdata
names(Rdata)
Rdata[-1] - Rdata
Rdata[1:7] - Rdata
panel.hist - function(x,
No, I know of no canned way to do it, as it's based on match(), which is
for two vectors only.
A brute force loop to do it is simple, as I suspect you realize. Perhaps a
faster way is to use tabulate() to find values that appear n times (for n
vectors). That is:
## assume your integer vectors
Dear All,
I've got into the habit of installing R from the precompiled Debian binaries,
including many of the packages from the r-cran-* Debian packages, and later
building from source (e.g., to link against Goto's BLAS, or to build patched
versions, etc). I install the newly built R to the
Hi,
I have been trying to do numerical differentiation using R.
I found some old S code using Richardson Extrapolation which I managed to get
to work.
I am posting it here in case anyone needs it.
richardson.grad -
Dear R users:
I have a quick question:
is there any way to get a subset of a data frame with the ordered indexes or
without these indexes?
example:
mydata- data.frame(A=seq(1,10), B=c(-5,2,-4,6,-8,-9,2,5,7,0))
subset(mydata, B0)
A B
2 2 2
4 4 6
7 7 2
8 8 5
9 9 7
I would like to obtain
But...
See ?numericDeriv which already does it via a C call and hence is much
faster (and probably more accurate,too).
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the statistician is to catalyze the scientific learning
process. - George E. P. Box
?row.names
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the statistician is to catalyze the scientific learning
process. - George E. P. Box
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of jose silva
Ah... I searched for half an hour for this function... you know, the help
function in R could really be a lot better...
But wait a minute... looking at this, it appears you have to pass in an
expression. What if it is an unknown function, where you only have a handle to
the function, but you
On May 5, 2005, at 6:31 PM, jose silva wrote:
rois there any way to get a subset of a data frame with the ordered
indexes or without these indexes?
example:
mydata- data.frame(A=seq(1,10), B=c(-5,2,-4,6,-8,-9,2,5,7,0))
subset(mydata, B0)
A B
2 2 2
4 4 6
7 7 2
8 8 5
9 9 7
I would like to obtain
Did you try format(x,digits=22). 22 digits is the maximum supported.
I hope this helps
Francisco
From: Huntsinger, Reid [EMAIL PROTECTED]
To: 'Francisco Fuentes' [EMAIL PROTECTED],
r-help@stat.math.ethz.ch
Subject: RE: [R] Precision in R
Date: Thu, 5 May 2005 14:18:50 -0400
The summands
Boxplot.stats seems to be somewhat helpful but not the full answer to my needs
for eliminating outliers. Any other suggestions?
In the first post I mentioned the Appendix A from
http://trb.org/publications/nchrp/nchrp_w71.pdf . They used X and Y varialbes
whereas boxplot.stats is using
Hi,
This is more of a general stat question. I am looking for a easily
computable measure of a distance between two empirical distributions.
Say I have two samples x and y drawn from X and Y. I want to compute a
statistics rho(x,y) which is zero if X = Y and grows as X and Y become
less similar.
Thanks Henric.
If I may I'd like to go a little further ... For example, Johnson and
Wichern's example 7.10:
ex7.10 -
+ data.frame(y1 = c(141.5, 168.9, 154.8, 146.5, 172.8, 160.1, 108.5),
+ y2 = c(301.8, 396.1, 328.2, 307.4, 362.4, 369.5, 229.1),
+ z1 = c(123.5,
Uzuner, Tolga [EMAIL PROTECTED] writes:
Ah... I searched for half an hour for this function... you know, the
help function in R could really be a lot better...
But wait a minute... looking at this, it appears you have to pass in
an expression. What if it is an unknown function, where you
74 matches
Mail list logo