Hi,
I want to control the plots in the output of the xyplot(). It is
easier to explain it through an example:
#-
library(lattice);
# months
months - c(Jan, Feb, Mar,
Apr, May, Jun,
Jul, Aug, Sep,
Oct,
Ravi Varadhan rvaradhan at jhmi.edu writes:
Try the package Rdonlp2, which can handle general, nonlinear, equality and
inequality constraints for
smooth optimization problems.
But be aware of the licencing.
Dieter
__
R-help@r-project.org mailing
Hi list,
has anyone succeeded in using pgmm() on any dataset besides
Arellano/Bond's EmplUK, as shown in the vignette?
Whatever I try, I eventually get a runtime error because of a singular
matrix at various points in pgmm.diff() (which gets called by pgmm()).
For example, when estimating a
Hi Anand,
The output of the xyplot() is sorted by the alphabetical order of
months (First Apr, then Aug so on). I need to output such that the
order of the months is maintained, i.e. plot of Jan first, then Feb, Mar
... so on.
Because the levels of your factor are sorted this way:
AR anandram at gmail.com writes:
I want to control the plots in the output of the xyplot(). It is
easier to explain it through an example:
.. modified example see below
The output of the xyplot() is sorted by the alphabetical order of
months (First Apr, then Aug so on). I need to output
Steven Lubitz slubitz1 at yahoo.com writes:
Thank you - this is very helpful. However I realized that with my real data
sets (not the example I have here),
I also have different numbers of columns in each data frame. rbind doesn't
seem to like this. Here's a
modified example:
x -
On Sat, 07 Mar 2009 18:04:24 -0500, David Winsemius wrote :
[ Snip ... ]
Nonetheless, I do think the relative numbers of package downloads might
be interpretable, or at the very least, the basis for discussions over
beer.
*Anything* might be the basis for discussions over beer (obvious
I do not understand, from a PCA point of view, the option center=F
of prcomp()
According to the help page, the calculation in prcomp() is done by a
singular value decomposition of the (centered and possibly scaled) data
matrix, not by using eigen on the covariance matrix (as it's done by
At 08:47 05/03/2009, herrdittm...@yahoo.co.uk wrote:
Dear useRs,
With the rise of netbooks and 'lifestyle laptops I am tempted to
get one of these to mainly run R on it. Processor power and hard
disk space seem to be ok. What I wonder is the handling and feel
with respect to R.
Has anyone
Dear list,
I'd like to ask for some advice in creating a wrapper function for
this C code,
http://www.fourmilab.ch/documents/specrend/
http://www.fourmilab.ch/documents/specrend/specrend.c
I could probably port it in R, but I've been hoping to use compiled
code for a while and this looks
Dear R users,
Is it possible to run a R script under a python script ?
And if yes ... how to do that ?
Thanks a lot,
Arnaud
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
Google for RPy2 and RSpython.
On Sun, Mar 8, 2009 at 8:38 AM, arnaud_mosn...@uqar.qc.ca wrote:
Dear R users,
Is it possible to run a R script under a python script ?
And if yes ... how to do that ?
Thanks a lot,
Arnaud
[[alternative HTML version deleted]]
Hi Agus,
But the rotation made with the eigenvectors of prcomp(X,center=F) yields
axes that are correlated. Therefore, prcomp(X,center=F) is not really a
PCA.
cor() is not an appropriate test of whether two vectors are orthogonal. The
definition that two vectors (in an inner product space)
Thanks for all the answers, they were really helpful. I noticed the
inflation of the variance, and that seems to be the problem indeed (thank
you, Thomas). When looking at the tables again (thank you Ted), I noticed
that this is actually a classical case of (quasi-)complete separation. Hence
the
Dear R-experts,
I have the following long data frame:
ID -
c(ID001,ID001,ID001,ID001,ID001,ID001,ID001,ID001,ID001,ID001,ID001,ID002,ID002,ID002,ID002,ID002,ID002,ID002)
Var -
Hi all,
i have a SpatialPixelsDataFrame object, which is a result of kriging in R.
What i would like to do is to convert this object to image (geoTiff) so that
i could use it in Gis application. I have tried using this image2Grid
function but it doesen't work. Here is the code :
im-image(pCir)
More seriously : I don't think relative numbers of package downloads can
be interpreted in any reasonable way, because reasons for package
download have a very wide range from curiosity (what's this ?), fun
(think fortunes...), to vital need tthink lme4 if/when a consensus on
denominator DFs
On Sun, Mar 8, 2009 at 10:49 AM, hadley wickham h.wick...@gmail.com wrote:
More seriously : I don't think relative numbers of package downloads can
be interpreted in any reasonable way, because reasons for package
download have a very wide range from curiosity (what's this ?), fun
(think
Dear Agustin the Listers,
Noncentred PCA is an old and establishes method. It is rarely used,
but still (methinks) it is used more often than it should be used.
There is nothing wrong in having noncentred PCA in R, and it is a real
PCA. Details will follow.
On 08/03/2009, at 11:07 AM,
On 08/03/2009 10:49 AM, hadley wickham wrote:
More seriously : I don't think relative numbers of package downloads can
be interpreted in any reasonable way, because reasons for package
download have a very wide range from curiosity (what's this ?), fun
(think fortunes...), to vital need tthink
Dear all,
Have a complete document about the Batch Execution of R
or nice website?
--
Nash - morri...@ibms.sinica.edu.tw
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Dear R-user,
I'm running a certain R script in DOS batch mode. Is there a way to prevent
R from saving the workspace once this script is finished? I'm asking this
because the resulting .RData file has the size of 70MB. I don't need this
file since my script already writes the required output and
Hi,
I have been using this site (
http://www.stat.pitt.edu/stoffer/tsa2/Rissues.htm) to help me with some
ARIMA modelling in R.
Unfortunately the methods mentioned do not appear to work with second order
differencing; arima(*, 2, *).
I have used some dummy data to illustrate my point.
When I
?quit # ,save=no
--
DW
On Mar 8, 2009, at 11:19 AM, Jim_S wrote:
Dear R-user,
I'm running a certain R script in DOS batch mode. Is there a way to
prevent
R from saving the workspace once this script is finished? I'm asking
this
because the resulting .RData file has the size of
http://stat.ethz.ch/R-manual/R-patched/library/utils/html/BATCH.html
HTH, Jim
Nash-16 wrote:
Dear all,
Have a complete document about the Batch Execution of R
or nice website?
--
Nash - morri...@ibms.sinica.edu.tw
__
Hi:
I am looking for style guides for larger R programs,
and would like suggestions of web pages or documents to look at.
I am particularly interested in guidelines for how to structure programs,
and such issues as managing scope, passing data around, data hiding,
and so forth.
I do have some
This works fine.
Thank you for the quick response,
Jim
David Winsemius wrote:
?quit # ,save=no
--
DW
On Mar 8, 2009, at 11:19 AM, Jim_S wrote:
Dear R-user,
I'm running a certain R script in DOS batch mode. Is there a way to
prevent
R from saving the workspace once this
This works fine.
Thank you for the quick response,
Jim
David Winsemius wrote:
?quit # ,save=no
--
DW
On Mar 8, 2009, at 11:19 AM, Jim_S wrote:
Dear R-user,
I'm running a certain R script in DOS batch mode. Is there a way to
prevent
R from saving the workspace once this
I think the situation is worse than messy. If a client comes in with data
that doesn't address the question they're interested in, I think they are
better served to be told that, than to be given an answer that is not
actually valid. They should also be told how to design a study that
On 08-Mar-09 15:14:03, Duncan Murdoch wrote:
On 08/03/2009 10:49 AM, hadley wickham wrote:
More seriously : I don't think relative numbers of package downloads
can be interpreted in any reasonable way, because reasons for
package download have a very wide range from curiosity (what's
this ?),
hi, all
I have two questions on using R console effectively (this is on Mac,
not sure if it applies to win platform):
First, I'd like to make the console font bigger, the default is too
small for my screen. There is a Show Fonts from Format menu where
you can adjust it, but it seems only for
Hi Ted,
Coming to think about your direction - another idea came to mind:
The next time a major release is made (there is one scheduled quite soon
actually), the core team could add a survey on the downloading page of the
R base package asking for just one question
please click here if this is
Is this another discussion of what data might be collected and
analyzed, and what could and could not be said if we only had such data?
Has anyone but me produced any actual data? If so, I missed it.
Hadly mentioned the 'fortunes' package. My earlier methodology,
hi,
the code below produces me two curved lines.
Now I want to fill the space between these curves with a color and I
don't have any idea how to do that with R.
I tried it with the polygon() function but it didn't worked out.
Can someone her give me a hint on that?
thanks!
x1 -
Hi,
Try:
polygon(c(x1, rev(x2)), c(y1, rev(y2)), col=grey)
In general you might need to make sure the data is well ordered.
Hope this helps,
baptiste
On 8 Mar 2009, at 16:52, Martin Batholdy wrote:
hi,
the code below produces me two curved lines.
Now I want to fill the space
A dataframe holds 3 vars, each checked true or false (1, 0). Another
var holds the grouping, r and s:
### start:example
set.seed(20)
d - data.frame(sample(c(0, 1), 20, replace=T), sample(c(0, 1), 20,
replace=T), sample(c(0, 1), 20, replace=T))
names(d) - c(A, B, C)
e - rep(c(r, s), 10)
###
Hi,
You could use the reshape package:
d$e - e
recast(d, variable~e, fun=sum)
The doBy package is another option.
baptiste
On 8 Mar 2009, at 17:14, soeren.vo...@eawag.ch wrote:
A dataframe holds 3 vars, each checked true or false (1, 0). Another
var holds the grouping, r and s:
###
On 08/03/2009 12:08 PM, Barry Rowlingson wrote:
I think the situation is worse than messy. If a client comes in with data
that doesn't address the question they're interested in, I think they are
better served to be told that, than to be given an answer that is not
actually valid. They should
On Sun, Mar 8, 2009 at 7:08 AM, Michael Dewey i...@aghmed.fsnet.co.uk wrote:
At 08:47 05/03/2009, herrdittm...@yahoo.co.uk wrote:
Dear useRs,
With the rise of netbooks and 'lifestyle laptops I am tempted to get one
of these to mainly run R on it. Processor power and hard disk space seem to
Thanks for the help, but ANOVA will give me a single pvalue, then how i can
make sure which group is showing the significant differences.
Tal Galili wrote:
Hello friend.
I believe anova might be a better solution for you.
You might have a look here:
Hey all,
I am trying to use the survreg function in R to estimate the mean and
standard deviation to come up with the MLE of alpha and lambda for the
weibull distribution. I am doing the following:
times-c(10,13,18,19,23,30,36,38,54,56,59,75,93,97,104,107,107,107)
Dear Barry,
As far as I understand, you're telling us that having a bit of data
mining does not harm whatever the data. Your example of pop music charts
might support your point (although my ears disagree ...) but I think it
is bad policy to indulge in white-noise analysis without a well-reasoned
Robert Biddle schrieb:
Hi:
I am looking for style guides for larger R programs,
and would like suggestions of web pages or documents to look at.
I am particularly interested in guidelines for how to structure programs,
and such issues as managing scope, passing data around, data hiding,
and so
David1234 wrote:
Hi everyone,
Im quite new to R an I have the following Question:
I would like to define Variables which I can add and multiply etc. and that
R can simplyfy the terms.
The variables should stand for integers.
For example I would like to have an entry in an array with
On 08-Mar-09 17:44:18, Douglas Bates wrote:
On Sun, Mar 8, 2009 at 7:08 AM, Michael Dewey i...@aghmed.fsnet.co.uk
wrote:
At 08:47 05/03/2009, herrdittm...@yahoo.co.uk wrote:
Dear useRs,
With the rise of netbooks and 'lifestyle laptops I am tempted
to get one of these to mainly run R on it.
On 8 March 2009 at 13:27, Duncan Murdoch wrote:
| But we don't even have that data, since CRAN is distributed across lots
| of mirrors.
On 8 March 2009 at 19:01, Emmanuel Charpentier wrote:
| As far as I can see (but I might be nearsighted), I see no model linking
| package download to package
1Rnwb schrieb:
Thanks for the help, but ANOVA will give me a single pvalue, then how i can
make sure which group is showing the significant differences.
Hi,
ANOVA is fine and please have a look on ?TukeyHSD and don't forget to
consult your statistics textbook about post-hoc testing. If you
I think R and survival took a look at what you provided for a formula
to survreg and rejected it because it had no lh-rh separator, ~, and
thus no rh terms. Turns out that it also expects a dataframe even if
both terms for the Surv function are already given and there is
nothing further to
If you prefer to do it with base functions, you could start by
thinking about one column at a time, for which the tapply function is
a logical choice:
tapply(d[,A], e, sum)
r s
3 6
Then wrap that in an apply call that handles each column sequentially:
apply(d, 2, function (x) tapply(x, e,
since the estimation is not done pairwise so i cannot use pairwise.t.test,
how do i apply tukeyHSD
Thomas Petzoldt-4 wrote:
1Rnwb schrieb:
Thanks for the help, but ANOVA will give me a single pvalue, then how i
can
make sure which group is showing the significant differences.
Hi,
1Rnwb wrote:
since the estimation is not done pairwise so i cannot use pairwise.t.test,
how do i apply tukeyHSD
Note correct capitalization: TukeyHSD and follow the examples on the
help page:
?TukeyHSD
You may also contact a statistics textbook. In addition I send you some
links off list.
Dirk Eddelbuettel wrote:
On 8 March 2009 at 13:27, Duncan Murdoch wrote:
| But we don't even have that data, since CRAN is distributed across lots
| of mirrors.
On 8 March 2009 at 19:01, Emmanuel Charpentier wrote:
| As far as I can see (but I might be nearsighted), I see no model linking
|
Suppose a house is cleaned by three people in succession. The time it takes
each cleaner in hours is an Exp(0.5). Assume that the time required for each
cleaner is independent. How would I use simulation in R to estimate the
probability that it takes longer than a total of 5 hours for the house
On 9/03/2009, at 4:14 AM, Duncan Murdoch wrote:
... analyzing bad data will just give bad conclusions.
Fortune?
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and
On 8/03/2009, at 3:54 AM, David Freedman wrote:
Hi - I'd like to construct and plot the percents by year in a small
data set
(d) that has values between 1988 and 2007. I'd like to have a
breakpoint
(buy no discontinuity) at 1996. Is there a better way to do this
than in
the code
On 9/03/2009, at 8:42 AM, BowlesMcCann wrote:
Suppose a house is cleaned by three people in succession. The time
it takes
each cleaner in hours is an Exp(0.5). Assume that the time required
for each
cleaner is independent. How would I use simulation in R to estimate
the
probability
Hi R users,
Can anyone share some example code using merge_all (from the reshape
package) to merge 10 data frames into 1 file.
Thanks in advance for any help!
--
View this message in context:
Hi,
Try this:
DF1 - data.frame(var1 = letters[1:5], x = rnorm(5), y =2)
DF2 - data.frame(var1 = letters[3:7], x = rnorm(5), y=3)
DF3 - data.frame(var1 = letters[6:10], x = rnorm(5), y=0)
# ... DF10 if you wish
( result - merge_all(list(DF1, DF2, DF3) ))
save( result, file =merged.rda)
I
Thanks for posting this. I've never used Eclipse before but this
document inspired me to give it a try. Unfortunately, it's a little
out of date, but I think that I finally got it to work. In case
anyone else would like to try this, I put a PDF of my notes online at:
Rolf Turner wrote:
On 9/03/2009, at 4:14 AM, Duncan Murdoch wrote:
... analyzing bad data will just give bad conclusions.
Fortune?
looking for fortunes? got one for you:
A key reason that R is a good thing is because it is a language
who/where is left as an (easy) exercise.
vQ
On 9/03/2009, at 10:23 AM, John Fox wrote:
Dear Rolf,
Tukey put it nicely: The combination of some data and an aching
desire for
an answer does not ensure that a reasonable answer can be extracted
from a
given body of data. Inasmuch as there are no current fortunes from
Tukey, I
Le dimanche 08 mars 2009 à 13:22 -0500, Dirk Eddelbuettel a écrit :
On 8 March 2009 at 13:27, Duncan Murdoch wrote:
| But we don't even have that data, since CRAN is distributed across lots
| of mirrors.
On 8 March 2009 at 19:01, Emmanuel Charpentier wrote:
| As far as I can see (but I
On 8 March 2009 at 23:45, Emmanuel Charpentier wrote:
| Le dimanche 08 mars 2009 13:22 -0500, Dirk Eddelbuettel a crit :
| Once you have data, you have an option of using or discarding it. But if you
| have no data, you have no option. How is that better?
|
| I question 1) the usefulness
2009/3/8 Emmanuel Charpentier charp...@bacbuc.dyndns.org:
I question 1) the usefulness of the effort necessary to get the data ;
and 2) the very concept of data mining, which seems to be the rationale
for this proposed effort.
Furthermore (but this is seriously off-topic), I seriously
I tried using merge_all as shown below but I am getting an error ... can
anyone tell me what I am doing wrong? The result table below is what I am
looking for.
DF1 - data.frame(var1 = letters[1:5], x = rnorm(5), y =2)
DF2 - data.frame(var1 = letters[1:5], t = rnorm(5), u =2)
DF3 -
One can just include it in the regression, but the potential problems
for interpretation are surely greater than those indicated. Inclusion
of
X1 = T1+E1 may cause X2 to appear significant when in fact it is having
no effect at all. Or the true effect can be reversed in sign. This
happens
The function expects a list of data.frames as a first argument but you
provided a data.frame instead, the others are interpreted as optional
arguments to merge_recurse(). Try this instead,
merge_recurse(list(DF1,DF2,DF3,DF4))
var1x yt u d ef o
1a
Perfect - many thanks!!!
baptiste auguie-2 wrote:
The function expects a list of data.frames as a first argument but you
provided a data.frame instead, the others are interpreted as optional
arguments to merge_recurse(). Try this instead,
merge_recurse(list(DF1,DF2,DF3,DF4))
Oliver,
Go and find the file named 'Rconsole' under ~/etc folder, then you can
change whatever you want, the font size, color etc. The settings will be
your default.
For your second question, you need to set it up in Rprofile.site. Refer to
the Rprofile help.
Jun
On Sun, Mar 8, 2009 at 11:20
I question 1) the usefulness of the effort necessary to get the data ;
and 2) the very concept of data mining, which seems to be the rationale
for this proposed effort.
Furthermore (but this is seriously off-topic), I seriously despise the
very idea of popularity in scientific debates...
hi,
this is my first time I use the sem package in R.
I made a simple path analysis.
Now I was wondering how to get the standardized solution.
How can I get the standardized estimates of the path coefficients?
__
R-help@r-project.org mailing list
Dear All,
I have a csv file which has total 510 records and 15 variables.
I am trying to write a code which will return me for each record which
variable is missing.
For e.g.
If for record 1, variables 8 9 are missing, the R code should return me
that.
In similar way I need to find the missing
hi,
I have following data and code;
cov -
c
(1.670028
,-1.197685
,-2.931445,-1.197685,1.765646,3.883839,-2.931445,3.883839,12.050816)
cov.matrix - matrix(cov, 3, 3, dimnames=list(c(y1,x1,x2),
c(y1,x1,x2)))
path.model - specify.model()
x1 - y1,x1-y1
x2 - x1, x2-x1
At 4:16 AM +0100 3/9/09, Martin Batholdy wrote:
hi,
this is my first time I use the sem package in R.
I made a simple path analysis.
Now I was wondering how to get the standardized solution.
How can I get the standardized estimates of the path coefficients?
?std.coef
Martin,
hi,
I have following data and code;
cov -
c
(1.670028
,-1.197685
,-2.931445,-1.197685,1.765646,3.883839,-2.931445,3.883839,12.050816)
cov.matrix - matrix(cov, 3, 3, dimnames=list(c(y1,x1,x2),
c(y1,x1,x2)))
path.model - specify.model()
x1 - y1, x1-y1
x2 - x1, x2-x1
x2
thank you very much!
I definitely need more theoretical background ...
but for now;
what does that mean for this dataset?
x1 should be the intermediate variable of x2 and y1
(x2 - x1 - y1)
Can I test that with this kind of analysis?
or do I see know that this kind of intermediate variable
If your data.frame after importing that csv file is DF, then this will
return a list of variables with the row numbers that have missing
values for each variables:
sapply(DF, function(x) which(is.na(x)))
But the overall tenor of you question suggests that you may need to
back up a step
Dear R Gurus:
What is the difference between a Primitive and a Generic, please?
Thanks,
Edna Bell
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
78 matches
Mail list logo
