Ulrich Leopold [EMAIL PROTECTED] writes:
Dear list,
could someone point me to the right command to subtract 2 columns in a
data.frame. Might be a bit embarrassing question. But I cannot figure
out how to do this simple command in R.
E.g.,
mydata$difference - mydata$x - mydata$y
--
Is the following what you want?
DF - data.frame(a=1:2, b=3:4)
DF$a-DF$b
[1] -2 -2
DF[,a]-DF[,b]
[1] -2 -2
hope this helps. spencer graves
Ulrich Leopold wrote:
Dear list,
could someone point me to the right command to subtract 2 columns in a
data.frame. Might be a bit embarrassing question.
Hi,
I'm using the boot package to bootstrap a linear model. The boot command
is:
blm01 - boot(mat, boot.fishpower, 1000, strata=boot.strata)
I'm having several problems with the linear model so I decided to check
the stratification. The first thing I did was checking the number of
observations
On 4 Sep 2003, Ulrich Leopold wrote:
could someone point me to the right command to subtract 2 columns in a
data.frame. Might be a bit embarrassing question. But I cannot figure
out how to do this simple command in R.
Suppose your data frame is called foo, and you want the first column minus
I think the Mozilla plugin site explicitly says that you cannot copy the
Java plugin into the Mozilla plugin directory. You have to make a symbolic
link to it in the plugin directory. Please try that.
Jason Turner [EMAIL PROTECTED]
Sent by: [EMAIL PROTECTED]
09/03/2003 07:32 PM
Vadim == Vadim Ogranovich [EMAIL PROTECTED]
on Wed, 3 Sep 2003 14:29:25 -0700 writes:
Vadim Hi, I thought it would be convenient if the
Vadim check.names argument to read.table, which currently
Vadim can only be TRUE/FALSE, could take a function value
Vadim as well. If the
Hi
I have two geostatistical models from geoR. An ordinary kriging model with
AIC=-148.6 and a universal kriging model with AIC=-156.7, there are 345
data points. The improvement shown by the AIC by adding a trend component
to the model seems quite small given the number of data points, is
E.g.,
mydata$difference - mydata$x - mydata$y
That's what I thought, but I get the following message:
propLSK.STONE.Pox0t30$Pox0t30STONE-propLSK.STONE.Pox0t30$Pox0t30
numeric(0)
Does it mean the resulting vector is empty? If yes, what could be the reason
for it? Both columns are available:
2*log(likelihood ratio) is approximately chi-square for nested
models. AIC = (-2)*(log(likelihood)-k), where k = number of parameters
in the model.
Thus, del(AIC) = 2*(log(likelihood ratio)-del(k)). If the trend is
strictly linear, then it involves only 1 parameter, so del(k) = 1.
Dear list,
sorry for having bothering you. I found the problem. It was the a stupid
error made by me. R could not recognise the variable as it is of course case
sensitive for names. So I specified the wrong variable name and R did
complain about it with numeric(0).
Ulrich
Peter Dalgaard BSA
Ulrich Leopold [EMAIL PROTECTED] writes:
E.g.,
mydata$difference - mydata$x - mydata$y
That's what I thought, but I get the following message:
propLSK.STONE.Pox0t30$Pox0t30STONE-propLSK.STONE.Pox0t30$Pox0t30
numeric(0)
Does it mean the resulting vector is empty? If yes, what
Hi
On 4 Sep 2003 at 14:35, Ulrich Leopold wrote:
E.g.,
mydata$difference - mydata$x - mydata$y
That's what I thought, but I get the following message:
propLSK.STONE.Pox0t30$Pox0t30STONE-propLSK.STONE.Pox0t30$Pox0t30
numeric(0)
Beware
The officially sanctioned way to put the expression lambda_1 = x in a title
is something like this:
title(substitute(lambda[1] == lamb, list(lamb = x)))
but suppose I have two lambdas and would like something like
lambda_1 = x_1 , lambda_2 = x_2
to appear. What then?
- Original Message -
From: Richard A. O'Keefe [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Sent: Thursday, September 04, 2003 2:56 AM
Subject: Re: [R] Overlaying graphs
I do not know how to overlay the curve graphic on top of hist graphic.
Do you know about the add=TRUE option for
On Thu, 4 Sep 2003, Roger Koenker wrote:
The officially sanctioned way to put the expression lambda_1 = x in a title
is something like this:
title(substitute(lambda[1] == lamb, list(lamb = x)))
but suppose I have two lambdas and would like something like
lambda_1 = x_1 ,
Hi people,
I have made some progress trying to work out how to solve this problem
but I have got a bit stuck - sorry if this turns out to be a simple
exercise . .
Allelic Differentiation (AD) in genetics measures the number of
different alleles between (say) two populations eg:
Organisms
Please note Student-Newman-Keuls is NOT a recommended multiple comparison procedure.
In the language of Hsu (1996), Student-Newman-Keuls is not even a confident
inequalities method. In other words, it does not control the probability of making at
least one incorrect assertion of inequality
Dear users,
is anybody of you aware of a R command to perform laplace transform or
even its inversion?
Thank you very much.
Luca
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Hi there,
I have four panels in a lattice bwplot. I want to have
two different ylim for the panels, for example panels
[1,1] and [1,2] with ylim=c(0,200) and panels [2,1]
and [2,2] with ylim=c(0,100).
Thanks for help in advance.
Mahbub.
__
[EMAIL
Your observation that qqnorm does not appear to be very general is
rebutted by Venables and Ripley (2002) Modern Applied Statistics with S,
4th ed. (Springer, p.108): One of the best ways to compare the
distribution of a sample x with a distribution is to use a Q-Q plot. ...
This idea can
Paul Meagher wrote:
2. Does R have a suite of best-fit tools for finding the best
fitting-probability distribution for any observed probability distribution?
I think that the best-fitting probability distribution for an observed
probability distribution is the empirical distribution of your
Hi,
I found the problem, now it is working fine. cheers Martin
On Thursday 04 September 2003 16:20, Martin Wegmann wrote:
Hello R user,
I have several data frames with 100 columns and I did a linear regression
over time of each column
df1.lm - lapply(df1, function(x) lm(x~year)$coeff[2])
Hello there fellow R-users,
I am using the RODBC functionality to query a database. I am trying to read
in a columns of strings which have a character field lengths greater than
255.
The data.frame that I recieve back from the RODBC query only contains the
first 255 characters (the rest having
On Thu, Sep 04, 2003 at 05:37:46PM +0200, Brunschwig, Hadassa {PDMM~Basel} wrote:
Is there an equivalent in R for the dos() command in S?
help(system)
Hth, Dirk
--
Those are my principles, and if you don't like them... well, I have others.
--
I have a personal rule that any function I write that is longer than
about one page should be modularized. This rule pushes me to think
carefully about how to structure the problem. It also helps me debug my
code and modify it as appropriate. If I can develop appropriate
primatives, I
Sorry Folks,
I'm sure I could suss out the answer myself but I need it
soon ... !
1. Given a set of 4 variables X,Y,Z,W in a dataframe DF, I make
a scatter-plot matrix using splom(DF).
2. I do all regressions of U on V using lm(U~V), where U and V
are all 12 different ordered pairs from
Marcel Vieira wrote:
Thanks a lot for your help.
But I am getting that error message when I am trying to
load my function. Then I can't use anything to minimise it.
Is there anything I can do to increase R's
capacity for loading very long functions?
I assume by 'load' you must be somehow using
XLSolutions Corporation ([1]www.xlsolutions-corp.com) is pleased to
announce a two-day course, Splus/R: Complementing and Extending
Statistical Computing for SAS Users
[2]www.xlsolutions-corp.com/Rsas.htm
Dates/City: October 9-10, 2003 in Washington DC
October 2-3,
Would S Poetry be the sort of thing you are looking for?
Patrick Burns
Burns Statistics
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and A Guide for the Unwilling S User)
Rado Bonk wrote:
Hi,
Does anybody know suitable documents (manuals) on writing user
Hello
When I run
scatter.smooth(jitter(weight), jitter(height2), span = .25, evaluation
= 50, pch = '.')
I get the type of graph I thought I would get, but also a warning.
k-d tree limited by memory. ncmax= 528
I always get concerned when there are warnings I don't understand.
What's
How does one get counts for grouped data
from ungrouped data. For example, how can
I get
X1 X2 Count
1 1 1
1 2 2
2 1 1
2 2 1
from
X1 X2
AA
AB
BB
AB
BA
For example (variations are possible),
y - 100 * runif(400)
a - gl(4, 100)
x - gl(3, 7, 400)
library(lattice)
bwplot(y ~ x | a, scales = list(y = free),
ylim = list(c(0, 100), c(0, 200), c(0, 100), c(0, 200)))
(Assuming you have up to date versions of R and lattice).
On Thursday 04
On Thu, 4 Sep 2003, Paul Green wrote:
How does one get counts for grouped data
from ungrouped data. For example, how can
I get
X1 X2 Count
1 1 1
1 2 2
2 1 1
2 2 1
from
X1 X2
AA
I admit I should have been more clear in my original posting. Let me try again (and I
do know that by deafulat read.table discards everything after '#' which is why I use
comment.char=, my bad not to mention this).
Here is a typical example of my data file:
#keyvalue
foo 1.2
boo
I am one of only 5 or 6 people in my organization making the
effort to include R/Splus as an analysis tool in everyday work -
the rest of my colleagues use SAS exclusively.
Today, one of them made the assertion that he believes the
numerical algorithms in SAS are superior to those in Splus
and R
Does any of you know where I can find an explanation of lapack errors codes?
I get error code 17 when using optim().
Is there a way to handle errors in R such that depending on the type of
error I can decide what to do next?
Thanks,
Haky
__
[EMAIL
On Thu, 2003-09-04 at 08:34, Frank E Harrell Jr wrote:
On Thu, 04 Sep 2003 14:50:25 -0400
Paul, David A [EMAIL PROTECTED] wrote:
I am one of only 5 or 6 people in my organization making the
effort to include R/Splus as an analysis tool in everyday work -
the rest of my colleagues use
On Thu, 4 Sep 2003, Spencer Graves wrote:
Have you considered the two books by Venables and Ripley (2002) Modern
Applied Statistics in S and (2000) S Programming (both Springer)? If
yes, I don't know what you mean by the usually available manuals.
In addition, John Chambers's Programming
I just did a Google search for lapack error codes and got several
hits. However, I don't know if a translation of error code 17 would
really help you.
Have you tried fn as a stand-alone function? If yes, have you
tried including a cat or print statement in it, so it prints the
value
Paul, David A wrote:
I am one of only 5 or 6 people in my organization making the
effort to include R/Splus as an analysis tool in everyday work -
the rest of my colleagues use SAS exclusively.
Today, one of them made the assertion that he believes the
numerical algorithms in SAS are superior to
yes: From ?optim, I get the following:
optim(par, fn, gr = NULL,
method = c(Nelder-Mead, BFGS, CG, L-BFGS-B, SANN),
lower = -Inf, upper = Inf,
control = list(), hessian = FALSE, ...)
By fn, I meant the function you supply for that argument.
spencer graves
Paul, David A [EMAIL PROTECTED] writes:
I am one of only 5 or 6 people in my organization making the
effort to include R/Splus as an analysis tool in everyday work -
the rest of my colleagues use SAS exclusively.
Today, one of them made the assertion that he believes the
numerical
Douglas Bates [EMAIL PROTECTED] writes:
McCullough, B. D. (1998), Assessing the reliability of statistical
software: Part I, The American Statistician, 52, 149-159.
McCullough, B. D. (1999), Assessing the reliability of statistical
software: Part II, The American Statistician, 53,
On Thu, 4 Sep 2003, Paul, David A wrote:
I am one of only 5 or 6 people in my organization making the
effort to include R/Splus as an analysis tool in everyday work -
the rest of my colleagues use SAS exclusively.
Today, one of them made the assertion that he believes the
numerical
On Thu, 2003-09-04 at 12:00, Vadim Ogranovich wrote:
I admit I should have been more clear in my original posting. Let me
try again (and I do know that by deafulat read.table discards
everything after '#' which is why I use comment.char=, my bad not to
mention this).
Here is a typical
My thanks to Drs. Armstrong, Bates, Harrell, Liaw, Lumley,
Prager, Schwartz, and Mr. Wang for their replies. I have
pasted my original message and their replies below.
After viewing http://www.itl.nist.gov/div898/strd/ as suggested
by Dr. Schwartz, it occurred to me that it might be educational
Thanks, Deepayan!
However, for me this has deepened the mystery (I don't really
understand in detail how lattice graphics works anyway!).
To clarify: The variables X,Y,Z,W in DF have some zero values,
and otherwise are positive. For U,V in X,Y,Z,W I plot log(1+V)
against log(1+U) for all the
By the way, Dong,
you shouldn't use the term mean in Splus to designate the mean of your trait
(mean-8000). mean' is an Splus function, and by doing mean-8000, you may experience
problems when you'll need the argument mean to designate the original Splus function
(as in the function
On Thursday 04 September 2003 19:52, Ted Harding wrote:
Thanks, Deepayan!
However, for me this has deepened the mystery (I don't really
understand in detail how lattice graphics works anyway!).
To clarify: The variables X,Y,Z,W in DF have some zero values,
and otherwise are positive. For
Hi,
I am a new user of R and am building an R add-on package. I followed the
Writing R Extension manual from cran website but still met some
problems that I cannot solve. I build it under redhat linux 9.0 R1.7.0.
Say foo is the name of the package.
First, after I wrote the Rd files for each
Dear Luca,
I don't think that R has a built-in function for doing
Laplace or inverse Laplace transforms. I remember having to
use an IMSL routine (INLP, I think) to do this many years
ago. When I looked at the article that the algorithm was
based on, I found that as an example the author showed
51 matches
Mail list logo