similar to: Help with modFit of FME package

* Apologies for multiple posting * I attached to my previous e-mail a .r file, and it was not permitted by the rules of the mailing lis. Again, please receive my sincere apologies for this. I re-send again the e-mail with .txt attachemnt in the hope someone an help me to solve my problem. I'm trying to fit a set an ODE to an experimental time series. In the attachment you find the R code I

Dear all, I used modFit of the package FME to fit a set of ODE to a ste of eperiemntal data. The summary of this fit give me the following error > summary(Fit) Residual standard error: 984.1 on 452 degrees of freedom Error in cov2cor(x$cov.unscaled) : 'V' is not a square numeric matrix In addition: Warning message: In summary.modFit(Fit) : Cannot estimate covariance; system is

Hello, I would like to use modFit and modCost from the package FME to find the optimal initial time t0 of a process. For simplicity, the process is either "off" (value 0) or "on" (value h). So I have a data vector with some zeros followed by some h's, e.g. > c(0,0,0,2,2,2,2,2,2,2) [1] 0 0 0 2 2 2 2 2 2 2 (hence h=2 here). I want to find the best guess for the initial

Dearl all, I have to fit a function y = f(x1, x2) to data experiemntal data describing the measured behavior of y. x1 and x2 are the independent variables. Could you suggest me wich R package can I use for this purpose? Thanks, Paola. -- *Paola Lecca, PhD* *The Microsoft Research - University of Trento* *Centre for Computational and Systems Biology* *Piazza Manci 17 38123 Povo/Trento,

I am attempting to calibrate a forest growth model using FME. I have several spatially separate plots with time series growth data for each(eg, height, basal area, etc.). Thus 2 independent variables: time and plot. Is there a way that I can calculate a model cost that takes into account these different plots? I notice in the documentation for the modCost function that there is a note

Dear all, I tried to use nls, but I got the following error Error in numericDeriv(form[[3L]], names(ind), env) : Missing value or an infinity produced when evaluating the model Any suggestion? Thanks, Paola. The code I wrote is Data_pp2_mrna <- data.frame( p1 = protein_1, p6 = protein_6, pp2_mrna

Dear All,   please help with some thoughts on overcoming the following issues, if possible:   #R Code   require(deSolve) require(FME) pars <- list(k = 0.06,v=18) intimes <- c(0,0.5,12,12.5,50) input   <- c(800,0,800,0,0) forc <- approxfun(intimes, input, method="constant") model <- function(pars, times=seq(0, 50, by = 1)) {   derivs <- function(t, state, pars) {     

Questions about the port of Xen to ARM with virt extensions are best posted to xen-devel, the xen-arm list focuses on the PV port. Adding xen-devel since it seems you are mainly asking about the w/-virt-extensions port. On Fri, 2012-05-18 at 12:38 +0100, Krishna Pavan wrote: > Please inform the status of Xen-ARM for Cortex-A8 CPU''s on FME from > ARM. AFAIK no one has tried the

Hi, all I can not figure this out, please have a look and help me out. thank you! Note: this is in SPLUS, not R. I have following code *********************************** modfit<-function(yir,yew, ft) { n<-length(yew) yew<-yew[1:(n-1)] yy<-yir-ft xx<-yew-ft n<-length(xx) xx0<-xx[2:n] yy0 <-yy [2:n] xx1<-xx[1:(n-1)] fit <- garch(yy0~xx0 + xx1+var.in.mean,

This is not so much a bug as an infelicity in the code that can easily be fixed. The initialize expression in the quasi family function is, (uniformly for all links and all variance functions): initialize <- expression({ n <- rep.int(1, nobs) mustart <- y + 0.1 * (y == 0) }) This is inappropriate (and often fails) for variance function "mu(1-mu)".

All- I am interested in estimating a parameter that is the starting value for an ODE model. That is, in the typical combined fitting procedure using nls and lsoda (alternatively rk4), I first defined the ODE model: minmod <- function(t, y, parms) { G <- y[1] X <- y[2] with(as.list(parms),{ I_t <- approx(time, I.input, t)$y dG <- -1*(p1 + X)*G +p1*G_b dX <-

On Fri, 6 Dec 2024 12:29:03 +0100 Mitja Tav?ar via samba <samba at lists.samba.org> wrote: > Il 06/12/24 10:45, Rowland Penny via samba ha scritto: > > On Fri, 6 Dec 2024 10:19:31 +0100 > > Mitja Tav?ar via samba <samba at lists.samba.org> wrote: > > > >> But the machine accoounts are not completely lost, they seem lost > >> if i query the Read

  Hi, All,   This is the first time I seriously use this package. However, I am confused that the result is quite unstable. Maybe I wrote something wrong in the code? So could anybody give me some hint? Many thanks.   My test model is really simple. Y_t = X_t * a_t + noise(V)，(no Intercept here) a_t = a_{t-1} + noise(W)   I first run the following code: (I shall provide data at the end of the

Hi, I am new in R and I have problem with reading this data file 0 TITLE Title 0 XLEGEND Legend -1 LABEL x 1 1 12 1 2 30 1 3 34 I want to read only lines starting with 1 (it indicates 1st plotting line) and create data set from second and third value on this row. Thank for advice Jakub Zlamal -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list

On Tue, 3 Dec 2024 09:15:36 +0100 Mitja Tav?ar via samba <samba at lists.samba.org> wrote: > Hi, i have some problems with a recently joined Read Only Domain > controller. > > I had 2 Domain Controllers based on Windows Server 2019 (hosts > vmw2srvdc1 an vmw2srvdc2). I and i recently added a new site (PSN) > and Read Only DC in this second site based on samba (host

  Hello All,   I am working with dlm for the purpose of estimating and forecasting with a Kalman filter model. I have succesfully set up the model and started generating results. Of course, I need to somehow be sure that the results make sense. Without any apparent target to compare with, my natural selection is the results by odinary least square. The idea being that if I choose a diffuse prior,

Hi everybody, I want to download the spatial library but I don't find it at the CRAN list package. Where can I find the package ? I have to download it to use the function correlogram. Thanks in advance, Paola. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Hello to everybody ! I don't have a question directly connected with R, but it takes starting point from the use of the R function Krig for kriging interpolation and from the idea of making a different implemetation of the function for my purposes. Doas anyone know where can I find free C or C++ codes of kriging ? As always, thanks a lot in advance, Paola.

Il 03/12/24 10:13, Rowland Penny via samba ha scritto: > On Tue, 3 Dec 2024 09:15:36 +0100 > Mitja Tav?ar via samba <samba at lists.samba.org> wrote: > >> Hi, i have some problems with a recently joined Read Only Domain >> controller. >> >> I had 2 Domain Controllers based on Windows Server 2019 (hosts >> vmw2srvdc1 an vmw2srvdc2). I and i recently

Author: jrutt Repository: /hg/zfs-crypto/gate Revision: 41d7d68d4ab013616c499ed0e497bffb54c8cc73 Log message: 6232309 need configurable limit on number of open FMEs Files: update: usr/src/cmd/fm/modules/common/eversholt/eft.c update: usr/src/cmd/fm/modules/common/eversholt/eft.conf update: usr/src/cmd/fm/modules/common/eversholt/fme.c update: usr/src/cmd/fm/modules/common/eversholt/fme.h