hi peter,
thanks for your reply. i was kind of hoping/thinking that others have
encountered the same problem too.
i try to provide a reproducable example.
my variables are the following (just the cases 1:50). the whole dataset has
about 2800 cases.
satlife[1:50]
[1] 9 10 11 8 11 10 11 6
Dear expeRts,
I would like to use expression() for creating labels in a splom, as shown in
the
first minimal example below. Is there any way I can simplify having to write
expression(italic(...)) several times?
The second example is what I tried so far, but I can't manage to get italic() to
I must be missing something.
I first state my g parameter with:
lamda-c(g=0.2)
However, when I do the next step R is telling me object g not found Here is
my next step:
Q-exp(g)
???
Any help would be greatly appreciated.
Mike
On Thu, Nov 25, 2010 at 09:23:17PM -0800, Stephen Liu wrote:
Hi David,
But you didn't try:
DNase# which was after all the name of the object you saved.
Sorry I don't follow.
He is telling you that it is not surprising that 'aaa' does not exist,
if the object you saved was called
Euu...
g - 0.2
Q - exp(g)
or
lamda - c(g=0.2)
exp(lamda)
Ivan
Le 11/26/2010 09:12, Mike Gibson a écrit :
I must be missing something.
I first state my g parameter with:
lamda-c(g=0.2)
However, when I do the next step R is telling me object g not found Here is
my next step:
Q-exp(g)
Hello Mike.
I'm not clear on why you would want to state your g parameter (In fact, I
don't even know for sure what you mean by that) with
lamda-c(g=0.2)
If you want a variable (vector) g containing 0.2, why don't you simply do:
g-0.2
If you need that lamda thing for some reason later on, you
Hi,
He is telling you that it is not surprising that 'aaa' does not exist,
if the object you saved was called DNase
No, such a file. Just rechecked it with Win7 search command. (Win7 has been
rebooted)
Again - you need to use the name of the object which happens to be
'DNase' - not
Thanks a lot, Prof. Ripley. the problem must be download.file() prior
R-patched 2.12.
for each loop, I tired a couple of candidate links and only one or
none will work, this must accumulate soon to an amount of unclosed
destination files beyond the tolerance of windows.
I updated R from 2.10 to
Dear R-experts,
Using persp3d and surface3d, I have two surfaces. I would like to plot the
projection of interesection of these two planes on X-Y plane -- such that I
could see a vertical 3D plane.
Any suggestions in this regard?
Thank you.
Regards,
Suresh
--
View this message in context:
Dear R-experts,
Using persp3d(), I plotted a 3d surface. I would like to colour this surface
based on the result values (dependent variable) grouped in, say, three to
five different ranges. Later, I would also like to paint the same surface
based on the range-values of one of the independent
Dear r users,
I wanted to perform contour figure from response surface model result. I typed
this:
model = rsm (Index ~ SO(NV, SKLON, pH), bonita)
summary(model)
par(mfrow = c(2, 3))
contour (model, ~NV+SKLON+pH)
However, I got the error answer as follows:
Error in lmobj$data[, nm] :
On Thu, 2010-11-25 at 22:56 +, Chris Mcowen wrote:
Hi,
I am not sure if it is more robust than a discriminant function but it
is certainly capable if differentiating between species based on
morphology. I used 12 measurements in my fish.
PCA extracts linear combinations of the data
He is telling you that it is not surprising that 'aaa' does not exist,
if the object you saved was called DNase
No, such a file. Just rechecked it with Win7 search command. (Win7 has been
rebooted)
You are confusing two things here:
a) Files living in the filesystem
b) Objects in your
Try this:
which.min(sapply(mylist, min))
On Thu, Nov 25, 2010 at 10:23 PM, Lorenzo Cattarino
l.cattar...@uq.edu.auwrote:
Hi R-users,
I have a list
mylist - list(c(0.79, 0.92, 0.91, 0.86, 0.96, 0.96, 0.95, 0.94, 0.99),
c(0.28, 0.45, 0.59, 0.69, 0.80, 0.87, 0.95, 0.94, 0.98), c(0.29,
On 26/11/2010 5:21 AM, suresh wrote:
Dear R-experts,
Using persp3d(), I plotted a 3d surface. I would like to colour this surface
based on the result values (dependent variable) grouped in, say, three to
five different ranges. Later, I would also like to paint the same surface
based on the
Hi! My regression is lm(Y ~ A + B + C + D + E + F + G), where each covariate
A--G is of the kind as.factor(column of the dataframe).
I receive an error for the first explanatory variable A: Error in column
of the dataframe: wrong number of dimensions. What does it mean?
--
View this message in
On 2010-11-25 23:59, Manderscheid Katharina wrote:
hi peter,
thanks for your reply. i was kind of hoping/thinking that others have
encountered the same problem too.
i try to provide a reproducable example.
my variables are the following (just the cases 1:50). the whole dataset has
about
Date: Thu, 25 Nov 2010 21:22:54 -0600
To: santosh.srini...@gmail.com
From: e...@debian.org
CC: r-help@r-project.org
Subject: Re: [R] Installing RQuantLib on Win 7 64 Bit
On 26 November 2010 at 07:05, Santosh Srinivas wrote:
| Hello Group,
|
Hi
I have 1800 response variables to regress on two factors (latitude and age),
what is the script to run all response variables at once instead of writing
1800 models?
Thanks
R.
Lähettäjä: r-help-boun...@r-project.org [r-help-boun...@r-project.org]
Hi,
It probably means that somehow your data is stored incorrectly. Try
using str(name of the dataframe) and see if everything is named what
you expect and is the class you think it is with equal numbers of
observations.
HTH,
Josh
On Fri, Nov 26, 2010 at 3:39 AM, effeesse
Hello!
I am trying to extract a p-value from an aov-object. I´ve done some research
on the internet and found a number of solutions, but somehow, none seems to
work for me.
This is the Object:
summary(rt.aov)
Error: pbn
Df Sum Sq Mean Sq F value Pr(F)
Residuals 9 378283 42031
great !
looks good. I am sure, it'd work.
Thanks a lot !
Regards,
Suresh
--
View this message in context:
http://r.789695.n4.nabble.com/3D-surface-colour-based-on-the-values-of-X-or-Y-data-points-tp3060018p3060160.html
Sent from the R help mailing list archive at Nabble.com.
Ok. So, I assume that colours based on the value-range of independent
variables would be similar.
Will try..
--
View this message in context:
http://r.789695.n4.nabble.com/3D-surface-colour-based-on-the-values-of-X-or-Y-data-points-tp3060018p3060162.html
Sent from the R help mailing list
Sorry, not very familar with the build process ... Can anyone with a
working 64 bit binary for rquantlib, pls email me a copy ...
Sent from my iPhone
On 26-Nov-2010, at 5:25 PM, Mike Marchywka marchy...@hotmail.com
wrote:
Date: Thu, 25 Nov
Worked fine Gabor, thanks!
I just removed the init argument from the Reduce function and I got it
without the first 0.
Henrique
-Mensagem original-
De: Gabor Grothendieck [mailto:ggrothendi...@gmail.com]
Enviada em: quinta-feira, 25 de novembro de 2010 22:13
Para: henrique
Cc:
On 2010-11-26 03:05, peter_petersen wrote:
Hello!
I am trying to extract a p-value from an aov-object. I´ve done some research
on the internet and found a number of solutions, but somehow, none seems to
work for me.
This is the Object:
summary(rt.aov)
Error: pbn
Df Sum Sq Mean Sq
Hi,
I've been a bit confused by different wyas we specify random effects
in lmer and MCMCglmm i just want to clear something. When I want to
look for intersexual genetic correlations in the trait, is it
equivalent to treat this trait for opposite sexes as separate traits
and include the term
thanks for your suggestions.
i think i am getting some sort of idea where the problem lies.
i am working with a dataset which i attached. the objects of my regression
model are derived from this dataset but, as far as i understand, only exist in
the r environment.
i have detached the main
# The result I am after is the result after a substitution in an expression,
such as
substitute(expression(a+b+c), list(a=1))
expression(1 + b + c)
# However, the way I want to do it is for a an expression stored as a
variable as
(expr - expression(a+b+c))
expression(a + b + c)
# a) The
On Fri, Nov 26, 2010 at 8:31 AM, Søren Højsgaard
soren.hojsga...@agrsci.dk wrote:
# The result I am after is the result after a substitution in an expression,
such as
substitute(expression(a+b+c), list(a=1))
expression(1 + b + c)
# However, the way I want to do it is for a an expression
it is working!
Thank you very much.
--
View this message in context:
http://r.789695.n4.nabble.com/extract-p-value-from-an-aov-oject-tp3060079p3060225.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
Just wondering if the proposed loop for the metadata on this conversion was
implemented by any chance...before I go trying the Long routine...
--
View this message in context:
http://r.789695.n4.nabble.com/Reading-sas7bdat-files-directly-tp1469515p3060319.html
Sent from the R help mailing list
Thanks, but I think there is a small mistake in your code:
expr - expression(a+b+c)
do.call(substitute, list(expr, list(a=1)))
expression(a + b + c)
I think it should be:
as.expression(do.call(substitute, list(expr[[1]], list(a=1
expression(1 + b + c)
- or maybe it can be done in a
I found one error in in the code, however, I still receive errors.
Error in storage.mode(y) - double :
invalid to change the storage mode of a factor
Moreover: Warning message:
In model.response(mf, numeric) :
using type=numeric with a factor response will be ignored
From my understanding
Hello,
I am using the mlogit package to model diplomatic representation and
stability.
Here is our sample model -
m2 = mlogit(y ~ 0|(time.curr + period), data=dip,
alt.levels=c(asym,end,sym,x), shape=wide,
subset=which(mlogdata$curr==sym))
From here, we would need to find the probabilities
On Fri, Nov 26, 2010 at 8:48 AM, Søren Højsgaard
soren.hojsga...@agrsci.dk wrote:
Thanks, but I think there is a small mistake in your code:
expr - expression(a+b+c)
do.call(substitute, list(expr, list(a=1)))
expression(a + b + c)
I think it should be:
as.expression(do.call(substitute,
Dear All,
I am using the code below to calculate error bars. I note that the length of
the error bars can be varied by varying the constant (0.975). It does
appear that any number can be substituted for 0.975, making it confusing for
me to know how to quantify the error bars.
I wish to quantify
Hi everybody,
I would like to know if there is any CPU scheduling package in R. Anybody knows?
Thanks,
Romildo
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Dear R users,
i would like to catch error message (coming after
unsuccessful database query) so that the script will
process further. How can I manage this?
Thank you in Advance,
Alla.
__
R-help@r-project.org mailing list
see below.
On Wed, Oct 27, 2010 at 3:32 AM, Prof Brian Ripley
rip...@stats.ox.ac.uk wrote:
On Tue, 26 Oct 2010, Kjetil Halvorsen wrote:
I have been searching the documentation for some time nmow, but cannot
find it. It must be possible to download packages (many), but only install
the help
?try
On Fri, Nov 26, 2010 at 9:26 AM, Alla Bulashevska
alla.bullashev...@fdm.uni-freiburg.de wrote:
Dear R users,
i would like to catch error message (coming after
unsuccessful database query) so that the script will
process further. How can I manage this?
Thank you in Advance,
Alla.
Dear list,
running a hierachical cluster analysis I want to define a number of objects
that build a cluster already. In other words: I want to force some of the cases
to be in the same cluster from the start of the algorithm.
Any hints? Thanks in advance!
Derik
[[alternative HTML
Dear R-gurus:
I want to convert a zoo object of daily stock prices
to a weekly one, based on the following 4 steps.
1. No problem in creating weekly data using Wednesday Obs.
wed.index lt;- ( weekdays(index(daily)) == Wed )
week.dat lt;- daily[wed.index]
week.dat
A10 A20 A30
On Fri, Nov 26, 2010 at 3:30 AM, Marius Hofert m_hof...@web.de wrote:
Dear expeRts,
I would like to use expression() for creating labels in a splom, as shown in
the
first minimal example below. Is there any way I can simplify having to write
expression(italic(...)) several times?
The
1. survreg() does NOT fit a proportional hazards model, a mistake
repeated multiple times in your post
2. The coxph function operates on the risk scale: large values of Xbeta
= large death rates = bad
The survreg operates on the time scale: large values of xbeta =
longer liftetime = good.
3.
Hi,
I'm currently trying desperately to get the nnet function for training a
neural network (with one hidden layer) to perform a regression task.
So I run it like the following:
trainednet - nnet(x=traindata, y=trainresponse, size = 30, linout = TRUE,
maxit=1000)
(where x is a matrix and y a
try.result - try(your.expr)
if(inherits(try.result, try-error)){
do.what.needs.to.be.done.here()
} else {
cat(OK, sep=\n)
}
Or have a look at ?tryCatch
Cheers,
Janko
-Ursprüngliche Nachricht-
Von: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
Im
I add an example , all the variables are mutually excluding dummy variables,
notice the different intercept: 5.627 vs 5.545:
survreg:
Value Std. Error zp
(Intercept) 5.6270.00887 634.3 0.00e+00
Var1.recR2 -0.1080.01026 -10.5 1.00e-25
Var1.recR3 -0.490
Hello,
To add more detail to my last post in which I asked for help regarding
a complex ANOVA, I've attached my data and have included my most
recent attempts at analysis. Although I tried to refer to the maiz
example as Rich suggested (and thank you, Rich, for suggesting it), I
didn't understand
I am writing a function that breaks the output if it gets bigger than the
available console width.
I noticed that a right margin always remains. The following code shows that the
actual console size
is always bigger than the one reported by options()$width. What might be the
reason?
(Software:
I asked my local java export and his reply was that the driver was compiled
with the later version of java that i need to install. I did install it in
parallel to existing java and changed JAVA_HOME variable to point to new
java directory. I also pre-pended PATH with path to java. Still, when i
http://r.789695.n4.nabble.com/file/n3060425/fig_1.png fig. 1
http://r.789695.n4.nabble.com/file/n3060425/fig_2.png fig. 2
I want to the picture like the above one, the origin crossover together,
while the following picture can be obtained by default and the origin is
detached, but throgut
On 2010-11-26 06:50, bbslover wrote:
http://r.789695.n4.nabble.com/file/n3060425/fig_1.png fig. 1
http://r.789695.n4.nabble.com/file/n3060425/fig_2.png fig. 2
I want to the picture like the above one, the origin crossover together,
while the following picture can be obtained by default and
Dear,
I would like to generate kml file with lines (LineString) of
different style (definitely color, but width would also be nice to
see). However with kmlLine from maptools package I am able only to
output the single first Lines object (first row of
SpatialLinesDataFrame)
I guess the option
On Fri, Nov 26, 2010 at 5:03 PM, fbielejec fbiele...@gmail.com wrote:
Dear,
I would like to generate kml file with lines (LineString) of
different style (definitely color, but width would also be nice to
see). However with kmlLine from maptools package I am able only to
output the single
Danielwc daniel.carstensen at gmail.com writes:
Im using the VEGAN package to do an RDA ordination. In my plot I get my
environmental scores as arrows/vectors, but my species scores as points. I
would like to get the species scores as arrows as well.
Is there not a way I can tell R to plot
On Fri, Nov 26, 2010 at 05:32:54PM +0100, Mark Heckmann wrote:
I am writing a function that breaks the output if it gets bigger than the
available console width.
I noticed that a right margin always remains. The following code shows that
the actual console size
is always bigger than the one
Dnia 2010-11-26, o godz. 17:11:22
Barry Rowlingson b.rowling...@lancaster.ac.uk napisał(a):
On Fri, Nov 26, 2010 at 5:03 PM, fbielejec fbiele...@gmail.com
wrote:
Dear,
I would like to generate kml file with lines (LineString) of
different style (definitely color, but width would also be
On Fri, Nov 26, 2010 at 6:55 AM, Derik Burgert derik2...@yahoo.de wrote:
Dear list,
running a hierachical cluster analysis I want to define a number of objects
that build a cluster already. In other words: I want to force some of the
cases to be in the same cluster from the start of the
-Original Message-
From: William Dunlap
Sent: Thursday, November 25, 2010 9:31 AM
To: 'randomcz'; r-help@r-project.org
Subject: RE: [R] help: program efficiency
If the input vector t is known to be ordered
(or if you only care about runs of duplicated
values, not all duplicated
Hello,
Somehow I cannot control my envelope plots on Spatstat package. I would like
to plot an envelope with no legend and also I would like to label the plot
as such:
plot(envelope(data, nsim=39, main=My K Plot, legend=FALSE,
xlab=Distance)
But somehow the main, legend and xlab calls do not
See if this works for you.
a - c(2,1,1,3,3,3,4)
a.fac - as.factor(a)
b - split(a, f = a.fac)
system.time(lapply(X = b, FUN = function(x) {
swn - seq(from = 0, to = 0 + 0.01*length(x), by = 0.01)
out - x + swn
return(out)
Oops, tiny mistake. Try
lapply(X = b, FUN = function(x) {
swn - seq(from = 0, to = (0 + 0.01*length(x))-0.01, by
= 0.01)
out - x + swn
return(out)
})
--
View this message in context:
see below.
On Fri, Nov 26, 2010 at 11:50 AM, Kjetil Halvorsen
kjetilbrinchmannhalvor...@gmail.com wrote:
see below.
On Wed, Oct 27, 2010 at 3:32 AM, Prof Brian Ripley
rip...@stats.ox.ac.uk wrote:
On Tue, 26 Oct 2010, Kjetil Halvorsen wrote:
I have been searching the documentation for some
Hi Terry, David, and Thomas,
Thank you for all your emails and the time you to took to clarify my
misunderstanding on survival analysis. I will need a bit of time to digest all
this information and to do some more reading.
Best regards,
Ben
From: Terry Therneau
1. survreg() does NOT
Hello,
Can we really make the assumption that the data is sorted. The original
example was not:
I am working on a function to make a duplicated value unique. For example,
the original vector would be like : a = c(2,1,1,3,3,3,4)
If we can make the assumption, here is a C++ based version:
Le 26/11/10 21:13, Romain Francois a écrit :
Hello,
Can we really make the assumption that the data is sorted. The original
example was not:
I am working on a function to make a duplicated value unique. For
example,
the original vector would be like : a = c(2,1,1,3,3,3,4)
If we can make
Can you create a small example that replicates your problem? If not,
at least send the code you are running right before R throws an error
(e.g., the exact call to lm) and some basic information on your data
like the class of each variable and the number of observations.
Is your outcome
On 2010-11-26 11:16, Neba Funwi-Gabga wrote:
Hello,
Somehow I cannot control my envelope plots on Spatstat package. I would like
to plot an envelope with no legend and also I would like to label the plot
as such:
plot(envelope(data, nsim=39, main=My K Plot, legend=FALSE,
xlab=Distance)
This
Date: Fri, 26 Nov 2010 11:25:26 -0800
From: roman.lust...@gmail.com
To: r-help@r-project.org
Subject: Re: [R] help: program efficiency
Oops, tiny mistake. Try
lapply(X = b, FUN = function(x) {
swn - seq(from = 0, to = (0 +
Dear expeRts,
I am not sure if I found a bug...
I would like to create a function that itself creates a lattice plot without
colors. Following
http://www.mail-archive.com/r-help@r-project.org/msg64699.html
I use trellis.device() to set the colors to FALSE. Whenever I call the minimal
example
Hi
I have 1800 response variables to regress on two factors (latitude and age),
what is the script to run all response variables at once instead of writing
1800 models?
Thanks
R.
__
R-help@r-project.org mailing list
On Nov 26, 2010, at 6:35 PM, Marius Hofert wrote:
Dear expeRts,
I am not sure if I found a bug...
I would like to create a function that itself creates a lattice plot
without
colors. Following http://www.mail-archive.com/r-help@r-project.org/msg64699.html
I use trellis.device() to set the
Hmm, for the sake of anyone who steps on the same mine. Basically, this class
error always say that you need to install latest java. It's not enough to
get new jre -- you need to get full latest jdk as well. After that you need
to do R reconfiguration, which is done by
R CMD javareconf
Then you
I want to change my parameter g to maximize the sum of my model. I keep
getting an unused argument error and I don't know why. Here are the details of
my problem.
g-0.2 #initial value for g
Qt-exp(-g*tagdat$t) #model building
PTT-Qt*Qt #model building
Thanks Fiona,
This is the right idea, but I would like to have the colors smoothed out
similar to what filled.contour does.
Jon
--
View this message in context:
http://r.789695.n4.nabble.com/Filled-contour-plot-showing-labeled-isolines-tp3056437p3060915.html
Sent from the R help mailing list
Thanks Steve. Could you please give me an example of what the command for
doing that would look like?
--
View this message in context:
http://r.789695.n4.nabble.com/Filled-contour-plot-showing-labeled-isolines-tp3056437p3060925.html
Sent from the R help mailing list archive at Nabble.com.
Thanks guys, the rle function works pretty well. Thank you all for the
efforts.
Zheng
--
View this message in context:
http://r.789695.n4.nabble.com/help-program-efficiency-tp3059079p3061103.html
Sent from the R help mailing list archive at Nabble.com.
?nlm
It does not seem to have a parameter 'data', so that is the cause of
your error message.
On Fri, Nov 26, 2010 at 4:36 PM, Mike Gibson megalop...@hotmail.com wrote:
I want to change my parameter g to maximize the sum of my model. I keep
getting an unused argument error and I don't know
Dear all,
I am fitting a glm to count data using poison errors with the log link. My
goal is to test for the significance of model terms by calling the anova
function on two nested models following the recommendation in Michael
Crawley's guide to Statistical Computing.
Without going into too
Hi folks
Run;
ToothGrowth
attach(ToothGrowth)
toothgrowth=lm(len~dose)
adding abline:
abline(toothgrowth)
I got it done adding single abline.
How to add more ablines on the same diagram?
I found following thread, applying mapply command;
Plotting multiple ablines
As Nick suggested, you're confusing the name of an object with the name(s)
of its elements.
Study this example:
lamda - 0.2
lamda
[1] 0.2
lamda - c(g=0.2)
lamda
g
0.2
lamda - c(1,3,4)
lamda
[1] 1 3 4
lamda - c(g=1, x=3, foo=4)
lamda
g x foo
1 3 4
lamda is now a
On Nov 26, 2010, at 9:30 PM, Jonathan Flowers wrote:
Dear all,
I am fitting a glm to count data using poison errors with the log
link. My
goal is to test for the significance of model terms by calling the
anova
function on two nested models following the recommendation in Michael
Hi David,
Thanks very much, that clears it up for me. I plan to report the result as
a typical F-test with numerator and denominator d.f., the value of F and the
significance. If you have other thoughts, I would appreciate it.
Thanks again.
Jonathan
On Fri, Nov 26, 2010 at 7:57 PM, David
Hi all,
It seems roxygen will create Rd files for any functions under the ./R/
directory, and I did not find a way to hide certain functions, even if
I specified a name space. E.g.
#' @export
f = function() {}
#' do nothing
g = function() {}
If I run roxygenize() on the package, both f and g
just add more 'abline' calls, or pass multiple values; e.g.,
abline(v=c(1,2,3,4), h=c(10,20))
On Fri, Nov 26, 2010 at 9:41 PM, Stephen Liu sati...@yahoo.com wrote:
Hi folks
Run;
ToothGrowth
attach(ToothGrowth)
toothgrowth=lm(len~dose)
adding abline:
abline(toothgrowth)
I got it done
Hi jim,
Tks for your advice. I got it.
1)
abline(v=c(1,2,3,4), h=c(10,20), toothgrowth)
The output of following command is the same.
abline(v=c(1,2), h=c(10,20), toothgrowth)
2)
Afterwards can I add further abline?
B.R.
Stephen L
From: jim holtman
thanks, I succeed.
kevin
--
View this message in context:
http://r.789695.n4.nabble.com/how-to-get-the-plot-like-the-attachment-tp3060425p3061217.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
Thank you, Joshua -- great name by the way; he was a Biblical champ! -- fro
your suggestion. Sorry about the belated response. Nope, even the script
editor will not allow pastes (after the R Console has also refused to allow
pastes). I do agree that it is probably not R but W XP.
Cheers
Troy
--
89 matches
Mail list logo