I am working on data that is in clusters (eg events from persons in
the family).So I have a clustered data set up and would like to build
a tree(regression /classification tree) taking into account this
clustered nature of the data.
Do anybody know how to do this or maybe the code to take into acc
Dear R users,
has anybody made an attempt to create the soil texture triangle graph in
R? For an example see here:
http://www.teachingkate.org/images/soiltria.gif
I would like to get the lines in black and texture labels in gray to
allow for plotting my texture results on top.
Any examples
Does anybody know of a R package to solve mixed-integer quadratic
optimisations with constraints?
Alternatively is anybody working with an R interface to ILOG CPLEX
which seems to be the 'best' commercial product in this area?
Finally, is there any relationship between the R-project and AMPL
The final version with the help of Gabor Grotendieck (thanks Gabor, very much!)
###
# EasieR - Package #
###
# Common function
er.make.table <- function(x,
start,
end,
h,
Thank you very much Ted! I have been looking at your simplification for
more then an hour, but I don't see how you did it.
Could you perhaps, if it is not to much work, explain me how you reduced
H? It would help me to understand what I am realy doing.
Looking at the result, it seems indeed tha
Dear all,
I'm interested in testing 2 nonlinear restrictions on coefficients of a nls
object. Is there a package for doing this? Something in the lines of `test(nls
object, res=c("res 1","res 2"),...)'
I only found the function delta.method in the alr3 library that calculates the
se of a single
Hello,
I'm a beginner to use R.
When use cancor to analysis two data set X,Y, for
example:
can<-cancor(X, Y)
then 5 components are returned, that is cor, xcoef,
ycoef, xcentre, ycentre;
The explained variance can be calculated by X %*%
can$xcoef or Y %*% can$ycoef, but how to alculate the
cano
I am using R 2.1 on Apple OS X.
When I get the ">" prompt, I find it works well with emacs commandline
editing. Keys like M-f C-k etc. work fine.
The one thing that I really yearn for, which is missing, is bracket
matching When I am doing something which ends in it is really
useful to have e
Hi Ajay
well ESS has such a facility.
However, I think Mathematica has a super scheme: unbalanced brackets
show up
in red, making them obvious.
This is particularly good for spotting wrongly interleaved brackets, as
in
([ blah di blah )]
in which case both opening braces are highligh
Laura Holt wrote:
Hi again!
I have a data.frame with the columns y, x1, x2, x3.
I would like to fit linear models with one variable at a time,
then 2 variables at a time, and then 3.
Makes me think of a power set.
Makes me think of irreproducible results if you use the output to select
a si
I am writing a function of weighted regression, as a procedure for
heteroskedasticity.
The function runs an auxiliary regression whose fitted values I assign
to fit, and then I go:
w <- 1/(exp(fit/2))
## Rerun the old regression ##
if(gls) {
wtd.model <- glm(model, weights=
> well ESS has such a facility.
>
> However, I think Mathematica has a super scheme: unbalanced brackets
> show up
> in red, making them obvious.
>
> This is particularly good for spotting wrongly interleaved brackets, as
> in
>
> ([ blah di blah )]
>
>
>
> in which case both opening brac
R-help,
I have lloked in the archives found no answer to how to round the line
joint.
I have usedthe arguments lnd, ljoin in par but I get no differences in
the plotting.
x=1:10
par(ljoin="round",lend="round")
plot(x,sin(x),type="l",lwd=2)
Any suggestions?
I run on a Windows XP machine.
> ve
Hi,
I would like to create a multi-page postscript file. How can I do that in R ? Is
it possible ?
Thanks for help.
Philippe
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http:
Hi, Laura:
Have you considered "regsubsets" in library(leaps)? Also, have you
done an R site search "www.r-project.org" -> search -> "R site search"
for something like "all subsets regression"?
Consistent with Frank Harrell's comment, comments, I suggest you
randomly permute your resp
Spencer Graves wrote:
Hi, Laura:
Have you considered "regsubsets" in library(leaps)? Also, have
you done an R site search "www.r-project.org" -> search -> "R site
search" for something like "all subsets regression"?
Or much simpler if you are running R-2.1.0, use RSiteSearch("all sub
Philippe Lamy wrote:
Hi,
I would like to create a multi-page postscript file. How can I do that in R ? Is
it possible ?
Yes, simply draw more than one plot. See ?postscript and its argument
"onefile", which already defaults to TRUE.
Uwe Ligges
Thanks for help.
Philippe
___
What kind of nonlinear restriction? Can you solve for one or more of
the parameters in terms of the other(s) [either directly or implicitly]?
If yes, then let
fit1 <- nls(... full model ... )
fit2 <- nls(... restricted model ...)
anova(fit1, fit2)
If my m
Luis Ridao Cruz wrote:
R-help,
I have lloked in the archives found no answer to how to round the line
joint.
I have usedthe arguments lnd, ljoin in par but I get no differences in
the plotting.
x=1:10
par(ljoin="round",lend="round")
plot(x,sin(x),type="l",lwd=2)
Any suggestions?
Well, ro
HI all
i have a seemingly simple question.
given a sequence of numbers say, 1,2,3,4,5.
i would like to get all of the possible two number arrangments
(combinations), all 3 number arrangents ... 5 number arrangements
(combinations).
i.e.
in the 2 number case:
12,13,14,15,23,24,25,34,35,45
hai, I want to perform a simulation like this:
Suppose that I have one population ,it's size is 500, is composed of x,y
and z. The probability of x, y and z is respectively is 0.3, 0.5, 0.2. I
wan to simulate a new same size population based ratio of x, y and z, how
can I get and assess the n
Hi Jim,
This looks impressive! It gives me the 'background' graph. However, I'm
not sure how I can use this function to plot my soil texture values! Can
you explain?
I would like to be able to plot my soil texture samples in the same
graph as the one your function plots.
Maybe I should try
Hi Jari,
I assume this has been superseded by the ternaryplot{vcd} function!?
Thanks,
Sander.
Jari Oksanen wrote:
> Sander,
>
> Just a quick note before I go to the field.
>
> I attach a tri.R file for drawing ternary plots. The base function was
> posted to R News someday. One thing that I ad
Clark Allan wrote:
HI all
i have a seemingly simple question.
given a sequence of numbers say, 1,2,3,4,5.
i would like to get all of the possible two number arrangments
(combinations), all 3 number arrangents ... 5 number arrangements
(combinations).
i.e.
in the 2 number case:
12,13,14,15
luan_sheng wrote:
>
> hai, I want to perform a simulation like this:
>
> Suppose that I have one population ,it's size is 500, is composed of x,y
> and z. The probability of x, y and z is respectively is 0.3, 0.5, 0.2. I
> wan to simulate a new same size population based ratio of x, y and z, h
here is a possible solution:
subsets <- function (n, r, s = 1:n) {
# from V&Rs (2000) S Programming, Springer, p.49
if (mode(n) != "numeric" || length(n) != 1 || n < 1 || (n%%1) !=
0) stop("bad value of n")
if (mode(r) != "numeric" || length(r) != 1 || r < 1 || (r%%1) !=
0) stop("bad v
On 5/26/05, Thomas Lumley <[EMAIL PROTECTED]> wrote:
> If you *want* to fit mixed models (eg because you are interested in
> estimating variance components, or perhaps to gain efficiency) then it's
> quite a bit trickier. You can't just use the sampling weights in lme().
> You can correct for the
rmultinom(n=1,size=500,prob=c(0.3,0.5,0.2))
to get n samples each of size 500 just use the "n=" argument.
Reid Huntsinger
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of luan_sheng
Sent: Friday, May 27, 2005 9:26 AM
To: r-help@stat.math.ethz.ch; R-help@
Dear List,
I have been trying to use library sm to do density estimation on a 3D
dataset. I am using the current MacOS X binary of sm from CRAN. If I do
this on a 2D dataset, sm.density returns a list including the component
"estimate" which contains the density estimate over a uniform grid. Wh
Thank you for your reply.
Does that mean that in order to take in account the repeated measures I denote
these as another cluster in R?
Dassy
Quoting Thomas Lumley <[EMAIL PROTECTED]>:
> On Thu, 26 May 2005 [EMAIL PROTECTED] wrote:
>
> >
> > Dear R-Users!
> >
> > Is there a possibility in R
Is there a proper function for transforming
a character to a number instead of using
i=match('c',letters)
# 3
Thanks.
Josef Eschgfäller
--
Josef Eschgfäller
Dipartimento Matematico
Universita' di Ferrara
http://felix.unife.it
__
R-help@stat.math.
Hi,
Consider the following code:
require(XML)
xmlFile = paste( "\n",
"\n",
"\n",
"\n")
cat(xmlFile)
a = function(x,...){
cat("Attributes of ", xmlName(x), ": ");
I am trying to construct a graph of 6 box plots of blood pressures. I
want them to be on a single set of axis and I want the SBP to be ontop
of the DBP. I have an array bp with the data in it and I tried
a[1,]<-c(145,60,147,62,140,57)
a[2,]<-c(160,75,160,74,160,70)
a[3,]<-c(140,55,140,65,142,55
"Luis Ridao Cruz" <[EMAIL PROTECTED]> wrote in message
news:[EMAIL PROTECTED]
> R-help,
>
> I have lloked in the archives found no answer to how to round the line
> joint.
>
> I have usedthe arguments lnd, ljoin in par but I get no differences in
> the plotting.
>
> x=1:10
> par(ljoin="round",lend=
I searched the archives, but couldnt find any way to do this with
boxplot. Is there a way? Thanks again ~BJ
__
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
Josef,
Not sure if this is exactly what you mean, but there is a generic function
as()
to which you can then specify "numeric" as an argument and which then
coerces stuff into numeric format.
Tobias
Josef Eschgfaeller wrote:
Is there a proper function for transforming
a character to a nu
as.numeric()
hih
__
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
Josef Eschgfaeller wrote:
Is there a proper function for transforming
a character to a number instead of using
i=match('c',letters)
# 3
I'd suggest to use the above if you really mean it. Note that
"transforming a character to a number" is not well defined, because you
have to think abo
BJ wrote:
I searched the archives, but couldnt find any way to do this with
boxplot. Is there a way? Thanks again ~BJ
See ?boxplot and its argument par which points you to ?bxp:
Now, you can do stuff like
boxplot(1:10, pars=list(staplelty=0, whisklty=0))
Uwe Ligges
_
BJ wrote:
I am trying to construct a graph of 6 box plots of blood pressures. I
want them to be on a single set of axis and I want the SBP to be ontop
of the DBP. I have an array bp with the data in it and I tried
Folks, please invest 1 minute of time to rethink whether other people
will und
Hello,
I do not have any statistical background, So I shall
apologise if I am asking trivial question or help.
I am trying to work with R.
The problem I have on hand is:
I have 2 sets of data(means & SD for each sample in
the group of both sets). The sample size is
massive(2000+ in each grp).
I ha
Hmm ok, that seems to work with that example, but when I add an "at="
option i get graphs that are only verticle lines. Sorry to be difficult. ~BJ
Uwe Ligges wrote:
BJ wrote:
I searched the archives, but couldnt find any way to do this with
boxplot. Is there a way? Thanks again ~BJ
See ?b
BJ wrote:
Hmm ok, that seems to work with that example, but when I add an "at="
option i get graphs that are only verticle lines. Sorry to be difficult.
He? So you need to pe much more precisely in this message as well!
boxplot(data.frame(1:10, 2:11), at = 1:2,
pars=list(staplelty=0, whisk
Does it matter what they are? they are just names. The six box plots are
from the array I created with columns. I forgot to add teh dim
statement, my mistake. But I thought it was obvious that i was using a
matrix from the data frame call and the assignments I provided. My
question was simply a
as.numeric
Some weeks ago one discussed on the list about how
to transform a hexadecimal representation to
decimal digits. Specifically one has to transform
'a' to 10 etc. In C one does this with 'a'-87,
but I did not find a function for this in R.
Using match on letterdigits as in the followi
On Fri, 2005-05-27 at 10:17 -0400, BJ wrote:
> I am trying to construct a graph of 6 box plots of blood pressures. I
> want them to be on a single set of axis and I want the SBP to be ontop
> of the DBP. I have an array bp with the data in it and I tried
>
> a[1,]<-c(145,60,147,62,140,57)
> a[2,
Read this book, "Multivariate Statistical Analysis: A Conceptual
Introduction" by Sam Kash Kachigan. I think it's *great*, and perfect
for someone without any statistical background.
-Original Message-
From: manav ram [mailto:[EMAIL PROTECTED]
Sent: Friday, May 27, 2005 10:56 AM
To: r-hel
--- Begin Message ---
> boxplot(data.frame(1:10, 2:11), at = 1:2,
+ pars=list(staplelty=0, whisklty=0), at=1:2)
Error in boxplot.default(data.frame(1:10, 2:11), at = 1:2, pars =
list(staplelty = 0, :
formal argument "at" matched by multiple actual arguments
a<-array(dim=c(3,8))
a[1,
BJ-
For the record as well I could not make heads or tails of your question.
The dedicated folks who day in, day out field poorly specified questions have
no obligation to do so, and have every right to get frustrated when it seems
that their time is being taken advantage of.
On Fri, 27 May 2005
I just asked a question. If I was too vague, then i am sorry. I dont
expect anyone to help me, but I thought that it was ok to put the
question out there in case someone wanted to help me. I didnt expect
abject hostility for it. Human decency was the only thing I did expect.
If my question was
On Fri, 27 May 2005 [EMAIL PROTECTED] wrote:
Thank you for your reply.
Does that mean that in order to take in account the repeated measures I denote
these as another cluster in R?
Yes, but unless you have multistage finite population corrections to put
in the design object only the first
Right,
Got the data points plotted on top of the soil texture background,
thanks to Jim and ternaryplot{vcd}! See code below.
Now there is some fine tuning to do, as it should really look like this
graph:
http://soil.scijournals.org/content/vol65/issue4/images/large/1038f2.jpeg
Things to do
Hi!
I want to convert S-Plus 6.2 code to R 2.1.0. Instead of the function nlminb I
use the function optmin
optmin(start,fn,gr,method="L-BFGS-B", lower, upper, hess,...)
But then I get the Error in optmin ...: L-BFGS-B needs finite values of fn
Then I used optmin(start,fn,gr,method="BFGS", hess
If you are running R in the *R* buffer inside emacs, then you automatically
have highlighting of mismatched parentheses and brackets.
Rich
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting g
Of course it should be in ESS. In fact, it is for me, anyway
(different colors depending on whether they match or not).
Check out the paren-hilit or paren-match (or something like that)
customize options.
I.e.
M-x customize-groups paren
(which ought to complete on names of groups/options st
Stefan Pohl wrote:
Hi!
I want to convert S-Plus 6.2 code to R 2.1.0. Instead of the function nlminb I
use the function optmin
optmin(start,fn,gr,method="L-BFGS-B", lower, upper, hess,...)
But then I get the Error in optmin ...: L-BFGS-B needs finite values of fn
Then I used optmin(start,fn
Stefan Pohl wrote:
> Hi!
>
> I want to convert S-Plus 6.2 code to R 2.1.0. Instead of the function nlminb
> I use the function optmin
>
> optmin(start,fn,gr,method="L-BFGS-B", lower, upper, hess,...)
>
> But then I get the Error in optmin ...: L-BFGS-B needs finite values of fn
>
> Then I used
Sorry, still confused. If I dont have fpc's ready in my dataset (calculate
myself?) that means that R will use the weight of an individual for each of his
repeated observations. But is that then still correct? The "cluster" individual
is ignored and each observation of an individual has the same we
Cleaned up the class divisions and created a full function.
Still to do:
- rotate axis labels;
- correct the partial covering of the bottom tick labels;
- rotate ticks in order to simplify viewing the graph.
See:
http://soil.scijournals.org/content/vol65/issue4/images/large/1038f2.jpeg
Wonder
In an earlier message today I mentioned that a source package with a
version of nlminb for R was available as
http://www.stat.wisc.edu/~bates/port_0.1-1.tar.gz
Thanks to Kjetil Halvorsen there is now a Windows binary version
available as
http://www.stat.wisc.edu/~bates/port_0.1-1.zip
__
Hello, Erithid:
I'm very sorry you feel you got yelled at first thing in the morning.
Please try not to take it personally, though I know that may be
difficult. After one superficially insulting reply, I just laughed and
told my manager, "Look what I learned for exposing myself to ridicul
Hi:
I have a question arising from my project.
A sample of the data is below. The first row stands
for the names of state in USA. The second row stand
for some numeric value in that state. Some of them are
NA. I can use the commands "data(stateMapEnv)" and
"map('state', fill = F)" in library "
Friends,
I am trying to install the current version of spatstat
on Mac OS 10.3.9, but the compilation fails with the
following messages at the end:
ld: warning -L: directory name
(/usr/local/lib/gcc/powerpc-apple-darwin6.8/3.4.2)
does not exist
ld: can't locate file for: -lg2c
make: *** [spa
This is is not difficult from the online help for the map() function.
Here is an example.
map('state', region = c('new york', 'new jersey', 'penn'),fill=TRUE,col=1:4)
There is also an example there in how to add text to the map.
Another way uses the maptools package.
require(maptools)
?plot.M
On Fri, 27 May 2005 [EMAIL PROTECTED] wrote:
Sorry, still confused. If I dont have fpc's ready in my dataset (calculate
myself?) that means that R will use the weight of an individual for each of his
repeated observations. But is that then still correct? The "cluster" individual
is ignored and e
To add to Spencer's reply, I have two things.
First, if you look back at your earlier list entries, Uwe asked you
politely at least twice to read the posting guide. The reason for this
is that the posting guide gives guidance to allow you to post questions
that are clear enough to answer, and poin
Perhaps the list would work better if posters who do not follow
the posting guide are nicely told that they will more likely get
helpful replies if they read and follow the posting guide and,
in particular, provide a reproducible example.
If the poster does not follow that advice then the poster c
67 matches
Mail list logo