On 1/9/07, rhk [EMAIL PROTECTED] wrote:
Hello, I am afraid I do not fully understand all intricacies of
programming in lattice plots. In the code below I try to plot an ROC
curve, following R-news 4(1). When I condition on the variable 'group' I
get the error message below, when I plot the
I did just the download of the pls package, but the NIR dataset is not
available
require(pls)
[1] TRUE
data(NIR)
Warning message:
data set 'NIR' not found in: data(NIR)
is there another package with the dataset for the examples?
With regards Carmen
NIR seems to be not dataset. It is a member of dataset yarn.
Try:
library(pls)
data(yarn)
str(yarn)
On 1/10/07, Carmen Meier [EMAIL PROTECTED] wrote:
I did just the download of the pls package, but the NIR dataset is not
available
require(pls)
[1] TRUE
data(NIR)
Warning message:
data
See changelog of pls 2.0-0
- The 'NIR' and 'sensory' data sets have been renamed to 'yarn' and 'oliveoil'.
you can see it in package source from CRAN
http://cran.r-project.org/src/contrib/Descriptions/pls.html
HTH
On 1/10/07, Carmen Meier [EMAIL PROTECTED] wrote:
talepanda schrieb:
NIR
On Tue, 9 Jan 2007, Bálint Czúcz wrote:
There is an improved version of the original random forest algorithm
available in the party package (you can find some additional
information on the details here:
http://www.stat.uni-muenchen.de/sfb386/papers/dsp/paper490.pdf ).
I do not know whether it
On Tue, 9 Jan 2007, Walter Durka wrote:
dear all,
Walter,
_please_ cc questions on contributed packages to the maintainer, since not
everybody follows r-help closely...
I want to perform a posthoc test for my ANCOVA:
a1-aov(seeds~treatment*length)
With
summary(glht(a1, linfct =
Dear WizaRds!
A trivial question indeed on selecting subsets in data frames. I am
sorry. Unfortunately, I did not find any helpful information on the
introduction, searched the help archive and read in introductory books.
Please help:
I want to select column KB which is read via read.csv2 as
Mark Hempelmann wrote:
Dear WizaRds!
A trivial question indeed on selecting subsets in data frames. I am
sorry. Unfortunately, I did not find any helpful information on the
introduction, searched the help archive and read in introductory books.
Please help:
I want to select column KB
Generally it is difficult to get an overview of what's there.
But the following function I acquired from (???) ages ago does a nice
job:
lls -
function (pos = 1, pat = )
{
dimx - function(dd) if (is.null(dim(dd)))
length(dd)
else dim(dd)
lll - ls(pos = pos, pat = pat)
John Lawrence Aspden wrote:
I've got a library (brainwaver), installed locally in ~/R/library, and
this information is recorded in the ~/.Renviron file.
In my script I load the library, but if I call it using
#!/usr/bin/r --vanilla, this stops working.
(Various private e-mails exchanged.
Dear R experts,
I have been encountering problems with the optim routine using for
loops. I am determining the optimal parameters of several nested models by
minimizing the negative Log-Likelihood (NLL) of a dataset.
The aim is to find the model which best describes the data. To this end,
Hi
I listened to all your advise and ran my data on a computer with a 64 bits
procesor but i still get the same error saying it cannot allocate a vector of
that size 1240 kb . I don't want to cut my data in smaller pieces because we
are looking at interaction. So are there any other options
I want to write something like (LaTeX style)
b_{norm}=\frac{F\prime(0)}{R\prime(0)}
how do I add the prime (first derivative) to a R-plot? The help of
plotmath just talks about partialdiff. Can you complete this
command?
text( 30,0.05,labels=expression(b[plain(norm)]==frac(F(0),R(0))) )
Thanks,
try:
d[d[,KB] == Eenie,]
or
d[d$KB==Eenie,]
What you were comparing is just two character strings (KB==Eenie) which
of course is FALSE and won't select anything. You have to explicitly
compare values from the KB column of your data
On 1/10/07, Mark Hempelmann [EMAIL PROTECTED] wrote:
Dear
how do I add the prime (first derivative) to a R-plot?
sorry for the noise, I found it myself:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/20984.html
I use now (works fine)
plot(1,1,xlab=expression(frac(F'(0),R'(0))),xaxt=n)
Thomas
__
Dear R-help community,
I have 2 different arrays of precipitation data each of the same dimensions of
[longitude, latitude, time] dim=[30,32,43], called array1 and array2. I need to
correlate them. This is the code I used to get one overall correlation value
for
the whole of the area of
Dear All;
I have used fbmSim to simulate a fbm sequence, however, when I tried to
estimate the Hurst effect, none of the nine procedures gave me an answer
close enough to the real value, which is 0.5 (n=1000). So, would you please
advice,
1. which is the best method to estimate the H among the 9
Hi,
I am downloading Bloomberg data from R. This data will be stored in a
zoo object by default. The command is
dat-blpGetData(con,c(NOK1V FH Equity,AUA AV
Equity),PX_OPEN,start=as.chron(as.Date(12/1/2006,
%m/%d/%Y)),end=as.chron(as.Date(12/28/2006, %m/%d/%Y)))
Here I am downloading the
[By now way off the subject line. Something like 'how to set the
libraries inside an R session'.]
On Wed, 10 Jan 2007, John Lawrence Aspden wrote:
John Lawrence Aspden wrote:
I've got a library (brainwaver), installed locally in ~/R/library, and
this information is recorded in the
I think the problem is in the BlpGetData function than in the zoo
object. Because instead of zoo object, the data was put into a data
frame, but then also the problem prevails.
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Shubha Vishwanath
Karanth
You have applied the Hurst parameter estimating function to the fBm, not
the fGn as described in the help page for Long Range Dependence
Modelling in fSeries.
The function aggvarFit computes the Hurst exponent from the variance of
an aggregated FGN or FARIMA time series process
Trying
Two possibilities are ' (which is the usual way to write \prime in TeX)
and minute (which is the nearest character in the Adobe symbol font that
plotmath uses). TeX has its own symbol fonts, and most R devices are
based rather on the Adobe one.
On Wed, 10 Jan 2007, Thomas Steiner wrote:
I
On Wed, 10 Jan 2007, Iris Kolder wrote:
Hi
I listened to all your advise and ran my data on a computer with a 64
bits procesor but i still get the same error saying it cannot allocate
a vector of that size 1240 kb . I don't want to cut my data in smaller
pieces because we are looking at
Date: Tue, 9 Jan 2007 11:13:41 + (GMT)
From: Mark Difford [EMAIL PROTECTED]
Subject: Re: [R] contingency table analysis; generalized linear model
Dear List,
I would appreciate help on the following matter:
I am aware that higher dimensional contingency tables can be
analysed using
Greetings -
Does anybody know if there is an R implementation of the Wei-Lachin
tests for incomplete multivariate observations (JASA 79: 653-661)?
The authors' example of data where this test can be applied is in a
study comparing treatment to placebo where the outcome is a series of
Hello,
I am fairly new to time series analysis, and I am wondering if anyone could
point me to some good references related to mechanics of what is going on
when you fit a model using arima that includes a set of external regressors
(i.e., use the xreg= feature).
Thank you,
Spencer
Hi All:
I'm testing a set of data classification algorithms in this paper
(www.stat.wisc.edu/~loh/treeprogs/quest1.7/mach1317.pdf )
I couldn't find such algorithms in R packages:
1. LOG: polytomous logistic regression (there was one in MASS
library: multinom. But after I
Some computer languages, including C, have a go or go to command
which can be used to shift control to a different part of a function.
Question: Does the R language have a corresponding command? (Yes, I am
aware that it can be abused; but, in the hands of a good programmer,
it is simple to
Hi Jenny!
So if i understand your datafile corect you have 960 case for a year. Any
you have 43 years.. Yes?
I'm not sure you should use correlation in this situation because of the
autocorrelation of the data. There are big autocorrelation on spatial data's
like what you use, and there are also
Dear all,
I try to first link data in a data.frame() to a (polygon) read.shape()
object and then to plot a polygon map showing the data in the
data.frame. The first linking is called link in ArcView and relate
in ArcMap. I use the code shown below, though without success.
Help with this would
Hi Zoltan,
Right, I have 30x32=960 data points per year (It is actually the mean febuary
precipitation total in case you were wondering) at each grid point over the
world, so I have 960 data points each of the 43 years. Therefore can I do
anything with a trend and residuals? I don't think I
Without more detail - a reproducible example - it is hard to give you
concrete advice.
I wonder if the functions NLL23 and NLL21 depend on numerical solutions
of a system of ODEs, since you invoke the odesolve package? If so, try
switching to the Nelder-Mead optimizer, enforcing the parameter
Looks like oneway.test has been changed for R 2.5.0.
Paste the code in this file:
https://svn.r-project.org/R/trunk/src/library/stats/R/oneway.test.R
into your session. Then fun.2 from your post will work without
the workarounds I posted:
fun.2(values ~ group)
On 1/9/07, Adrian Dusa
On 1/10/2007 12:13 AM, Thomas L Jones wrote:
Some computer languages, including C, have a go or go to command
which can be used to shift control to a different part of a function.
Question: Does the R language have a corresponding command? (Yes, I am
aware that it can be abused; but, in
1. multinom is is the nnet package
2. There is a polyclass function in package polspline
On 10/01/07, Feng Qiu [EMAIL PROTECTED] wrote:
Hi All:
I'm testing a set of data classification algorithms in this paper
(www.stat.wisc.edu/~loh/treeprogs/quest1.7/mach1317.pdf )
I couldn't
Greetings,
I am new to R and have been comparing CPH survival analysis hazard ratios
between R and SAS PhReg. The binary covariates' HRs are the same, however
the continuous variables, for example age, have quite different HRs
although in the same direction. SAS PhReg produces HRs which are
On Wednesday 10 January 2007 19:03, Gabor Grothendieck wrote:
Looks like oneway.test has been changed for R 2.5.0.
Paste the code in this file:
https://svn.r-project.org/R/trunk/src/library/stats/R/oneway.test.R
into your session. Then fun.2 from your post will work without
the
Brian Ripley wrote:
Exactly as documented. The argument is named 'new' and not 'add', BTW.
Please 'be careful' in what you say about the work of others.
Agreed, no criticism intended. I really like R. Sorry.
Cheers, John.
--
Contractor in Cambridge UK -- http://www.aspden.com
Dear all,
When using latex.table from the quantreg package, I don't seem to be able to set
table.env=FALSE: when I don't specify caption (as I think I should, when
understanding the R help rightly(?)), I get an error message, and when I
do so, of course I get one, as well.
The funny thing is,
can anybody provide a reference to a document on using DBI and what is
needed to make it work?
I try: drv=dbDriver(Oracle)
and it complains about not finding a function.
The DBI.pdf does not require anything else installed, but looking at the
archive I discovered there is more stuff to be
The usual R-help etiquette recommends:
1. questions about packages go to the maintainer, not to R-help.
2. examples should be reproducible: ie self contained.
if you look carefully at the function latex.summary.rqs you will see
that there is a failure to pass the argument ... on to
*An Introduction to Data Analysis Using R*
by Paul Mathews
for The Quality Managers Network, Lakeland Community College, Kirtland, Ohio
26 January 2007, 7:30-9:00 AM, Room T136.
R is free software used for graphical and statistical data analysis that
can be downloaded from the Internet. Because
The way MySQL works, I use RMySQL to contact, which in turn uses DBI.
There is a library ROracle which, from the manual, works the same. Hence I
would start by looking at Roracle for the connection and proceed from there.
Kees
On Wednesday 10 January 2007 19:39, Bond, Stephen wrote:
can
On Wed, 10 Jan 2007, Tord Snäll wrote:
Dear all,
I try to first link data in a data.frame() to a (polygon) read.shape()
object and then to plot a polygon map showing the data in the
data.frame. The first linking is called link in ArcView and relate
in ArcMap. I use the code shown below,
I am trying to install mixtools on Debian Etch and get the following error
dell-2 /usr/lib # R CMD INSTALL mixtools_0.1.0.tar.gz
* Installing *source* package 'mixtools' ...
** libs
gcc -I/usr/share/R/include -I/usr/share/R/include -fpic -g -O2 -std=gnu99
-c new_svalues.c -o new_svalues.o
Hi David:
Thanks for you information. 2 further questions:
1. I found out that multinom is not doing politomous logistic
regression, do you know which function does this?
2. the polyclass in polspline does polychotomous regression,
while I'm looking for
Thomas L Jones wrote:
Some computer languages, including C, have a go or go to command
which can be used to shift control to a different part of a function.
Question: Does the R language have a corresponding command? (Yes, I am
aware that it can be abused; but, in the hands of a good
[EMAIL PROTECTED] wrote:
Greetings,
I am new to R and have been comparing CPH survival analysis hazard ratios
between R and SAS PhReg. The binary covariates' HRs are the same, however
the continuous variables, for example age, have quite different HRs
although in the same direction. SAS
On Wed, 10 Jan 2007, [EMAIL PROTECTED] wrote:
I am new to R and have been comparing CPH survival analysis hazard ratios
between R and SAS PhReg. The binary covariates' HRs are the same, however
the continuous variables, for example age, have quite different HRs
although in the same direction.
I am hoping something has changed since I last asked about this.
Is there any good source of documentation, including examples, of using
tcl/tk as a gui for R applications?
__
R-help@stat.math.ethz.ch mailing list
For example, this works:
x = seq(-100, 1000, 25)
y = x * x
plot(x,y, xaxt=n)
axis(1,x,FALSE,tcl=-0.3)
axis(1,x[x %% 100 ==0])
It creates two axis objects and the values of the x-axis are the
labels. The following scenario is more difficult, because it uses
'image' to plot a grid of values:
a =
Two possibilities for why your 7 parameter model fits worse than the 6
are that you are finding a local maximum, which might suggest using a
different parameterisation or the functions are accessing some global
data and so aren't behaving as expected. This could be why they work
properly when run
Michael Prager wrote:
I am hoping something has changed since I last asked about this.
Is there any good source of documentation, including examples, of
using tcl/tk as a gui for R applications?
Well, this probably hasn't changed, but did you see it in the first place?
[Philippe Grosjean]
Please, don't reinvent the wheel: putting functions in a dedicated
environment is one of the things done by R packages (together with a
good documentation of the function, and making them easily installable
on any R implementation).
[...] this is probably the time for you
Hello list,
plotting the following example 1 in lattice only labels the x-axis with one
tick mark at the beginning of February.
In contrast the standard plot() function labels the x-axis with equally
spaced tick marks.
#example 1 ---
library(lattice)
date - as.Date(c(1/10/2006,
The 'Right Thing' is for oneway.test() to allow a variable for the first
argument, and I have altered it in R-patched and R-devel to do so. So if
your students can make use of R-patched that would be the best solution.
If not, perhaps you could make a copy of oneway.test from R-patched
## The WoodEnergy example in package HH (available on CRAN) is similar
## to the example you asked about. I have two factors and a
## covariate, and the interaction between the factors and covariate is
## significant. I construct the multiple comparisons of factor Stove
## at each level of
Hi,
I'd like to produce a marginless or zero margin plot so that the pixel
coordinates represent the mathematics.
xy-data.frame(x=c(0,1,1,0,0),y=c(0,1,0,0,1))
png('junk.png',width=300,height=300)
par(mar=c(0,0,0,0))
plot(xy$x,xy$y,xlim=c(0,1),ylim=c(,1))
dev.off()
The resultant file has about
On Wed, 2007-01-10 at 21:18 -0600, David Forrest wrote:
Hi,
I'd like to produce a marginless or zero margin plot so that the pixel
coordinates represent the mathematics.
xy-data.frame(x=c(0,1,1,0,0),y=c(0,1,0,0,1))
png('junk.png',width=300,height=300)
par(mar=c(0,0,0,0))
59 matches
Mail list logo