On 11/17/06, Turgut Durduran <[EMAIL PROTECTED]> wrote:
> On 11/17/06, Turgut Durduran <[EMAIL PROTECTED]> wrote:
> >> Hello everyone,
> >>
> >> I am ploting a groupeddata object with formula:
> >>
> >> formula(mydatausegroup)
> >> BF ~ HO | ID/Infar/Day
> >>
> >> Using this command:
> >> plot(na.o
Sorry, a slight glitch got in because I subsetted the original data frame.
The last part, that works should read:
dd.nls.d2 <- nls(Lum ~ dd.deriv2(Blev, beta, gamm, GL), data = dd,
start = list(Blev = B[2], beta = B[3:5], gamm = B[1]))
#But not here
confint(dd.nls.d2)
Thank you for your rapid response.
This is reproducible on my system. Here it is again, with, I hope,
sufficient detail to properly document what does not work and what
does on my system,
But my original question, properly motivated or not, concerns whether
there is a way to use or adapt deriv()
Not quite:
> library(lme4)
> sessionInfo()
R version 2.4.0 (2006-10-03)
i386-pc-mingw32
locale:
LC_COLLATE=French_France.1252;LC_CTYPE=French_France.1252;LC_MONETARY=French_France.1252;LC_NUMERIC=C;LC_TIME=French_France.1252
attached base packages:
[1] "methods" "stats" "graphics" "grDevi
I meant that it demonstrates that the number is not an integer to the
reader, regardless of how it is stored. In the reporting of
measurements,
10
is interpreted as an integer, whereas
10.
is interpreted as a real, rounded to two significant digits.
Cheers
Andrew
On Fri, Nov 17, 2006 at
You must be thinking of some other language:
> is.integer(1)
[1] FALSE
On 11/17/06, Andrew Robinson <[EMAIL PROTECTED]> wrote:
> Nice solution.
>
> I believe that showing the decimal point is correct. It demonstrates
> that the number is not an integer.
>
> Cheers
>
> Andrew
>
> On Fri, Nov 17,
Use
distribute.type = TRUE
or
panel = panel.superpose.2.
as an argument to xyuplot. e.g. Using the builtin anscombe data set:
xyplot(y1 + y2 ~ x1, anscombe, type = c("p", "l"), distribute.type = TRUE)
# same
xyplot(y1 + y2 ~ x1, anscombe, type = c("p", "l"), panel = panel.superpose.2)
On
Nice solution.
I believe that showing the decimal point is correct. It demonstrates
that the number is not an integer.
Cheers
Andrew
On Fri, Nov 17, 2006 at 02:43:10PM -0800, RMan54 wrote:
>
> I think that I answered my own question.
>
> Since formatC is an implementation of the C-style form
I would appreciate help trying
to set different plot types in a lattice plot with multiple responses e.g.
xyplot(y+fval~x)
but have "y" as points and "fval" as a line.
I have tried
xyplot(y+fval~x,type=c("p","l"))
but this results in both plots having both types.
Thanks
Ross Darnell
Univer
Hi,
I am using SAM (from siggenes_1.2.11 package) method to select genes from
a microarray data set. After installing the latest R2.4.0 on my computer, to
my surprise the results are totally different from that calculated using
R2.1.1. Even the example code doesn't work the same way under these
Please provide reproducible code which shows the error.
> dd <- structure(list(Lum = c(0.15, 0.07, 0.1, 0.19, 0.4, 0.73, 1.2,
+ 1.85, 2.91, 3.74, 5.08, 6.43, 8.06, 9.84, 12, 14.2, 16.6, 0.1,
+ 0.1, 0.17, 0.46, 1.08, 2.22, 3.74, 5.79, 8.36, 11.6, 15.4, 19.9,
+ 24.6, 30.4, 36.1, 43, 49.9, 0.06, 0.06
Hi,
I'm fitting a standard nonlinear model to the luminances measured
from the red, green and blue guns of a TV display, using nls.
The call is:
dd.nls <- nls(Lum ~ Blev + beta[Gun] * GL^gamm,
data = dd, start = st)
where st was initally estimated using optim()
st
$Blev
Hi R community,
I'd like to automatically create a number of matrices named for
elements of another list or matrix.
#I'm planning to deposit data into matrices named for plates in a platelist
platelist<- c("G2625462", "G2625464", "G2625466", "G2625468",
"G2625470", "G2625472", "G2625474", "G26254
Thanks, Tim
To get the time difference in the proper units, it's as simple as this:
difftime(myfileinfo$mtime[2], myfileinfo$mtime[1], units="hours")
Cheers,
Warren
On 11/13/06, Tim Calkins <[EMAIL PROTECTED]> wrote:
> hardly the most efficient way to go, but consider using a substring
> functi
Thank you, Martin, for pointing this out (not to mention the missing
sessionInfo() piece). Indeed the resid() call works as expected with
the most recent update of lme4.
Cheers, Joel.
Martin Maechler wrote:
>> "Joel" == Joel Dubin <[EMAIL PROTECTED]>
>> on Fri, 17 Nov 2006 16:19:1
Hi there. Thanks for your time in advance.
My final goal is to calculate 1/2*integral of
(f1(x)^1/2-f2(x)^(1/2))^2dx (Latex codes:
$\frac{1}{2}\int^{{\infty}}_{\infty}
(\sqrt{f_1(x)}-\sqrt{f_2(x)})^2dx $.) where f1(x) and f2(x) are two
marginal densities.
My problem:
I have the following R codes
Yes, you can always tweak one number but it may mess up the others. That's
why my example included two numbers. When you use sprintf this way:
v <- c(9.6996, 99.99)
sprintf("%.2f", v)
you get:
"9.70" "99.99"
Now I have 4 significant digits for the 2nd number! The ".2" in the format
string al
I have found that sprintf gets this right, although the syntax of the
command itself is a little less clear.
> sprintf("%.2f", 9.6996)
[1] "9.70"
I hope that this helps,
Andrew
On Fri, Nov 17, 2006 at 02:14:58PM -0800, RMan54 wrote:
>
> This, for example:
>
> v <- c(9.6996, 99.99)
> formatC
Try this where vf is the output of your format statement:
gsub("[.]\\B", "", vf)
Actually I was surprised at this as my reading of ?regex suggests it
should have been \\b but when that did not work I tried \\B and that
did work.I am using
"R version 2.4.0 Patched (2006-10-24 r39722)"
on Wi
I think that I answered my own question.
Since formatC is an implementation of the C-style formatting, I thought that
a "#" as flag could work (for g and G conversions, trailing zeros are not
removed from the result as they would otherwise be). Although not in the
online help, this worked in
Hi there. Thanks for your time in advance.
My final goal is to calculate 1/2*integral of
(f1(x)^1/2-f2(x)^(1/2))^2dx (Latex codes:
$\frac{1}{2}\int^{{\infty}}_{\infty} (\sqrt{f_1(x)}-\sqrt{f_2(x)})^2dx
$.) where f1(x) and f2(x) are two estimated marginal densities.
My problem:
I have the followi
cleaning up my overflowing R mail box...
> "Tamas" == Tamas K Papp <[EMAIL PROTECTED]>
> on Fri, 27 Oct 2006 16:11:46 -0400 writes:
Tamas> Hi, I wonder if it would make sense to make uniroot
Tamas> detect zeros at the endpoints, eg
Tamas> if f(lower)==0, return lower as
This is interesting. I am on Win2000 running R version 2.4.0 Patched
(2006-11-15 r39915) and plotrix 2.1.5. I get the plot that scaled a little
differently that the 2.3.1 example but the dates look fine. I ran the
original code posted by John cane with no changes (copy and paste). I used
the wi
> "Joel" == Joel Dubin <[EMAIL PROTECTED]>
> on Fri, 17 Nov 2006 16:19:16 -0500 writes:
Joel> Hello, Indeed, implementing the fitted and resid calls
Joel> on the fm1 object from example(lmer), the fitted
Joel> function worked, but the resid did not:
>> resid(fm1)
J
This, for example:
v <- c(9.6996, 99.99)
formatC(v, digits=3, format="g")
shows:
" 9.7" " 100"
This is scientifically incorrect for the first number in the sense that I
like to show all 3 significant digits, including trailing zero's.
Is there a way that the first number would show as " 9.70"
I don't know how to use the fitted
> values function with a given function and given input-variables but yet
> unknown result-values.
Take a look at the predict function, ?predict.
Mark Lyman
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz
Gabriele Stocco units.it> writes:
>
> Dear R help,
> I have two dataframes with the same number of rows and columns. Each
> column is for a patient, each row for a variable. The columns and rows
> are matched in the two dataframes, so that each patient (column) and
> each variable (row) has the
I am on ubuntu linux with R 2.4.0 compiled and I get the same picture
as with yours with R 2.4.0.
__
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-guid
Hello,
Indeed, implementing the fitted and resid calls on the fm1 object from
example(lmer), the fitted function worked, but the resid did not:
> resid(fm1)
Error in resid(fm1) : no slot of name "family" for this object of class "lmer"
I am using R 2.4.0 (on Windows XP), with the update of lme4
Hi,
how to derive an estimate of skewness and kurtosis out of a predicted
distribution by quantile regression?
Example:
library(quantreg)
data(airquality)
airq <- airquality[143,]
f <- rq(Ozone ~ ., data=airquality,tau=seq(0.01,0.99,0.01))
predict(f,newdata=airq)
Any suggestions?
Kind regards,
Jeff Miller adsam.com> writes:
> I thought more about this last night after my email about putting an end to
> "Read the Manual" replies.
...
> A potential solution is to offer two lists.one for newbies and one
> for.umm.not-newbies. Some not-newbies may prefer subscribing to both and
> setting up
Dear R help,
I have two dataframes with the same number of rows and columns. Each
column is for a patient, each row for a variable. The columns and rows
are matched in the two dataframes, so that each patient (column) and
each variable (row) has the same position in the two dataframes. The
values a
On 11/17/06, Turgut Durduran <[EMAIL PROTECTED]> wrote:
>> Hello everyone,
>>
>> I am ploting a groupeddata object with formula:
>>
>> formula(mydatausegroup)
>> BF ~ HO | ID/Infar/Day
>>
>> Using this command:
>> plot(na.omit(mydatausegroup), displayLevel=2,layout=c(10,2),aspect=2)
>>
>>
>> This t
On 11/17/06, Turgut Durduran <[EMAIL PROTECTED]> wrote:
> Hello everyone,
>
> I am ploting a groupeddata object with formula:
>
> formula(mydatausegroup)
> BF ~ HO | ID/Infar/Day
>
> Using this command:
> plot(na.omit(mydatausegroup), displayLevel=2,layout=c(10,2),aspect=2)
>
>
> This trellis plot
Hello everyone,
I am ploting a groupeddata object with formula:
formula(mydatausegroup)
BF ~ HO | ID/Infar/Day
Using this command:
plot(na.omit(mydatausegroup), displayLevel=2,layout=c(10,2),aspect=2)
This trellis plot does almost what I want and produces a 10x2 trellis plot,
each panel is la
Hi:
I think I found an error in plot.lm with the option which=5, of course I
can be wrong , as usually happen, but I had work on it for a while and show it
to some other people that work with R, and so far I don't see what I can be
interpreting wrong. I also worked over the plot.lm's co
That did it. However the file "R for windows GUI front-end" is named
"R.cmd" in R's bin folder. When it is associated with the .Rdata
extension, it changes to "R for windows GUI front-end". Weird, but
thanks.
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf
Tomas Goicoa wrote:
> Dear R users,
>
> I am trying to fit the following ANCOVA model in R2.4.0
>
> Y_ij=mu+alpha_i+beta*(X_ij-X..)+epsilon_ij
>
> Particularly I am interested in obtaining estimates for mu, and the effects
> alpha_i
>
>
> I have this data (from the book Applied Linear Statist
> How can write my data-frame object into a SAS data set using R? I know
> write.foreign does this...But I need to know more with examples...Could
> somebody help on this?
I find the easiest way is to write out a csv file, and then import
that into SAS.
Hadley
___
Dear R users,
I am trying to fit the following ANCOVA model in R2.4.0
Y_ij=mu+alpha_i+beta*(X_ij-X..)+epsilon_ij
Particularly I am interested in obtaining estimates for mu, and the effects
alpha_i
I have this data (from the book Applied Linear Statistical Models by Neter
et al (1996), page
I thought more about this last night after my email about putting an end to
"Read the Manual" replies.
It seems that there are a handful of R Super-Programmers. If I were one, I
would get tired of 100's of questions like "How can I get my Excel file into
R!!!". Concurrently, R is experiencing a
Dear HelpeRs,
I wonder if anyone knows of ways to read EViews file types.
I did not find a function in the package 'foreign' and a search query
submitted to http://search.r-project.org was not successful.
Any hint is very much welcome.
Dietrich Trenkler
--
Dietrich Trenkler c/o Universit
Did you try to right click on the file, select "Properties", and in the
"Opens with" dialog select "change" and then select "R for windows GUI
front-end". Apply the changes, and now double (left) click on the
.Rdata file to see if it opens in R's GUI.
Cheers,
Francisco
Lucke, Joseph F wrote:
Somehow, I've managed to have my .Rdata files become `disassociated'
from the R program. I am running Windows XP Pro. I have re-installed R
2.4 in an attempt to have it re-associate itself with .Rdata files, but
to no avail. .Rdata files are now associated with a file compression
program or alter
Hi all,
I am looking for a way to display venn diagrams of a handful of
vector, so that they are connected to the previous and the next.
Something like:
a-b-c-d-e-f-g
where a and b are connected, b is also connected to c, etc.
Any pointers?
Thanks,
Pfizer Global Research and Development is seeking non-clinical
statisticians to support several areas of the drug development pipeline.
Positions available in Connecticut are:
- a non-clinical statistician to support non-clinical areas of R&D,
including: Discovery, Drug Safety R&D, and the
Ph
--- Jim Lemon <[EMAIL PROTECTED]> wrote:
> John Kane wrote:
> > I am having a problem with a gantt chart since
> > moving to R2.4.0. from 2.3.1
> >
> > I made some adaptations to the code from
> >
>
http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=74
> >
> > and successful produc
On Fri, 2006-11-17 at 03:57 -0800, Tobias wrote:
> I just recently purchased a Macbook Pro, Intel Core Duo 2GHz, 512 MB RAM.
>
> After looking through several articles in this forum I found the link to the
> R Mac Wiki. I am not that tech-savy but it appeared to me that this link,
>
> http://wik
Hello,
I'd like to be able to set options for grid and lattice globally, once
for all subsequent plots, just like what ps.options() does. I walked
through the functions of Grid and didn't find it. The following are
related things that are available:
get.gpar()
trellis.par.set()
any ideas? tha
"Gabor Grothendieck" <[EMAIL PROTECTED]> writes:
> > > Aha, the dreaded "..1" argument. Where do I look for documentation
> > > on this?
> > >
> > > [It is mentioned twice in The R Language Definition, but I'm no wiser]
> >
> > It's just the first of the ... arguments, essentially list(...)[[1]
I just recently purchased a Macbook Pro, Intel Core Duo 2GHz, 512 MB RAM.
After looking through several articles in this forum I found the link to the
R Mac Wiki. I am not that tech-savy but it appeared to me that this link,
http://wiki.urbanek.info/index.cgi?IntelMacR
describes how to install
"Rahul Thathoo" <[EMAIL PROTECTED]> writes:
> Hi.
>
> I am trying to get data from mysql database using a couple of queries.
> I do one query to find out the indexes. Then i need to use these
> indexes in another query, but i keep getting errors.
>
> Here is something:
>
> numb <- dbSendQuery(con2
Following suggestions from Prof. Ripley and several others to use gzfile,
here's rough code that will unzip a tgz into your working directory and
return a list of the files. (It doesn't warn you that it is overwriting
files!)
The magic numbers refer to the current tar header specification; the blo
"Gabor Grothendieck" <[EMAIL PROTECTED]> writes:
> I did not realize that ++ was available. Is there a comprehensive list
> somewhere of which operators are available for definition? I searched
> the R Language manual for ++ but that only came up with a reference to C++ .
It doesn't work infix.
I did not realize that ++ was available. Is there a comprehensive list
somewhere of which operators are available for definition? I searched
the R Language manual for ++ but that only came up with a reference to C++ .
On 17 Nov 2006 11:38:43 +0100, Peter Dalgaard <[EMAIL PROTECTED]> wrote:
> [EM
On 11/17/2006 9:15 AM, Robin Hankin wrote:
> Thanks for this Duncan
>
>>>
>
>
>>> Why is Peter Dalgaard's suggestion necessary? Why can't "+"
>>> take more than two arguments?
>>
>> One reason is that it's a binary operator, so that's all it needs to
>> take. We have the sum function for mult
On 17 Nov 2006 15:59:24 +0100, Peter Dalgaard <[EMAIL PROTECTED]> wrote:
> Robin Hankin <[EMAIL PROTECTED]> writes:
>
> > On 17 Nov 2006, at 14:14, Peter Dalgaard wrote:
> >
> > [snip]
> >
> > >
> > > so a+b+c is really (a+b)+c and I was calculating a+(b+c). That's
> > > actually a little bit harde
Thomas Lumley <[EMAIL PROTECTED]> writes:
> It's not so long ago that we had a proposal for "+" to implement a
> non-associative and non-commutative operation: "+.character" was going to
> be string concatenation.
Actually, that *is* associative, but not commutative...
-p
(Who suspects that
Robin Hankin <[EMAIL PROTECTED]> writes:
> On 17 Nov 2006, at 14:14, Peter Dalgaard wrote:
>
> [snip]
>
> >
> > so a+b+c is really (a+b)+c and I was calculating a+(b+c). That's
> > actually a little bit harder because you don't get help from argument
> > matching:
> >
> > "++" <- function(...) i
On 17 Nov 2006, at 14:14, Peter Dalgaard wrote:
[snip]
>
> so a+b+c is really (a+b)+c and I was calculating a+(b+c). That's
> actually a little bit harder because you don't get help from argument
> matching:
>
> "++" <- function(...) if ((n <- nargs()) == 1) ..1 else {
> l <- list(...)
>
On Fri, 17 Nov 2006, Robin Hankin wrote:
>
> OK, I see. But in algebra the "+" symbol is special: it is reserved
> exclusively for associative and commutative operations [thus a+b+c is
> always well-defined]; perhaps the parser could fall in with this
> convention?
>
It's not so long ago that we
Thanks for this Duncan
>>
>> Why is Peter Dalgaard's suggestion necessary? Why can't "+"
>> take more than two arguments?
>
> One reason is that it's a binary operator, so that's all it needs to
> take. We have the sum function for multiple operands.
>
> I would guess the historical reason is
Duncan Murdoch <[EMAIL PROTECTED]> writes:
> By the way, another complaint is that sum() is supposed to be generic,
> but you can't define a sum.matrix() method so that sum(a,b,c) does the
> same as a+b+c when a is a matrix.
Ah, yes. This also means that my suggestion doesn't work unless "+" is
On Fri, 17 Nov 2006, Duncan Murdoch wrote:
[...]
> By the way, another complaint is that sum() is supposed to be generic,
> but you can't define a sum.matrix() method so that sum(a,b,c) does the
> same as a+b+c when a is a matrix. (This would probably be a bad idea
> because people may be relyin
Hi all
I don't know much about R internals and what goes on behind the scenes but
to me it would seem that
lapply(dim(a), function(i) seq(length=i))
accomplishes the same as
jj <- function(i) seq_len(dim(a)[i])
lapply(seq_len(length(dim(a))), jj)
and might be easier to read. Is there a perform
On Fri, 2006-11-17 at 12:26 +, Gavin Simpson wrote:
> On Fri, 2006-11-17 at 12:18 +0100, Alicia Amadoz wrote:
> > Hello,
> >
> > Thank you for your help. I have tried to perform the analysis I wanted
> > with data of example, I mean not real data because I can't provide it
> > here. So, what I
Thank you for your quick help.
Regards
Jaci
Original-Nachricht
Datum: Fri, 17 Nov 2006 13:15:13 + (GMT)
Von: Prof Brian Ripley <[EMAIL PROTECTED]>
An: [EMAIL PROTECTED]
Betreff: Re: [R] Kolmogorov Smirnov Distribution
> On Fri, 17 Nov 2006, [EMAIL PROTECTED] wrote:
>
> > H
You might use the idea below together with sum.exact in the caTools package.
On 11/17/06, Martin Maechler <[EMAIL PROTECTED]> wrote:
> > "Paul" == Paul Smith <[EMAIL PROTECTED]>
> > on Fri, 17 Nov 2006 12:12:52 + writes:
>
>Paul> On 11/16/06, Prof Brian Ripley <[EMAIL PROTECTED
On 11/17/2006 6:57 AM, Robin Hankin wrote:
> Hello everyone
>
> thanks for the replies. My application (predictably) involves
> arbitrary dimensioned arrays, so I will need to generalize
> the suggestions slightly (except Peter Dalgaard's, which
> works out-of-the-box).
>
> At least I wasn't mis
> "Paul" == Paul Smith <[EMAIL PROTECTED]>
> on Fri, 17 Nov 2006 12:12:52 + writes:
Paul> On 11/16/06, Prof Brian Ripley <[EMAIL PROTECTED]>
Paul> wrote:
>> > For my calculations, I am needing to use more
>> floating-point precision > than the default one of R. Is
On Fri, 2006-11-17 at 14:50 +0200, Roy Spitz wrote:
> Hello
>
>
>
> I have 2 rows in a matrix and I want to calculate the Gower Distance between
> the 2 , how can I do it?
> I searched and found nothing that can help me, and my program doesn't know
> the gdist function and I couldn't find it on
For a CSV file I would use 'read.csv()'. I've compiled some notes on reading
in
large data frames into R at
http://www.biostat.jhsph.edu/~rpeng/docs/R-large-tables.html.
-roger
Sumanta Basak wrote:
> Hi R-Experts,
>
>
>
> I'm having a problem with reading a large data file whi
At 15:59 16/11/2006, Peter Bolton wrote:
>Hello!
>I have some data stored into 2 separate csv file. 1 file (called
>A.csv) (12 results named Group1, Group2, Group3, etc...) odds
>ratios, 2 file (called B.csv) 12 corresponded errors.
>How to import that data into R and make forest plot like I saw
Hello Gavin,
Thank you very much for your help. I'm sorry I forgot to include all
commands that I used but next time I will try to write all of them. I
will try with my real data and see how it goes. I think I finally have
understood how capscale works with this kind of data. Thank you.
Regards,
On Fri, 17 Nov 2006, [EMAIL PROTECTED] wrote:
> Hi,
>
> i need to calculate the p-value given the kolmogorv smirnov test
> statistic (two sided). When I look at the function ks.test I found the
> following function which calls C code:
>
> pkstwo <- function(x, tol = 1e-06) {
>if (is.numeric(
Roy Spitz pravi:
> Hello
>
>
>
> I have 2 rows in a matrix and I want to calculate the Gower Distance between
> the 2 , how can I do it?
> I searched and found nothing that can help me, and my program doesn't know
> the gdist function and I couldn't find it on the R help site.
>
>
>
> Can anyon
Hello
I have 2 rows in a matrix and I want to calculate the Gower Distance between
the 2 , how can I do it?
I searched and found nothing that can help me, and my program doesn't know
the gdist function and I couldn't find it on the R help site.
Can anyone help me plz
Thank u all
On Fri, 2006-11-17 at 12:18 +0100, Alicia Amadoz wrote:
> Hello,
>
> Thank you for your help. I have tried to perform the analysis I wanted
> with data of example, I mean not real data because I can't provide it
> here. So, what I have tried is this,
Hi Alicia,
It would have been more helpful if
On 11/17/06, Paul Smith <[EMAIL PROTECTED]> wrote:
> On 11/16/06, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
> > > For my calculations, I am needing to use more floating-point precision
> > > than the default one of R. Is that possible? And, if yes, how?
> >
> > See package gmp (but that will be
On Fri, 2006-11-17 at 05:21 -0500, Chuck Cleland wrote:
> Murray Pung wrote:
> > Is it possible to add standard error bars to the means on interaction plots?
>
> Not with interaction.plot(), as far as I know. You might consider the
> effects package by John Fox. For example:
Another possibili
On 11/16/06, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
> > For my calculations, I am needing to use more floating-point precision
> > than the default one of R. Is that possible? And, if yes, how?
>
> See package gmp (but that will be slow and cumbersome for all but simple
> calculations).
>
> T
Hello again everyone.
I've further added to Martin and Gabor's suggestion an ellipsis to
pass additional arguments to f(). Cut-n-paste below.
BUT.do.index() comes with a Warning: function arow() of the magic
package is much much much faster; use it if at all possible:
a <- array(0,
Hello everyone
thanks for the replies. My application (predictably) involves
arbitrary dimensioned arrays, so I will need to generalize
the suggestions slightly (except Peter Dalgaard's, which
works out-of-the-box).
At least I wasn't missing anything obvious about do.call().
Why is Peter Dalgaa
Hi,
i need to calculate the p-value given the kolmogorv smirnov test statistic (two
sided). When I look at the function ks.test I found the following function
which calls C code:
pkstwo <- function(x, tol = 1e-06) {
if (is.numeric(x))
x <- as.vector(x)
else stop("Argument x mus
Hi All,
I'm using Microsoft Office 2007 Beta version and I have
one file in .xlsx format with size 168MB. How can I use odbcConnectExcel
in this as this function only supports .xls format? I have huge data
containing 603470 rows in a .csv file. Thanks in advance.
S
I am trying to use nlmeODE. But I don't know what is ObsEq. How you can get it
from your ODEs system? Could someone help me out?
Best,
Ron
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Hello,
Thank you for your help. I have tried to perform the analysis I wanted
with data of example, I mean not real data because I can't provide it
here. So, what I have tried is this,
> matrix
[,1] [,2] [,3]
[1,] 0.00 0.13 0.59
[2,] 0.13 0.00 0.55
[3,] 0.59 0.55 0.00
> dist_mat
12
John Kane wrote:
> I am having a problem with a gantt chart since
> moving to R2.4.0. from 2.3.1
>
> I made some adaptations to the code from
> http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=74
>
> and successful produced a simple gantt chart. However
> when I upgraded to 2.4.0
[EMAIL PROTECTED] writes:
> Hi,
>
> You could stack your list in an array and then use apply :
>
> myArray <- array( unlist(lapply(1:3, f) ) , dim =c(2, 3, 3))
> apply(myArray, c(1,2), sum)
Or fixup "+" to take more than two arguments, e.g.
"++" <- function(x, ...) if (nargs() == 1) x else x +
We are happy to present you our new project!Today our online auction bidding
company SmartBill has finally begun successful cooperation with American
financial oganizations.Our online auction bidding system is secured with
many ameican major banks and this is why we, as the official af
We are happy to present you our new project!Today our online auction bidding
company SmartBill has finally begun successful cooperation with American
financial oganizations.Our online auction bidding system is secured with
many ameican major banks and this is why we, as the official af
Murray Pung wrote:
> Is it possible to add standard error bars to the means on interaction plots?
Not with interaction.plot(), as far as I know. You might consider the
effects package by John Fox. For example:
with(ToothGrowth, interaction.plot(dose, supp, len))
ToothGrowth$dose <- as.factor
Hi,
You could stack your list in an array and then use apply :
myArray <- array( unlist(lapply(1:3, f) ) , dim =c(2, 3, 3))
apply(myArray, c(1,2), sum)
Cheers,
Romain
Quoting Robin Hankin <[EMAIL PROTECTED]>:
> Hi
>
> How do I make do.call() take "+" as a function for a list of more
> than t
Marc,
This is an exercise of the lesson STAT2430 at the UCL. Please, do your
homework alone
Marc Feuerstein a écrit :
> Hey listmembers,
> I am desperately trying to write a data frame to a file. Not in CSV, but as
> they "appear" on the screen (nice, easy to read tables). I've read tha
Who can help me in finding the right fitted-value-function?Based upon a
linear model (2-dimensional, 1 result-variable and 1 input-variable).
I found the perfect regression model (not always simply linear) by the
function lm(). Based on the given input and output tuples. For these exact
input valu
Hi Experts,
How can write my data-frame object into a SAS data set using R? I know
write.foreign does this...But I need to know more with examples...Could
somebody help on this?
Thanx in advance...
[[alternative HTML version deleted]]
___
Marc Feuerstein <[EMAIL PROTECTED]> writes:
> Hey listmembers,
> I am desperately trying to write a data frame to a file. Not in CSV, but as
> they "appear" on the screen (nice, easy to read tables). I've read that the
> sink function is the way to go.
>
> I have tried the following code inside
On Fri, 2006-11-17 at 01:39 -0800, Marc Feuerstein wrote:
> Hey listmembers,
> I am desperately trying to write a data frame to a file. Not in CSV,
> but as they "appear" on the screen (nice, easy to read tables). I've
> read that the sink function is the way to go.
>
> I have tried the following
Hey listmembers,
I am desperately trying to write a data frame to a file. Not in CSV, but as
they "appear" on the screen (nice, easy to read tables). I've read that the
sink function is the way to go.
I have tried the following code inside a function.
sink("ABC.txt")
MyFrameA
MyFrameB
sink()
try the following:
mat.lis <- lapply(1:4, f)
matFun(mat.lis, sum)
where
matFun <- function(lis, FUN){
if(!is.list(lis) || !all(sapply(lis, is.matrix)))
stop("'lis' must be a list containing 2-dimensional arrays")
dims <- sapply(lis, dim)
n <- dims[1, 1]
p <- dims[2, 1]
1 - 100 of 103 matches
Mail list logo