[R] Help with modFit of FME package

[David Winsemius]

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