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
You can do it interactively with a script. If you know how large your
output will be, you can preallocate an object and then just store the
results in it. The loop would look something like this:
files <- Sys.glob("*.csv") # get names of files to process
result <- numeric(length(files)) # prea
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 research
Wells, Brian wrote:
> I have the following code.
>
>
>
>
>
>> f <- cph(formula = Surv(TimeToDeath, Dead == "Yes")
> ~1,data=single.dat, x=T, y=T, surv=T)
>
>> for(i in c('A', 'B', 'C', 'D', 'E', 'F')){
>
>> f <-update(f,as.formula(paste('Surv(TimeToDeath, Dead ==
> "Yes")~',i,sep='')))
>
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
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
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 packa
"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.
> ***
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 source
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 |
> -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 Lev
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 t
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:
>
//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 determ
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
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:
> length
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
* usi
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
[[alterna
I have the following code.
> f <- cph(formula = Surv(TimeToDeath, Dead == "Yes")
~1,data=single.dat, x=T, y=T, surv=T)
> for(i in c('A', 'B', 'C', 'D', 'E', 'F')){
> f <-update(f,as.formula(paste('Surv(TimeToDeath, Dead ==
"Yes")~',i,sep='')))
> print(summary(f, paste(i,"=1st Quartile",
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 c
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. W
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 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.
The
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]]
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
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
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
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 p
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
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
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) proble
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 = as.numeric(
> 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/
Principal Post-Doctoral Fellow
PKPD modelling platform (full time)
Leiden University, Leiden, The Netherlands
Job description
TOP Institute Pharma (TI Pharma) has granted our proposal to set up a
mechanism-based PK-PD
Post-Doctoral Fellowship
Modelling of Cardiovascular Safety (full time)
Leiden University, Leiden, The Netherlands
Job description
TOP Institute Pharma (TI Pharma) has granted our proposal to set up a
mechanism-based PKPD modelling platform. This platform focuses on the
transf
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
> ha
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 GN
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?
http://www.nabble.com/file/
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 hoo
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
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
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 the
> > 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 you
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 ~ SSm
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
ht
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 wo
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,
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
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.
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
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
On 3/28/2008 5:33 AM, Lila86 wrote:
> 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,
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
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
try
master_list <- apply(x, 1, as.list)
HTH
> Dear group,
> I have a dataframe (x). (4 rows, 3 columns for
> example)
>
> I want to create a master list (my_list) from data
> frame.
>
> I want each row of x to be a list.
>
>
> I tried 3 different ways: writing a loop, using apply,
> and go
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 da
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 or
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
summary(
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 <-
data.frame(V1=round(runif(2
59 matches
Mail list logo