Agnieszka,
Package 'mvpart' is documented. In this case, ?rpart.object explains
*where* in the rpart object is the membership vector.
cheers, jari oksanen
On Mon, 2005-08-08 at 16:02 +0200, [EMAIL PROTECTED] wrote:
> Hi,
>
> I'd like to perform Dufrene-Legendre Indicator Species Analysis for
>
Ok, thanks,
I'll try with a simplier example:
I have a vector with 4 levels
dataframe 1
station temp
aaa12
aaa13
bbb12
bbb20
aaa23
bbb21
ccc30
ccc18
ddd15
aaa11
ddd15
ddd10
and a thresholds vector
You can also use the RODBC package to hold the data in a database, say MySQL
and only import it when you do the modelling, e.g.
> library(RODBC)
> library(sspir)
> con <- odbcConnect("MySQL Test")
> data(vandrivers)
> sqlSave(con,dat=vandrivers,append=FALSE)
> rm(vandrivers)
> gc()
> van.call <-
you could use something like this:
dat1 <- data.frame(station = rep(letters[1:5], 4), temp =
round(rnorm(20, 15, 3)))
dat2 <- data.frame(station = letters[1:5], temp = round(rnorm(5, 15,
4)))
dat <- merge(dat1, dat2, by = "station")
do.call("rbind", lapply(split(dat, dat$station
On 09-Aug-05 Prof Brian Ripley wrote:
> On Mon, 8 Aug 2005, Deming Mi wrote:
>
>> Does anybody know if there is an R function (package) to fit a
>> two-terms exponential model like y = a*exp(bx) + c*exp(dx)
>> where y is dependent variable and x is independent variable.
>> MATLAB has a Curve Fitti
Dear list,
Sorry for answering my own post, but I have had partial sucess in
this. With the panel function below, I get separate histograms in each
panel using the group argument.
histogram(~vot | agem, nint=50,data=work,groups=Type, subset=agem > 24
& agem < 30, panel=panel.grouphist,type="count
I'm wondering if there's a possibility of having a
bunch of models that have already been assembled and
run using R, that might be available to users? My own
interest is in AR and msm models, but examples would
seem to be of great use to all kinds of statistical
analyses with R-- surely "a picture
If 'thr' were a vector with the stations as names,
then you could do (untested):
above <- dataframe1[, 'temp'] > thr[as.character(dataframe1[, 'station'])]
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
ales
Hello Jari,
I've tried to do the analysis this way:
> duleg(data.matrix(dane[,9:28]),predict(fit,type="vector"),numitr=1000)
Error in "row.names<-.data.frame"(`*tmp*`, value = NULL) :
invalid 'row.names' length
I haven o idea what's going on with the row names
cheers :)
Agnieszka
True. But from the "Note" section of help(step) :
This is a minimal implementation. Use 'stepAIC' in package 'MASS'
for a wider range of object classes.
And from author(s) section :
B. D. Ripley: 'step' is a slightly simplified version of 'stepAIC'
in package 'MASS' (Venable
Nice one. But I think you could replace the last line (the one with
do.call) with the simpler
w <- which( dat[ ,2] > dat[ ,3] )
w
[1] 6 11 13 14 16 18 20
dat[ w, ]
station temp.x temp.y
6b 18 16
11 c 17 15
13 d 16 14
14 d 17 1
Hi,
I'm a R newbie and want to accelerate the following pre-filtering
step of a data set with more than 115,000 rows :
#-
# Function to filter out constant data columns
filter.const<-function(X, vectors=c('column', 'row'), tol=0){
realdata=c()
filteredX<-matrix()
if( ve
Le 05.08.2005 22:41, Globe Trotter a écrit :
>Hi,
>
>I have a grouping of some observations. I want to cluster them using
>hierarchical clustering and compare how the hierarchical clustering shows up
>vis-a-vis the groupings. Is it possible to do this in color? I guess what I am
>looking for is a
You are right, but unfortunately this is not the limiting step or
bottleneck in the code below.
The filter.const() function is only used to get the non-constant
columns in the
training data set, which is initially small (49 rows and 525 columns).
And this function is only applied for filtering
I do not fully comprehend the codes below. But if I usually want to
check if all the elements in a row/column are the same, then I would
check the variance or range and see if they are nearly zero.
v.row <- apply( mat, 1, var )
v.col <- apply( mat, 2, var )
tol <- 0
good.row <- which( v.
Dear all,
Could someone please explain the following perculiarity?
> 2 == 0.2/0.1
[1] TRUE
> 3 == 0.3/0.1
[1] FALSE
Similarly,
floor(0.2/0.1) = 2
floor(0.3/0.1) = 2
Thank you,
Simon
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mai
On 8/9/05 7:42 AM, "Simon Woodhead" <[EMAIL PROTECTED]> wrote:
> Dear all,
>
> Could someone please explain the following perculiarity?
>
>> 2 == 0.2/0.1
> [1] TRUE
>> 3 == 0.3/0.1
> [1] FALSE
>
> Similarly,
>
> floor(0.2/0.1) = 2
> floor(0.3/0.1) = 2
This is a FAQ
(http://cran.r-project.or
Simon Woodhead wrote:
> Dear all,
>
> Could someone please explain the following perculiarity?
Please read the FAQ "Why doesn't R think these numbers are equal" (as
the posting gude asks you to do).
Uwe Ligges
> > 2 == 0.2/0.1
> [1] TRUE
> > 3 == 0.3/0.1
> [1] FALSE
>
> Similarly,
>
> f
Look at ?Comparison, especially in the "Note" section, i.e.,
3 == 0.3/0.1
identical(all.equal(3, 0.3/0.1), TRUE)
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +3
Hi,
This is just an FYI documenting a conflict between R and
Google-Desktop. The solution was to click on the Google-Desktop icon
in the systray and click on "Pause Indexing". I also temporarily
suspended anti-virus scanning before successfully
install.packages("VR") many times without getting a
**
This email and any files transmitted with it are confidential and
intended solely for the use of the individual or entity to whom they
are addressed. If you have received this email in error please notify
the system manager.
*
dear R wizards:
plot( 1, 1, ylim=(2,10), xlim=(2,10), type="n");
rect( -1, -1, 12, 12, col=gray(0.99) );
unfortunately wipes out the border axes around the plot. how do I keep
this?
regards,
/ivo
__
R-help@stat.math.ethz.ch mailing list
https://st
Dear all,
Is it possible to calculate the first and second derivatives of a cubic spline
curve (interpolating through the points) as opposed to a smoothing spline?
If so, how is this done? Is it also possible to find the equation of the spline
curve that is fitted?
Any help is much appreciat
On Tue, 9 Aug 2005 [EMAIL PROTECTED] wrote:
>
> dear R wizards:
>
> plot( 1, 1, ylim=(2,10), xlim=(2,10), type="n");
> rect( -1, -1, 12, 12, col=gray(0.99) );
>
> unfortunately wipes out the border axes around the plot. how do I keep
> this?
I think you meant
plot( 1, 1, ylim=c(2,10), xlim=c(
If there is a helpfile open (f.e ?glm) and it is the top window, then
an exception error occurs (closing RGUI)
when I hit the open file button.
If the helpfile is not the top window (of the RGUI) I am able to open a
new script without any error.
The RGUI is not closing complete there is a blank
Knut Krueger schrieb:
>If there you need any further details please ass for.
>
>
ask for :-(
AppName: rgui.exeAppVer: 2.10.50418.0ModName: r.dll
ModVer: 2.10.50418.0 Offset: 77cc
with regards
Knut Krueger
http://www.biostatistic.de
Dear all
I am quite confused by rpart plotting. Here is example.
set.seed(1)
y <- (c(rnorm(10), rnorm(10)+2, rnorm(10)+5))
x <- c(rep(c(1,2,5), c(10,10,10))
fit <- rpart(x~y)
plot(fit)
text(fit)
Text on first split says x < 3.5 and on the second split x < 1.5 what
I understand:
If x < 3.5 so y
Dear list,
when I execute
library(mda)
x1 <- rnorm(10,mean=1,sd=.12)
x2 <- rnorm(10,mean=.6,sd=.12)
y <- pp(x1,1)+pp(x1,1)*pp(x2,.6)+.12*rnorm(length(x1))
mymars <- mars(cbind(x1,x2),y,degree=2)
mymars$cuts
mymars$factor
I get what I expected. Instead, when I execute
x1 <- rnorm(1
Hello,
I'm performing mvpart analysis, I've chosen to interactively pick the tree
(xv="pick") and to do multiple cross-validation (xvmult=10). On the plot x-val
relative error=f(cp/size of the tree) there are some green lines, which didn't
appear during normal cross-validation. What do they dep
I have forgotten to list the definition of pp. Here it is:
pp <- function(x,a) {ifelse(x>a, x-a, 0)}
__
Dear list,
when I execute
library(mda)
x1 <- rnorm(10,mean=1,sd=.12)
x2 <- rnorm(10,mean=.6,sd=.12)
y <- pp(x1,1)+pp(x1,1)*pp(x2,.6)+.12*rnorm(len
With much help from Chuck Cleland, I was able to get xyplot to
generate the plot I wanted. I'm trying to write a script (that can be
read using source("file")) to create the plots I want and save them to
files.
If I type the following lines into the R (in interactive mode), the
correct plot gets
I have a large dataset (about 1 million data points from a
68-dimensional state space, result of an MCMC simulation) which won't
fit in memory. I think that the only solution for analyzing this is
saving it in relational database (when generated) and then reading
back only portions of this data.
Hello all,
X is a list of 20 lists, and each individual list has 65 elements.
Y is a list of 65 elements.
WANT: subtract Y from each of the 20 lists in X.
Here's what I tried and the error messages:
> X - rep(Y, 20)
Error in X - rep(Y 20) : non-numeric argument to binary operator
I tried se
This is in the FAQ: You need to wrap any trellis calls in functions or
scripts in print().
Andy
> From: Jamieson Cobleigh
>
> With much help from Chuck Cleland, I was able to get xyplot to
> generate the plot I wanted. I'm trying to write a script (that can be
> read using source("file")) to c
On 8/9/05 12:21 PM, "Afshartous, David" <[EMAIL PROTECTED]>
wrote:
>
> Hello all,
>
> X is a list of 20 lists, and each individual list has 65 elements.
>
> Y is a list of 65 elements.
>
> WANT: subtract Y from each of the 20 lists in X.
>
> Here's what I tried and the error messages:
>
>>
Suggestions:
1. Read the relevant sections of "An Introduction to R" and the R language
definition to learn how to work with lists.
3. ?lapply to learn how to use the apply family on lists
Hint: sapply(X,"-",y=Y)
The above references will tell you what these hieroglyphics mean and what
you get.
The stat.table function in the Epi package won't do standard deviations.
It didn't seem that it would be difficult to add an "sd" function to the
stat.table function. Following the example for the mean, I set up a
similar function for the sd (and included it as an options) but it just
won't work. (
After you copy stat.table to stat.table2 and modify stat.table2
try:
> environment(stat.table2) <- environment(stat.table)
(you should only need to do that 1 time after creating/editing
stat.table2).
hope this helps,
Greg Snow, Ph.D.
Statistical Data Center, LDS Hospital
Intermountain Health Ca
Like many, I am sure, I get R-Help in digest form. Its easy enough to
browse the
subject lines, but then if an entry interests you, you have to embark
on this tedious search or scroll to find it.
It would be great to have a "clickable" digest, where the topics list
is a set of pointers, and clic
Thanks, that worked. I see now where it says that in the help for
xyplot, however, it wasn't obvious to me before.
Jamie
On 8/9/05, Liaw, Andy <[EMAIL PROTECTED]> wrote:
> This is in the FAQ: You need to wrap any trellis calls in functions or
> scripts in print().
>
> Andy
>
> > From: Jamieso
> Like many, I am sure, I get R-Help in digest form. Its easy enough to
> browse the
> subject lines, but then if an entry interests you, you have to embark
> on this tedious search or scroll to find it.
> It would be great to have a "clickable" digest, where the topics list
> is a set of pointers,
Hi, I am trying to write a package(A) for myself and need to use a
function from another package(B) which is in R already(need to install
it before use). Could anyone tell me how to implement that? Also I
hope that my package gives an ERROR message(something like "STOP,
please install package B fir
Trevor -
There's a wonderful feature in gnus (the emacs information (including
mail) reader), that "bursts" digests into subparts for reading (and
hence, easy access). I think there are other similar tools as well for
other mail readers.
best,
-tony
On 8/9/05, Trevor Hastie <[EMAIL PROTECTED]>
> After you copy stat.table to stat.table2 and modify stat.table2
> try:
>
>> environment(stat.table2) <- environment(stat.table)
>
> (you should only need to do that 1 time after creating/editing
> stat.table2).
>
> hope this helps,
>
> Greg Snow, Ph.D.
> Statistical Data Center, LDS Hospital
> In
See "Writing R Extensions" / Creating R package / Description file /
"depends" field (section 1.1.1).
Jarek
\===
Jarek Tuszynski, PhD. o / \
Science Applications International Corporation <\__,|
(703) 676-419
On Tue, 9 Aug 2005, Peter Yang wrote:
> I am trying to write a package(A) for myself and need to use a
> function from another package(B) which is in R already(need to install
> it before use). Could anyone tell me how to implement that? Also I
> hope that my package gives an ERROR message(somethi
On 09-Aug-05 Trevor Hastie wrote:
> Like many, I am sure, I get R-Help in digest form. Its easy enough
> to browse the subject lines, but then if an entry interests you,
> you have to embark on this tedious search or scroll to find it.
> It would be great to have a "clickable" digest, where the top
A.J. Rossini wrote:
> Trevor -
>
> There's a wonderful feature in gnus (the emacs information (including
> mail) reader), that "bursts" digests into subparts for reading (and
> hence, easy access). I think there are other similar tools as well for
> other mail readers.
If you are using linux you
> Thanks very much for your help. Do I still have to write library( ) to
load the library in my code after I add a "depends" field in the DESCRIPTION
file?
No,
First when user of your package downloads it (at least on Windows machines)
all required packages are also downloaded (this mechanism s
Jeff Gentry <[EMAIL PROTECTED]> writes:
> > Like many, I am sure, I get R-Help in digest form. Its easy enough to
> > browse the
> > subject lines, but then if an entry interests you, you have to embark
> > on this tedious search or scroll to find it.
> > It would be great to have a "clickable" di
> "Trevor" == Trevor Hastie <[EMAIL PROTECTED]>
> on Tue, 9 Aug 2005 10:27:32 -0700 writes:
Trevor> Like many, I am sure, I get R-Help in digest
Trevor> form. Its easy enough to browse the subject lines,
Trevor> but then if an entry interests you, you have to
Trevor> em
[EMAIL PROTECTED] wrote:
> **
> This email and any files transmitted with it are confidential and
> intended solely for the use of the individual or entity to whom they
> are addressed. If you have received this email in error ple
Brief correction: it should read
>> van.call <- call('sqlQuery',con,query='select * from vandrivers;')
rather than
>> van.call <- sqlQuery(con,'select * from vandrivers;')
The latter statement would load the data into memory as usual.
Best wishes,
Andreas
- Original Message -
From
Is there any way to be prompted for a password when using a proxy server
to down load files to linux. The only approach I have seen in the
documentation involves putting a clear text password into an environment
variable, which would get me into trouble.
Thanks for the help,
George Williams
_
I can duplicate this problem under WinXP Professional 5.1.2600 Service Pack
2 Build 2600 on a Toshiba Satellite A10. Thanks to Knut for identifying just
what was causing this problem - I was having this error repeatedly but
couldn't track down the cause.
Duncan
> version
_
Any prospects of getting R-news delivered as an RSS feed???
Duncan
*
Dr. Duncan Mackay
School of Biological Sciences
Flinders University
GPO Box 2100
Adelaide
S.A.5001
AUSTRALIA
Ph (08) 8201 2627FAX (08) 8201 3015
http://www.scieng.flinders.edu.au
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
I can duplicate this problem under WinXP Professional 5.1.2600 Service Pack
2 Build 2600 on a Toshiba Satellite A10. Thanks to Knut for identifying just
what was causing this problem - I was having this error repeatedly but
couldn't track down the cause.
Duncan
> version
_
Any prospects of getting R-news delivered as an RSS feed???
Duncan
*
Dr. Duncan Mackay
School of Biological Sciences
Flinders University
GPO Box 2100
Adelaide
S.A.5001
AUSTRALIA
Ph (08) 8201 2627FAX (08) 8201 3015
http://www.scieng.flinders.edu.au
On 10 August 2005 at 11:14, Duncan Mackay wrote:
| Any prospects of getting R-news delivered as an RSS feed???
The R-help (!!) list is available via Gmane.org, a fabulous service that
archives, redirects, displays, ... a gazillion mailing lists. Among them are
are r-help, r-devel, and several r-s
Hi everyone
The suggestion I received from Mr. Dimitris Rizopoulos and Mr. Gabor
Grothendieck seems to work well. But unfortunately when I am trying to
load nlme package, the following error message is generated. Please
help me on this, at the earliest.
> library(nlme)
Error in loadNamespace(i, c
Could this be a web-proxy problem?
http://cran.r-project.org/bin/windows/base/rw-FAQ.html#The-internet-download-functions-fail_002e
Have you tried doing update.packages() ? does this work?
I presume you're on some sort of Windows?
Sean
On 09/08/05, Frank E Harrell Jr <[EMAIL PROTECTED]> wrote
63 matches
Mail list logo