useRs,
Is there a way to schedule R scripts? I would like to run certain
scripts three times a day. I'm running R on Windows XP.
sessionInfo()
R version 2.6.0 (2007-10-03)
i386-pc-mingw32
locale:
Thanks for your detailed explanation.
You are right, a set of boxplots done with bwplot
is a much better graphic for this type of data:
bwplot(V1~VAR|f,data=datos2)
This was not a good example. The barplot would be suited
for counts, ie. species composition:
datos4 -
glenn andrews wrote:
What should I be looking for in the output of the nls() routine that
alerts me to the fact that the Hessian is potentially ill-conditioned?
Unreasonably large s.e.'s and large correlations in the
variance-covariance matrix of estimates (cov2cor(vcov(nlmod)) or
On Thu, 27 Mar 2008, Peter Dalgaard wrote:
glenn andrews wrote:
Thanks for the response. I was not very clear in my original request.
What I am asking is if in a non-linear estimation problem using nls(),
as the condition number of the Hessian matrix becomes larger, will the
t-values of one
Adrian Johnson wrote:
Hello:
I have a tab delim file with 100 rows and 100 columns.
I have numerical values in this table. What I want is to create an image
color map with color gradation.
my values range from -5 to 0. max value is 0.
to acheive this, I need to convert my data.frame
Under WIndows just write a .bat file that is called by the Scheduled
Tasks to execute when you want it.
On Fri, Mar 28, 2008 at 1:54 AM, Lauri Nikkinen [EMAIL PROTECTED] wrote:
useRs,
Is there a way to schedule R scripts? I would like to run certain
scripts three times a day. I'm running R on
Hello there,
I have two groups (men and women) and I know per group how many of them
smoke or don't smoke (women 40 of 200; men 100 of 300). I would like to know
how I can compare in R if men and women differ significantly in their
smoking. However, because there are more men in the sample than
On Wed, Mar 26, 2008 at 05:27:22PM +1300, Frank Scherr wrote:
Hello together,
Is there a tool to test the statistical differences between parameter
estimates of a nlsList fit?
I fitted degradation data using the nlsList method and want to find out
whether derived rate constants are
Lila86 wrote:
I have two groups (men and women) and I know per group how many of them
smoke or don't smoke (women 40 of 200; men 100 of 300). I would like to
know how I can compare in R if men and women differ significantly in their
smoking. However, because there are more men in the
Hi
[EMAIL PROTECTED] napsal dne 28.03.2008 10:18:20:
Adrian Johnson wrote:
Hello:
I have a tab delim file with 100 rows and 100 columns.
I have numerical values in this table. What I want is to create an
image
Probably not. See below.
color map with color gradation.
my
Okay, I will do that, thank you.
prop.test seems to work, too, someone wrote me?!
prop.test(x=c(40,100), n=c(200,300))
I tried it and it seems fine.. is it the same or is chisq.test better? I
will read a bit about chi-square tests now.. Thanks a lot,
Lila
Richard Cotton wrote:
Lila86
On Mar 28, 2008, at 2:42 AM, Agustin Lobo wrote:
Thanks for your detailed explanation.
You are right, a set of boxplots done with bwplot
is a much better graphic for this type of data:
bwplot(V1~VAR|f,data=datos2)
This was not a good example. The barplot would be suited
for counts, ie.
On Fri, Mar 28, 2008 at 1:42 AM, Agustin Lobo [EMAIL PROTECTED] wrote:
Thanks for your detailed explanation.
You are right, a set of boxplots done with bwplot
is a much better graphic for this type of data:
bwplot(V1~VAR|f,data=datos2)
This was not a good example. The barplot would be
Why in arima output there is no significance tests of estimated parameters as
in regression models? I know that I have s.e. and I could build these test by
myself but are there any reasonable reasons?
Best,
Rob
__
R-help@r-project.org mailing list
Use nlme. Here is an example, fitting the Michaelis-Menten model to some enzyme
data:
library(nlme)
data(Puromycin)
fit.nls - nlsList(rate ~ SSmicmen(conc, Vm, K)|state, data=Puromycin)
summary(fit.nls)
# Use the output to calculate starting values for the nlme fit.
fit.nlme - nlme(rate ~
Is there any other way I can do this.
say row names are city schools
colnames are student names
Hi Adrian,
Have you looked at color2d.matplot in the plotrix package?
Jim
You might also want to check out the seriation package, which provides
tools to reorder your matrix to
Dear R-Users,
I will be really grateful if you help me here. This is my script. I want to
extract the vif values for x1 and x2 respectively as well as the regression
coefficients for x1 and x2. Because I then want to plug them in a next step
of the analysis. Thank you in advance! I am using
http://www.nabble.com/file/p16349649/at-length.png
Hi,
I produced this graph with R and I am very happy with it. But now I wonder:
Is this some standard way of displaying data and if so, does it have a name?
(like Histogram, Box plot etc.)
I need to explain the plot to people so knowing a bit
On 3/28/2008 8:37 AM, mika03 wrote:
http://www.nabble.com/file/p16349649/at-length.png
Hi,
I produced this graph with R and I am very happy with it. But now I wonder:
Is this some standard way of displaying data and if so, does it have a name?
(like Histogram, Box plot etc.)
It's a
Hi Spencer and David,
Thanks for your answers first - yes it is deviance but just before i just
spoke and explain that it is the equivalent of r square from the normal
regression.
I hope i can do the comparison and show that the model is significant and
hopefully i am off the
Hi,
I think it looks like a spineplot.
Best regards,
Peter.
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of mika03
Sent: 28. marts 2008 13:37
To: r-help@r-project.org
Subject: [R] Quick question: Does this graph have a name?
I updated my modified GNU Emacs for Windows release to the now current
version 22.2 of Emacs. The distribution is available at
http://vgoulet.act.ulaval.ca/en/emacs
For those hearing about this distribution for the first time, here's a
short description: ready-to-install version of
I have always heard this referred to as a stack bar graph.
Ken
mika03 wrote:
http://www.nabble.com/file/p16349649/at-length.png
Hi,
I produced this graph with R and I am very happy with it. But now I
wonder: Is this some standard way of displaying data and if so, does it
have a
Post-Doctoral Fellowship
Modelling of Cardiovascular Safety (full time)
Leiden University, Leiden, The Netherlands
javascript:doredirect();
Job description
TOP Institute Pharma (TI Pharma) has granted our proposal to set up a
mechanism-based PKPD modelling platform. This
Principal Post-Doctoral Fellow
PKPD modelling platform (full time)
Leiden University, Leiden, The Netherlands
javascript:doredirect();
Job description
TOP Institute Pharma (TI Pharma) has granted our proposal to set
I think it looks like a spineplot.
Which is a special case of a mosiac plot with only two variables.
Hadley
PS. I suspect that the plot could be improved by ordering the x and y
variables by their marginal counts. (Although maybe the y axis
already has an ordering?)
--
http://had.co.nz/
I have tried the following code to plot and save many box plots in files. The
code works by itself, but does not run in a loop. The error indicates that
there is no device. The error is printed below the code.
Code:
for(i in 1:len)
{
A = as.numeric(delta[Delta[i,1]])
B =
A new version of the Matrix
(well, actually the R package named Matrix) has become
available on the CRAN mirrors.
As some of you have noticed, the version numbers (current is version
0.999375-8) are converging to one, and we feel that we have solved
enough of the many (mostly small) problems
Hi All,
Just hoping some one can give me a hand with a problem...
I have a dataframe (DF) with about 5 million entries that looks something
like the following:
DF
ID Cl Co BrdInd A AB AB
1 S-3 IND A BR_F BR_F01 1 0 0
2 S-3 IND A BR_F BR_F01 1 0 0
3 S-3 IND A BR_F BR_F01 1
Hi,
this is probably quite stupid but I have no clue
what's wrong. Let's say I write the function
hline - function() {
cat(\\hline \n)
}
and call hline() from within a Sweave chunk. Why is
there no carriage return after the \hline in the
resulting tex file?
if I call hline() hline() in the
mika03 wrote:
http://www.nabble.com/file/p16349649/at-length.png
Hi,
I produced this graph with R and I am very happy with it. But now I wonder:
Is this some standard way of displaying data and if so, does it have a name?
(like Histogram, Box plot etc.)
I need to explain the plot to
Do you think p values are bad? That's not my understanding. P
values may not be reported by some software, because the algorithm
developers didn't know how to efficiently compute a reasonably accurate
p value. And if you do a thousand or a million statistical tests and
report only the
This is an instance of FAQ Q7.22: you haven't plotted anything and so the
device is not yet open.
On Fri, 28 Mar 2008, Todd Remund wrote:
I have tried the following code to plot and save many box plots in files.
The code works by itself, but does not run in a loop. The error indicates
Hi Sean,
is this roughly what you are looking for (please note that in the
example data you provided there is only one level of ID given, no S-4,
...) ?
DF
ID Cl Co BrdInd A AB AB.1 frq
1 S-3 IND A BR_F BR_F01 1 00 1.0
2 S-3 IND A BR_F BR_F01 1 00 1.0
3 S-3 IND A
Is there a way of changing the line type for the wiskers in bwplot? I have
been able to do pretty much everything I would like to do in bwplot, just can't
change these lines from dashed to solid. Thanks.
Todd Remund
[[alternative HTML version deleted]]
On Fri, Mar 28, 2008 at 9:23 AM, Todd Remund [EMAIL PROTECTED] wrote:
Is there a way of changing the line type for the wiskers in bwplot? I have
been able to
do pretty much everything I would like to do in bwplot, just can't change
these lines from
dashed to solid. Thanks.
These are
Hi,
1)
par(ask=TRUE)
histogram(~Y|F1 * F2 * F3,data = ...)
Successively displays many plots as the # of levels in F3 w/o asking for
input between each.
In other words I only get to see the plot for the last level. What's the
right commands?
2)
Also I am wondering if there are plot objects i.e.
On 28/03/2008, [EMAIL PROTECTED]
[EMAIL PROTECTED] wrote:
Hi,
1)
par(ask=TRUE)
histogram(~Y|F1 * F2 * F3,data = ...)
Successively displays many plots as the # of levels in F3 w/o asking for
input between each.
In other words I only get to see the plot for the last level. What's the
On Fri, 28 Mar 2008, [EMAIL PROTECTED] wrote:
Hi,
1)
par(ask=TRUE)
histogram(~Y|F1 * F2 * F3,data = ...)
Successively displays many plots as the # of levels in F3 w/o asking for
input between each.
In other words I only get to see the plot for the last level. What's the
right commands?
Dear R community,I have a scatterplot with multiple vertical ablines. I
wish to color each interval between two ablines in a different color...
Could you please indicate me how to do this efficiently?
Thank you!
Georg.
Georg Ehret
Johns Hopkins
Baltimore
Dear all,
I have prepared a new package to install in R, I tried to check my package
and I got the following message
232sub141:~ zahra$ R CMD check /Users/zahra/Desktop/HighProbability
* checking for working latex ... OK
*
Hello,
I have read recent posts on this topic (Dr. Ronnen Levinson's Monday 02:39:55
pm), but before I install a 64 bit system, and purchase more RAM, I want to
make sure I understand this interesting issue.
I was attempting to plot a data frame containing Dow Jones stock information:
Looks like you attach() the data frame before you try t/o plot. Note
that in the Details section of ?attach, it says:
... The database is not actually attached. Rather, a new environment is
created on the search path and the elements of a list (including columns
of a data frame) or objects in a
//Referring to the response posted many years ago, copied below, what
is the specific criterium used for singularity of the gradient matrix?
Is a Singular Value Decomposition used to determine the singular
values? Is it the gradient matrix condition number or some other
criterion for
On 3/27/08, Karin Lagesen [EMAIL PROTECTED] wrote:
Hello.
I have made many normal boxplots where I have added a new boxplot to
an existing one. When I have done this, I have used the at command to
move the boxplots a bit so that they could fit next to eachother, like
this:
On 3/27/08, Weidong Gu [EMAIL PROTECTED] wrote:
I have two data sets with locations, X, Y of houses (df.house) and
habitats(df.habitat), respectively. In each dataset, there are 3
replicates (Repeat). Because each replicate has different locations of
houses and habitats, I would like to
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf
Of John
Sent: Friday, March 28, 2008 12:04 PM
To: r-help@r-project.org
Subject: [R] Error: cannot allocate vector of size 3.0 Gb
Hello,
I have read recent posts on this topic (Dr. Ronnen Levinson's
you haven't escaped the \ for the \n, I think. Your line should be
cat(\\hline \\n). You did escape the \ for hline, though.
Abhijit Dasgupta, Ph.D
Assistant Professor | Division of Biostatistics
Dept of Pharmacology and Experimental Therapeutics | Thomas Jefferson
University
1015 Chestnut St |
A QR decomposition is done on the weighted gradient matrix; if the
estimate of the rank that results is less than the number of columns in
the gradient (the number of nonlinear parameters), or less than the number
of rows (the number of observations), nls stops.
You can see the calls in the
Georg Ehret [EMAIL PROTECTED] wrote:
Dear R community,I have a scatterplot with multiple vertical ablines. I
wish to color each interval between two ablines in a different color...
Could you please indicate me how to do this efficiently?
Thank you!
Georg.
Dear
On Fri, Mar 28, 2008 at 7:52 PM, [EMAIL PROTECTED] wrote:
* checking whether the package can be loaded ... ERROR
Error in library(HighProbability) :
'HighProbability' is not a valid package -- installed 2.0.0?
Execution halted
Although it might be obvious, can you load this package
At the risk of oversimplifying, there are three steps to developing
and using a package:
1) build
2) install
3) load
Evidently your package does steps 1 and 2 but not 3.
Furthermore, it tells you, not a valid package, and says something
that is suggestive of a version problem:
installed
I have a sparse contingency table (most cells are 0):
xtabs(~.,data[,idx:(idx+4)])
, , x3 = 1, x4 = 1, x5 = 1
x2
x11 2 3
1 0 0 31
2 0 0 112
3 0 0 94
, , x3 = 2, x4 = 1, x5 = 1
x2
x11 2 3
1 0 0 0
2 0 0 0
3 0 0 0
, , x3 = 3, x4
R users,
I've been searching around for a solution to this for a few days, any help
would be appreciated.
I have working R code (a custom function) that transforms 3D coordinate data
into angles (between any three sets of coordinates). The code outputs a
single scalar (an angle). For my
Dear 'Born',
There was thread on this recently, but I cannot seem to find it.
The best suggestion (IMHO) was along these lines:
aggregate( rep(1,40), as.data.frame(diag(4)[sample(1:4,40,repl=T),]), sum )
See also
http://thread.gmane.org/gmane.comp.lang.r.general/104798/focus=104841
and
55 matches
Mail list logo