Search Archives

>as.numeric(format(x, f="%j")) which is the right code, works perfectly, too. Thomas P.

Hello Andy, does mtext(c("SPP1", "SPP2", "SPP3"), side=1, at=c(1,2,3)) complete your script? Thomas P.

Dear Ekaterina, thanks for your report, can you please give us a minimal reproducible example? Thomas Petzoldt

mike.campana at freesurf.ch wrote: > Question: is there a function that average in a single step over the 3 > columns? You may look for ?aggregate Thomas P.

Dear Eike, all solvers of deSolve are able to solve nonlinear ODE systems. Regards, Thomas

Hi, removing names from a vector can also be forced by using: unname(x) see ?unname that internally assigns NULL to names. Thomas

Omar Lakkis schrieb: > Using system() is theer a way to make the R interpreter not wait for > the command to finish? system("cmd", wait=FALSE) see ?system in online help. Thomas P.

Hi, this is not a problem of R it is a problem of the pdf viewer. Solution: just use an alternative pdf viewer like gsview or (even simpler) SumatraPDF. Hope it helps Thomas P.

Hello, (1) your example is incomplete, so we are unable to identify your problem. (2) please use always plain text and not HTML for posting to the mailing list. Thomas

NEWS, rev. 47094 now says: o nlminb() copies names from 'start' to the parameter vector used (for consistency with optim()). Dear Prof. Ripley, thank you very much for doing this. Thomas Petzoldt

We have been able to reproduce the reported issue on another Linux system: Fedora 19, and the solution was quite simple: The deSolve package must always to be loaded *before* loading the shared library of the compiled model. Thomas

Jacques VESLOT wrote: > I am sorry to ask such question, but I can't find a solution... > > I have a dataframe 'd2004' and I want to remove two columns: > 'd2004$concentration' and 'd2004$stade". d2004$concentration <- NULL d2004$stade <- NULL...

Hello, the following may work: t <- "sample text" x <- 333 plot(1, 1, main = expression(main[x1]), xlab = substitute(xxx ~ xlab[x2], list(xxx=t)), sub = parse(text = paste("sub[", x, "]")) ) Some more examples are on the ?substitute help page...

Hi, I have no example at hand, but the usual way could be indeed to calculate a time series of an adequate statistic (e.g. spatial statistic) from both, the observed and the simulated data and then to apply standard...

Hi Dabing Chen, thanks for the upate, but please note that ?advection.1D says that: "When solved dynamically, the euler method should be used, unless the first-order upstream method is used." Thomas P.

We have been able to reproduce the reported issue on another Linux system: Fedora 19, and the solution was quite simple: The deSolve package must always to be loaded *before* loading the shared library of the compiled model. Thomas

Hi, your code cannot work, because: - the objective function (cost) does not call the ode solver - names between data and model output do not match The docs (for example Vignette 2, p14) shows you how to write an objective function...

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...

TAPO (Thomas Agersten Poulsen) wrote: > Is there a way to convert the data frames into arrays, or am I > going the wrong way about this? This is possible, but the behaviour depends on the datatype, e.g. numeric or character...

Hello Frank, does the following example, what you want? Thomas P. ## some test data x <- data.frame(matrix(rnorm(72*10, mean=50, sd=20), ncol=10)) me <- colMeans(x) sd <- apply(x, 2, sd) coff <- me + 2*sd # see...