Uwe Ligges ligges at statistik.uni-dortmund.de writes:
Well, in principle you cannot easily (it's a warning, not an errror, BTW).
But you can fake a little bit by, e.g., renaming z to pch:
panel.hist - function(x, ...)
{
usr - par(usr); on.exit(par(usr))
On Sat, 4 Dec 2004 18:37:43 -0700, Nathan Leon Pace, MD, MStat
[EMAIL PROTECTED] wrote:
Hi,
I using the boot package 1.2-20 on R 2.0.1.
My statistics function estimates 6 parameters.
In a small percentage of resampled data sets my statistics function
doesn't produce an estimate for one
Wensui Liu wrote:
Is it possible use SAS/DDE to communicate with R?
I don't know of any DDE implementation for R.
Uwe Ligges
Sorry for offending useR. ^_^
Thanks a lot.
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Rolf Turner rolf at math.unb.ca writes:
:
: I gather from reading the back-issues of r-help that it should be
: possible (modulo a number of caveats) to read an excel (yuck!) file
: into R using RODBC. I have obtained and installed ODBC and the RODBC
: package, but cannot for the life of me
On 05-Dec-04 Patrick Foley wrote:
It is easy to spot response nonlinearity in normal linear
models using plot(something.lm).
However plot(something.glm) produces artifactual peculiarities
since the diagnostic residuals are constrained by the fact
that y can only take values 0 or 1.
What do
(Ted Harding) wrote:
On 05-Dec-04 Patrick Foley wrote:
It is easy to spot response nonlinearity in normal linear
models using plot(something.lm).
However plot(something.glm) produces artifactual peculiarities
since the diagnostic residuals are constrained by the fact
that y can only take values 0
On 05-Dec-04 Ted Harding wrote:
[...]
For example, adopting the ritual sigificant == P0.05,
power = 80%, you can see a histogram of the p-values
over the conventional significance breaks with
hist(pvals,breaks=c(0,0.01,0.03,0.1,0.5,0.9,0.95,0.99,1),freq=TRUE)
Sorry for the typo! That
On Sun, 2004-12-05 at 01:37, Nathan Leon Pace, MD, MStat wrote:
Hi,
I using the boot package 1.2-20 on R 2.0.1.
My statistics function estimates 6 parameters.
In a small percentage of resampled data sets my statistics function
doesn't produce an estimate for one parameter and the boot
Peter Dalgaard [EMAIL PROTECTED] writes:
Re. the smoothed residuals, you do need to be careful about the
smoother. Some of the robust ones will do precisely the wrong thing
in this context: You really are interested in the mean, not some
trimmed mean (which can easily amount to throwing away
Hello,
Does anyone know of a slick way to get R to convert a string which is all
upper case to a string where the first letter in each word is upper case and
all others are lower case?
I suspect the solution is to begin by parsing the string, convert the
appropriate letters to upper and
Damian Betebenner a écrit :
Hello,
Does anyone know of a slick way to get R to convert a string which
is all upper case to a string where the first letter in each word is
upper case and all others are lower case?
I suspect the solution is to begin by parsing the string, convert the
appropriate
Damian Betebenner damian.betebenner at bc.edu writes:
:
: Hello,
:
: Does anyone know of a slick way to get R to convert a string which is all
upper case to a string where the first
: letter in each word is upper case and all others are lower case?
:
: I suspect the solution is to begin by
a data set like this:
year monthdaycount
2001 11 10
2001 12 11
2004 7178
basically it is a count of of some numbers everyday through a few years
now I'd like to get the mean of the count for every day.
Renaud Lancelot [EMAIL PROTECTED] writes:
This question was posted a while ago and answered by Christian
Hoffmann. I wrapped it in a small function:
CapLeading - function (string){
...
Version of same:
s - the quick red fox jumps over the lazy brown dog
ss - strsplit(s, )[[1]]
ss1 -
Terry Mu muster at gmail.com writes:
:
: a data set like this:
:
: year monthdaycount
: 2001 11 10
: 2001 12 11
:
: 2004 7178
:
:
: basically it is a count of of some numbers everyday through a few
On Sun, Dec 05, 2004 at 12:58:19PM +, Gabor Grothendieck wrote:
Rolf Turner rolf at math.unb.ca writes:
:
: I gather from reading the back-issues of r-help that it should be
: possible (modulo a number of caveats) to read an excel (yuck!) file
: into R using RODBC. I have obtained and
thank you,
one more question:
How can I list only values for given condition?
for example, I want to see only data from febuary, i.e. month=2?
On Sun, 5 Dec 2004 16:26:11 + (UTC), Gabor Grothendieck
[EMAIL PROTECTED] wrote:
Terry Mu muster at gmail.com writes:
:
: a data set like this:
Sorry, I didn't express my problem clearly. Your answers surely help,
but it's not what I meant to ask.
I'd like mean for every day in a year, i.e,
mean for January 1st, 2nd, 3rd,
Febuary, 1st, ..., 28th, 29th(if leap year)
...
There is only one count per day, so there is no need to
Do you mean the mean for all January firsts, the mean for all January 2nds,
etc.?
Look up ?format in addition to the commands you looked up before:
tapply(DF$count, format(Dates, %Y-%m), mean)
Terry Mu muster at gmail.com writes:
:
: Sorry, I didn't express my problem clearly. Your
Place Dates in DF and use subset(DF, month == 2) in place of DF.
I think you need to reread the Introduction to R manual since you
should find this sort of question answered there. If not the pointers
I gave you should get you started.
Terry Mu muster at gmail.com writes:
:
: thank you,
:
:
In my previous post the format should have been %m %d.
Terry Mu muster at gmail.com writes:
:
: thank you,
:
: one more question:
:
: How can I list only values for given condition?
: for example, I want to see only data from febuary, i.e. month=2?
:
: On Sun, 5 Dec 2004 16:26:11 +
thanks, this time it works perfectly.
I use
tapply(DF$count, format(Dates, '%m-%d'), mean).
as Patrick mentioned. It works without as.date().
with as.date() on, it gives a result one day after the actual date.
I also got subset().
Thank you all very much. I like R and will learn more.
Terry
Dear all,
Thanks to the help from the R mailing list we now have a script performing
the required tasks, i.e. applying a mask and filter to a 3D array. See
script below. Any further suggestions to simplify the code are very
welcome. We are currently debugging the larger script of which this
: Date: Sun, 5 Dec 2004 10:35:46 -0600
: From: Dirk Eddelbuettel [EMAIL PROTECTED]
: To: Gabor Grothendieck [EMAIL PROTECTED]
: Cc: [EMAIL PROTECTED]
: Subject: Re: [R] Excel *.xls files, RODBC
:
:
: On Sun, Dec 05, 2004 at 12:58:19PM +, Gabor Grothendieck wrote:
: Rolf
On 05-Dec-04 Peter Dalgaard wrote:
Peter Dalgaard [EMAIL PROTECTED] writes:
Re. the smoothed residuals, you do need to be careful about the
smoother. Some of the robust ones will do precisely the wrong thing
in this context: You really are interested in the mean, not some
trimmed mean (which
I would like to use R to generate a histogram which has bars of variable bin
width with each bar having an equal number of counts. For example, if the
bin limits are the quartiles, each bar would represent 1/4 of the total
probability in the distribution.
An example of such an equal-probability
That's great, Peter.
For pedestrians like me who are not quite as facile with regular
expressions, the following seems slightly more readable:
s - the quick red fox jumps over the lazy brown dog
ss - strsplit(s, )[[1]]
ss1 - substring(ss, 1,1)
ss2 - substring(ss, 2)
(Ted Harding) [EMAIL PROTECTED] writes:
x - runif(500)
y - rbinom(500,size=1,p=plogis(x))
xx - predict(loess(resid(glm(y~x,binomial))~x),se=T)
matplot(x,cbind(xx$fit, 2*xx$se.fit, -2*xx$se.fit),pch=20)
Not sure my money isn't still on the splines, though.
.
Serves me right for
On 05-Dec-04 Peter Dalgaard wrote:
(Ted Harding) [EMAIL PROTECTED] writes:
x - runif(500)
y - rbinom(500,size=1,p=plogis(x))
xx - predict(loess(resid(glm(y~x,binomial))~x),se=T)
matplot(x,cbind(xx$fit, 2*xx$se.fit, -2*xx$se.fit),pch=20)
Not sure my money isn't still on the splines,
Dear lists,
I want to construct a loop in R, but don't know how to do it. I can do it
in SAS, but I prefer in R (which I am hoping I will off SAS for good
soon). Could anyone help me to convert the SAS codes to equivalent R codes.
Basically, the following codes were written to establish the
Thanks Peter.
Basically, I am trying to establish the line of marker allele the progeny
has inherited from either the sire or the dam. The progeny will share the
sires but each have different dam (ie., daughter design or half-sib model).
The SAS code I have written to identify the lines of
We've recently had a thread on beginners and FAQs and the like.
I decided that it might be a good idea to offer a short list of
on-line help pages that beginners should read. Goals:
- the list short be short (I think mine is too long)
- many of the commonest questions should find answers there
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of
[EMAIL PROTECTED]
Sent: Sunday, December 05, 2004 7:14 PM
To: [EMAIL PROTECTED]
Subject: Re: [R] What is the most useful way to detect
nonlinearity in lo
On 05-Dec-04 Peter Dalgaard wrote:
I'll really be brief this time...
From: Tony Plate
[snip]
Also, I think that John Maindonald's idea of a I am new to
R, where do I
start? page, with a link from the posting guide, is an
excellent idea.
Someone mentioned a tip shown at R startup (a la S-PLUS for Windows, I
guess). I
I am trying to store a couple numbers for each partition, in a subset of the
partitions
of my data set. Of course, one can accomplish this using a binary tree. (The
first split is on
inclusion/exclusion of the first object, and so on.) I can probably simulate
a tree using
vectors. (One
Have you had a look at the rpart package? If you haven't installed it it may be
worth doing so. Then you can type
require(rpart)
?rpart.object
Tom
-Original Message-
From: DFARRAR [mailto:[EMAIL PROTECTED]
Sent: Monday, 6 December 2004 11:09 AM
To: [EMAIL PROTECTED]
Subject: [R]
Hi all,
I'm a new user of R gam() function. I am wondering how do we decide on the
smooth function to use?
The general form is gam(y~s(x1,df=i)+s(x2,df=j)...) , how do we decide
on the degree freedom to use for each smoother, and if we shold apply
smoother to each attribute?
Thanks!!
Unfortunately that's not really an R question. I recommend that you read up
on the statistical methods underneath. One that I'd wholeheartedly
recommend is Prof. Harrell's `Regression Modeling Strategies'.
[BTW, there are now two implementations of gam() in R: one in `mgcv', which
is fairly
DFARRAR DFARRAR at vt.edu writes:
:
: I am trying to store a couple numbers for each partition, in a subset of the
: partitions
: of my data set. Of course, one can accomplish this using a binary tree.
(The
: first split is on
: inclusion/exclusion of the first object, and so on.) I can
Again, thanks for all the replies and suggestions.
At this point, I guess the only fix is NOT to use libsunperf from sun
studio 8. If anyone really wants to use libsunperf, use sun one studio
7 or 9's instead.
I will also try and send a bug report to the sun folks (in case there
are people like
On Mon, 6 Dec 2004 13:36:37 +1300
[EMAIL PROTECTED] wrote:
Dear lists,
I want to construct a loop in R, but don't know how to do it. I can do it
in SAS, but I prefer in R (which I am hoping I will off SAS for good
soon). Could anyone help me to convert the SAS codes to equivalent R codes.
Thanks for the email. I will check that out
However when I was doing this :gam(y~s(x1)+s(x2,3), family=gaussian,
data=mydata )it gives me the error :
Error in terms.formula(formula, data = data) :
invalid model formula in ExtractVars
What does it mean ?
Thanks
-Janice
42 matches
Mail list logo