On my trials, after eliminating all the extra matrix<->dgeMatrix conversions,
using expm() and the method below were equally fast.
Michael
On 19 Aug 2011, at 1:32AM, Ravi Varadhan wrote:
> Which is why I said it applies when the system is "diagonalizable". It won't
> work for non-diagonaliza
Which is why I said it applies when the system is "diagonalizable". It won't
work for non-diagonalizable matrix A, because T (eigenvector matrix) is
singular.
Ravi.
From: peter dalgaard [pda...@gmail.com]
Sent: Thursday, August 18, 2011 6:37 PM
To: Ravi
On Aug 17, 2011, at 23:24 , Ravi Varadhan wrote:
> A principled way to solve this system of ODEs is to use the idea of
> "fundamental matrix", which is the same idea as that of diagonalization of a
> matrix (see Boyce and DiPrima or any ODE text).
>
> Here is the code for that:
>
>
> nlogL2
Sorry, embarrassing stupid error: of course '.libPaths(path.r.lib)' will
do the trick for me.
Regards,
Janko
Original Message
Subject:Bug: argument 'lib.loc' not passed on appropriately in
'library()'?
Date: Thu, 18 Aug 2011 19:31:02 +0200
From: Janko Thyson
To:
Dear list,
I'm experimenting with setting up custom 'lib' and 'destdir' directories
for my R packages. It seems to me that 'library()' does handle a custom
'lib.loc' argument the way it should for an arbitrary package but NOT
for its dependencies on other packages. The latter are looked for in
2011/8/18 Krystian Radlak
>
> Mayby it's possible to do the same thing easier...
>
> You shouldn't be disappointed ;-)
> a=ff(1, length=15*10^6)
> k=bit(15*10^6)
> a=ff(2,15*10^6)
> system.time(k<-a[,chunk(1,length(a),1e4)]<1)
user system elapsed
0.40 0.23 0.72
--
Mi³ego dnia
[[altern
W dniu 18 sierpnia 2011 16:26 u¿ytkownik £ukasz Rêc³awowicz <
lukasz.reclawow...@gmail.com> napisa³:
>
> It should be ok.
>
It wasn't, sorry! This seems to do right job.
library(ff)
a=ff(1, length=15*10^2)
k=bit(15*10^2)
a=ff(-1,15*10^2)
system.time(k<-a[,i]<1)
--
Mi³ego dnia
[[alter
Recently on stackoverflow following problem came up :
http://stackoverflow.com/questions/7103429/all-the-connections-are-in-use-execution-halted/7108799#7108799
A reproducible example: When trying
replicate(200,tryCatch(read.table("this.is.no.file"),
warning=function(w){print(showConne
2011/8/18 Krystian Radlak
> Yes, it's necessery
>
Really? How about something like this:
system.time(k[i]<-a[i]<1)
It should be ok.
--
Mi³ego dnia
[[alternative HTML version deleted]]
__
R-devel@r-project.org mailing list
https://stat.eth
Yes, it's necessery. Operator < is not definied in ff packages.
str(a<1)
logi(0)
Mayby it's possible to do the same thing easier...
--
View this message in context:
http://r.789695.n4.nabble.com/Use-logical-expresion-on-ff-object-tp3752127p3752673.html
Sent from the R devel mailing list archiv
2011/8/18 Krystian Radlak
> Any suggestions?
>
Do you really, really need that loop and chunk...?!
--
Mi³ego dnia
[[alternative HTML version deleted]]
__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
I want to do operation with ff vector from packages ff like simple operation
in R like this. (I have read data from file)
Example
a=1:10
b=a[a>5]
how to do this operation with ff packages:
Example :
a=ff(1, length=15*10^7)
I know that I could get logical vector like this:
k=bit(15*10^7) # from
12 matches
Mail list logo