Hi
most probably testError has different range then testVar, but only
you know if it is true.
What says
str(testError) and str(testVar)
HTH
Petr
On 30 Jan 2006 at 23:01, Michael wrote:
Date sent: Mon, 30 Jan 2006 23:01:04 -0800
From: Michael [EMAIL PROTECTED]
If you're not interested in fitting caribou-specific responses, you
can use beta-binomial logistic models. There are several package
available for this purpose on CRAN, among which aod. Because these
models are fitted using maximum-likelihood methods, you can use AIC
(or other information
Seth Falcon wrote:
On 29 Jan 2006, [EMAIL PROTECTED] wrote:
On Sun, 29 Jan 2006, Elizabeth Purdom wrote:
I came across the following behavior, which seems illogical to me.
What did you expect and why?
I don't know if it is a bug or if I'm missing something:
all(logical(0))
[1] TRUE
Hi
?identical
The safe and reliable way to test two objects for being _exactly_
^
equal. It returns 'TRUE' in this case, 'FALSE' in every other
case.
On Mon, 30 Jan 2006, Patricia J. Hawkins wrote:
ASA == Alexandre Santos Aguiar [EMAIL PROTECTED] writes:
ASA I am new to R and read this list to learn. It is amazing how
ASA frequently new functions pop in messages. Useful and timesaving
ASA functions like subset (above) must be documented
Hello R users
I ran more than 100 logistic regression analyses. Some of the analyses gave
me this kind warning below.
###
Warning messages:
1: algorithm did not converge in: glm.fit(x = X, y = Y, weights = weights,
start = start, etastart
On Mon, 2006-01-30 at 18:53 +0100, Peter Dalgaard wrote:
Gavin Simpson [EMAIL PROTECTED] writes:
On Mon, 2006-01-30 at 08:38 -0800, Waichler, Scott R wrote:
R-2.2.1 is still not available for Redhat Linux as an RPM on CRAN. It
is available as an SRPM. Can someone fill me in on why it
Hi Spencer, yes thanks it does help. In short I have resolved the issue.
In length, not directly, as I have not found a way to ask R to actually plot
frequency vs. period directly, I have scoured the R help pages and Venebles
Ripley but the call to R spectrum() or periodogram() generates error
Hi,
Apology for this question being off the topic (OT) of
R, though I expect
this list might be the best place on the net to ask
this question.
In brief, the question is: what classification
algorithm
can one use if the features are histograms?
I have a classification problem, and believe
Hello,
Is it possible to change the default repositories?
I've already changed the list in the repositories file under etc/
this is how my repositories file under R/etc looks like:
row namemenu_nameURLdefaultsourcewin.binary
mac.binary
CRANCRAN@CRAN@TRUETRUE
Greg Tarpinian [EMAIL PROTECTED] writes:
R2.2 for WinXP, Splus 6.2.1 for HP 9000 Series, HP-UX 11.0.
I am trying to get a handle on why the same lme( ) code gives
such different answers. My output makes me wonder if the
fact that the UNIX box is 64 bits is the reason. The estimated
Nice. May I incorporate these into the collection I am packaging up for
redistribution?
-Eric
-Original Message-
From: Vincent Zoonekynd [mailto:[EMAIL PROTECTED]
Sent: Mon 1/30/2006 6:48 PM
To: R-help@stat.math.ethz.ch
Cc: [EMAIL PROTECTED]; Kort, Eric; [EMAIL PROTECTED]; Thomas
Hi,
just run the Kolmogorov-Smirnov test on R.
Is there any detailed documentation available for the options for the KS
test, esp. with regard to the hypotheses.
The help file is rather thin.
Many thanks.
Bernd
__
R-help@stat.math.ethz.ch mailing
Dear list,
I have some difficulties fine-tuning a lattice conditional histogram plot and
found little help in the documentation.
My dataset consists of plant flowering ages from 3 different altitudes and 3
successional sites at each altitudinal level.
My questions are:
How do I avoid that
What's wrong with
z - spec.pgram(ldeaths, plot = FALSE)
plot(1/z$freq, z$spec, log=y, xlab=period, ylab=power)
?
On Tue, 31 Jan 2006, Tom C Cameron wrote:
Hi Spencer, yes thanks it does help. In short I have resolved the issue.
In length, not directly, as I have not found a way to ask R to
Hi,
I am looking for information how to compute TURF-Analysis and with which
softwares this can be done.
Can you please help me?
Thank you in advance.
Best regards
Afiwa Kuzeawu
Project Manager
76 Prozent der Verbraucher suchen bei Textilien
und Schuhen gezielt
Hi,
very likely your data exhibit quasi-separation which cause (log)Lik to
be monotone and thus ML estimate do not exist. However you can rely on
point estimate and use LRT to test for its significance.
Or Better: have a look to brlr or logistf packages which bypass the
monotone-likelihood
it looks like a barchart rather than an histogram (?).
if it is, you can do something like that:
CT.flow$age - factor(CT.flow$age)
CT.flow$succession -
factor(CT.flow$succession,levels=c(early,mid,late))
ct - expand.grid( succession=levels(CT.flow$succession),
Assa Yeroslaviz wrote:
Hello,
Is it possible to change the default repositories?
I've already changed the list in the repositories file under etc/
this is how my repositories file under R/etc looks like:
row namemenu_nameURLdefaultsourcewin.binary
mac.binary
CRAN
Dear Julien
Have a look at the type argument in ?anova.lme
Regards,
Christoph Buser
--
Christoph Buser [EMAIL PROTECTED]
Seminar fuer Statistik, LEO C13
ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND
phone:
The comment relating R to Stata, SAS and SPSS is now
on the Burns Statistics website at
http://www.burns-stat.com/pages/Tutor/R_relative_statpack.pdf
And also on the UCLA website where the original technical
report can be found as well.
http://www.ats.ucla.edu/stat/technicalreports/
Thank you
Dear Dr. Bates,
Thank you very much for your response.
--- Douglas Bates [EMAIL PROTECTED] escribió:
On 1/27/06, gabriela escati peñaloza
[EMAIL PROTECTED] wrote:
Hi, I' m having a hard time understanding the
computation of degrees of freedom
So do I and I'm one of the authors of the
Hello,
package arules is very well suited for mining association rules,
interfacing implementations of apriori and eclat.
What is not realised is an implementation of the Frequent Pattern Tree (Han,
J., Pei, J. and Yiwen, Y. (2000).
There are several free implementations available
Dear all
I want to have vector of length 24 which is filled by zeroes except
that on position
c(8,16,19,20) is 2
and
c(3,4,12,24) is 4
instead.
I am able to do
(1:24)%in%c(8,16,19,20)*2
[1] 0 0 0 0 0 0 0 2 0 0 0 0 0 0 0 2 0 0 2 2 0 0 0 0
and
(1:24)%in%c(3,4,12,24)*4
[1] 0 0 4 4 0 0 0
consider the following:
x - numeric(24)
x[c(8,16,19,20)] - 2
x[c(3,4,12,24)] - 4
x
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel:
Patricia,
If I understand correctly what you need, this is already available through
the main help page. If you type help.start() to get to the html help, then
follow the link Packages. This will give you an index of all available
packages installed on your system (base + contributed).
Hi,
I would like to perform computations on some variables belonging to the same
dataframe. For instance my data frame has the following shape:
toto1 toto2 toto3 toto4 toto5
1 1 2 3 4 5
2 6 7 8 9 10
I would like to
In addition, the search page at
http://finzi.psych.upenn.edu
can search all functions of all CRAN packages.
This is also available through
RSiteSearch(string,restrict=functions).
See the help page for RSiteSearch.
Jon
--
Jonathan Baron, Professor of Psychology, University of Pennsylvania
Received some excellent responses to my query. The problem was with the
mysql-devel libraries.
I upgraded/installed the mysql client/server/devel libraries to 4.1.16 and
issued the well known fixes of :
export PKG_LIBS='-L/usr/lib64/mysql -lmysqlclient'
export
Your post seems to be messed up but I will assume you have a
5 column data frame and the questino is how to run f on the first
three columns and separately on the last two. I think the
easiest is just the following where I have used the builtin iris
data set where I have assumed that the
Thanks to all who responded. I definitly need some fresh air to clean
my head. :-)
Petr
On 31 Jan 2006 at 14:56, Petr Pikal wrote:
From: Petr Pikal [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Date sent: Tue, 31 Jan 2006 14:56:47 +0100
Hi,
I am using the xyplot function in the lattice package to generate
multiple plots, but I would like to have them plotted on the same page. I
would like to set something equivalent to the command: par(mfrow=c(2,2)),
in order that I can plot 4 xyplots on the same page. How can I do this in
[EMAIL PROTECTED] wrote:
I am creating habitat selection models for caribou and other species with
data collected from GPS collars. In my current situation the
radio-collars
recorded the locations of 30 caribou every 6 hours. I am then comparing
resources used at caribou locations to
On Tuesday 31 January 2006 06:55, Gabor Grothendieck wrote:
On 1/30/06, Patricia J. Hawkins [EMAIL PROTECTED] wrote:
[...snip...]
#which generalizes to:
bb - matrix(1:50, ncol=10, nrow=5, byrow=TRUE)
bv - as.vector(bb)
ai - as.vector(aa) + rep((1:nrow(aa)-1)*10, each=3)
bv[ai] -
Dear R-users
Thanks for all the help/suggestions, they have been most helpful.
Best regards
Pryseley
-
[[alternative HTML version deleted]]
__
sudoku_2.0 is now available on CRAN and mirrors.
Thanks to some terrific contributions and suggestions, especially by
new co-author Greg Snow, this version has a full GUI interface (native
on Windows, or using tkrplot), a puzzle generator, and a function to
fetch recently published puzzles. The
B Dittmann herrdittmann at yahoo.co.uk writes:
Hi,
just run the Kolmogorov-Smirnov test on R.
Is there any detailed documentation available for the options for the KS
test, esp. with regard to the hypotheses.
The help file is rather thin.
Many thanks.
Bernd
Hmmm. Does
Hello,
I'm actually using the Step procedure in R for multiple regression analysis.
I'm using the stepwise selection which alternates between forward selection
and backward elimination (direction both in the step procedure).
I would like to know which F-levels R is using to enter and then to
On 1/29/06, Søren Højsgaard [EMAIL PROTECTED] wrote:
In connection with calculating Monte Carlo p-values based on sampled data
sets: The calculations involve something like
update(lmer.model, data=newdata)
where newdata is a simulated dataset comming from simulate(lmer.model). I
guess
If this is step() (not Step), it does not use tests. Please consult the
help page which says
step package:stats R Documentation
Choose a model by AIC in a Stepwise Algorithm
Description:
Select a formula-based model by AIC.
...
It seems these days we
See ?print.trellis.
Andy
From: Ravi Varadhan
Hi,
I am using the xyplot function in the lattice package to generate
multiple plots, but I would like to have them plotted on the
same page. I
would like to set something equivalent to the command:
par(mfrow=c(2,2)),
in order that
Dear R-list,
maybe someone can help me with the following mixed-effects models
problem, as I am unable to figure it out with the 'nlme-bible'.
I would like to fit (in R, obviously) a so-called animal model (google
e. g. Heritability and genetic constraints of life-history by Pettay
et al.) to
Hi everyone,
I have performed the Kolmogorov-Smirnov test (ks.test) with R the first
time.
What I am not sure about is the exact alternative hypotheses (H1) given
any H0.
According to Conover (1971) Practical Nonparametric Statistics, chapter
6, the following one-sample tests can be
Hi All,
I have a sample x={x1,x2,..,xn} fom a distribution with density f. I
wish to estimate the density. I know a priori that the density is
monotonically decreasing. Is there a way to do this in R?
Thanks
Debayan
__
R-help@stat.math.ethz.ch
I'm running SVM from e1071 package on a data with ~150 columns (variables)
and 5 lines of data (it takes a bit of time) for radial kernel for
different gamma and cost values.
I get a very large models with at least
3 vectors and the prediction I get is not the best one. What does it
On 30 Jan 2006, at 22:01, David Reitter wrote:
I'm trying to use predict with a linear mixed-effects logistic
regression model fitted with nlmmPQL from the MASS library.
Unfortunately, I'm getting an error non-conformable arguments in
predict.lme, and I would like to understand why.
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to
announce 2-day R/S-plus Fundamentals and Programming
Techniques in San Francisco: www.xlsolutions-corp.com/Rfund.htm
San Francisco, February 16-17
Seattle,February 20-21
Boston,February
I have done a fair bit of spectral analysis, and hadn't finished collecting my
thoughts for a reply, so hadn't replied yet.
What exactly do you mean by normalize?
I have not used the functons periodogram or spectrum, however from the
description for periodogram it appears that it returns the
a:b means - all the element in the vector from a to b
a[,-1] means for the matrix a keep all the rows but not the last or
first -can't remember column.
When in doubt do what I do make a small matrix and apply the command see
what it does.
After this:
a=matrix(c(1:9),3,3)
a[,-1]
[,1]
This post has some details:
http://tolstoy.newcastle.edu.au/~rking/R/help/04/10/5581.html
On 1/31/06, Leif Kirschenbaum [EMAIL PROTECTED] wrote:
I have done a fair bit of spectral analysis, and hadn't finished collecting
my thoughts for a reply, so hadn't replied yet.
What exactly do you
: in a formula is not the same as : otherwise!
Andy
From: Ionut Florescu
a:b means - all the element in the vector from a to b
a[,-1] means for the matrix a keep all the rows but not the last or
first -can't remember column.
When in doubt do what I do make a small matrix and apply the
I am using pnorm() with the log.p=T argument to get approximations to ln
\Phi(x) and qnorm with the log.p=T argument to get estimates of
\Phi^{-1}(exp(x)). What approximations are used in these two functions (I
noticed in the source pnorm.c it doesn't look like Abramowitz and Stegen) and
where
Dear List:
I'm trying to use the boot function to estimate some standard errors. I
actually programmed a bootstrap using some homebrew code and it worked
fine. But, I am trying to use the more efficient boot function. I have
placed some sample data for replication of my problem at the bottom of
Hi
I have a question concerning how to match word boundaries which I bet has a
very simple answer, but I haven't found it with trial and error nor by
searching the help archives for the terms in the subject line. The problem is
this: I have a vector of two character strings.
text-c(This is a
When I tried it on Windows XP there was a grinding sound, probably memory
being swapped and it just seemed to go on forever and I finally had to kill R.
I am using R version 2.2.1, 2005-12-20. What did seem to work was this:
gregexpr(X, gsub(\\b\\w|\\w\\b, X, text))
where X should be replaced
Suppose I have a formula
mfTest = (y ~ a + b + c + d1 + d2 + d3);
When I update the formula
mfTest2nd = update(mfTest, .~.^2);
I get all the combinations a:b, a:c, etc. But I don't want the d1:d2, d2:d3,
etc. interactions. I've tried
mfTest2nd = update(fmTest, .~.^2
Hi R users
I have a simple data for calculating goodness-of-fit statistics (e.g., X2 by
Pearson, G2 by Wilks)
#
observed-c(424,174,0,402)
expected-c(282.7174, 314.2972, 142.3142, 260.6712)
2*sum(observed*log(observed/expected)) # for X2
The trouble is log(0/anything) = log(0) = NaN. If you want them to evaulate to
zero you might try zeroing out the values you know will be NaN:
X2 = 2*sum(observed*log(observed/expected));
X2[observed==0] = 0;
--Brett
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL
Hi,
I have a data set with a continuous predictor X, a factor A and a continuous
dependent
variable Y.
I am trying to build a linear model of the form:
Y = (b0 + b1*X1)*B(A)
where B(A) is a constant for each level of the factor A.
I am not quite sure how to formulate the appropriate model
I have question (curiosity) regarding returned values of R's qcauchy
() function,
for nonexceedance probability (F). It seems the ideal returned range
of cauchy distribution should be [-Inf,Inf].
For F=0
qcauchy(0)
[1] -Inf
but for F=1
qcauchy(1)
[1] 8.16562e+15
It seems to me that the
I would like to open a console application from within R and then send input
to the newly opened console. I can use the following to open the application
...
setwd(path)
system(paste(path,icl.exe,sep=))
This allows me to type in what want, but I would like to do something like
this
setwd(path)
Don't know how you tell if its working but you need to write
intern = TRUE
in order to get the output back. For example, this works for me:
system(findstr [a-m], input = letters, intern = TRUE)
On 1/31/06, Jonathan Weeks [EMAIL PROTECTED] wrote:
I would like to open a console
On Tue, 31 Jan 2006, Morey, Richard D (UMC-Student) wrote:
I am using pnorm() with the log.p=T argument to get approximations to ln
\Phi(x) and qnorm with the log.p=T argument to get estimates of
\Phi^{-1}(exp(x)). What approximations are used in these two functions
(I noticed in the
Georges Orlowski wrote:
I'm running SVM from e1071 package on a data with ~150 columns (variables)
and 5 lines of data (it takes a bit of time) for radial kernel for
different gamma and cost values.
I get a very large models with at least
3 vectors and the prediction I get is
64 matches
Mail list logo