Help with modFit of FME package
On Aug 1, 2011, at 3:41 AM, Paola Lecca wrote:
Dear R users, I'm trying to fit a set an ODE to an experimental time series. In the attachment you find the R code I wrote using modFit and modCost of FME package and the file of the time series.
This is getting a bit tiresome. None of the three duplicate such messages have had successful attachment of any code. Why don't you look at what got distributed to the list? The rule I have developed is that I should assume that any file not ending in .pdf or .txt will get scrubbed by the mail server. I realize it is not an exact rule, but it keeps me from submitting files ending in .r or .rdata because I know they will get scrubbed. For some reason a recent rewrite of the Posting Guide appears to have left out this information which my memory tells me used to be there last year.
David. > > When I run summary(Fit) I obtain this error message, and the values > of the > parameters are equal to the initial guesses I gave to them. > > The problem is not due to the fact that I have only one equation (I > tried > also with more equations, but I still obtain this error). > > I would appreciate if someone could help me in understanding the > reason of > the error and in fixing it. > > Thanks for your attention, > Paola Lecca. > > Here the error: > >> summary(Fit) > > Parameters: > Estimate Std. Error t value Pr(>|t|) > pro1_strength 1 NA NA NA > > Residual standard error: 2.124 on 10 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 > singular > > > > > ______________________________________________ > R-help at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > > > > > -- > *Paola Lecca, PhD* > *The Microsoft Research - University of Trento* > *Centre for Computational and Systems Biology* > *Piazza Manci 17 38123 Povo/Trento, Italy* > *Phome: +39 0461282843* > *Fax: +39 0461282814* > <wild_pp1_mrna.txt>______________________________________________ > R-help at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. David Winsemius, MD West Hartford, CT