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

  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

  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",

Greetings, Trying to find a solution to creating home directories by using winbind. Have this setup: /etc/pam.d/system-auth: session required /lib/security/pam_mkhomedir.so skel=/etc/skel/ umask=0077 With smb.conf: # Global parameters [global] workgroup = domain server string = Monarch Server security = DOMAIN map to guest = Bad User obey pam

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.

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

hi I am a student in FME in Greece and with my friends, we try to make a project in computer lesson. We try to make a grid with some Pentium III and we use Knoppix Cluster. But we have a surius problem.We can not make the other pc to boot over the network because the bios of the pc do not support the network boot order.So we thought to make some boot loaders but the think is that we dont know

Dear R-listers, what is the best way to draw an XY scatterplot with symbol size proportional to the value taken on by a third variable? Thank you very much for any help, Maurizio Pisati