On Mon, 8 Aug 2005, Haibo Huang wrote:
I had a logit regression, but don't really know how to
handle the Warning message: non-integer #successes in
a binomial glm! in: eval(expr, envir, enclos)
problem. I had the same logit regression without
weights and it worked out without the warning,
First, AIC is `An Information Criterion' according to Akaike.
Second, help.search(AIC) turns up
stepAIC(MASS) Choose a model by AIC in a Stepwise Algorithm
multinom(nnet) Fit Multinomial Log-linear Models
mosaicplot(graphics)Mosaic Plots
AIC(stats) Akaike's An
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 Fitting Toolbox
to implement this fitting, but I
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
a
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
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 Fitting
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,
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)
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' (Venables
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 14
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(
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 way
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.row
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
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
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,
floor(0.2/0.1)
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:
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 an
**
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
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
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(2,10),
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 is
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 -
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
I have forgotten to list the definition of pp. Here it is:
pp - function(x,a) {ifelse(xa, 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 -
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
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 create
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:
X - rep(Y, 20)
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. (I
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
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
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: Jamieson
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
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
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
Intermountain
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-4192
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(something
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 topics
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 can use
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
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 digest, where
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 embark on this
[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 please
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
__
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
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
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,
65 matches
Mail list logo
