Skip to content
Prev 176827 / 398503 Next

Curve fitting,FDA for biological data

Dear Dr Gkikopoulos: 


      1.  Have you looked at "bioconductor.org"?  They have substantive 
extensions to R specifically for "genomic data", which I assume would 
include chromosome. 


      2.  To "identify periodicities at different timescales", I agree 
with Stephen that "spectrum" would likely help. 


      3.  The best software to "fit data into discrete number of curves" 
depends on the particular "discrete number of curves" you want to 
consider and how you want to "fit data into" them.  A finite Fourier 
series could be the best tool IF the the multiple periodicities are all 
integer fractions of a common scale.  In that case, using a "fourier" 
base in the "fda" package could be your method of choice.  Otherwise, 
you might consider Bayesian Model Averaging.  RSiteSearch("Bayesian 
Model Averaging") produced 80 hits for me just now, and 
RSiteSearch("Bayesian Model Averaging", "function") produced 60.  
"RSiteSearch.function" in the "RSiteSearch" package [available via 
install.packages("RSiteSearch",repos="http://r-forge.r-project.org")] 
told me that 27 of the 60 were in the "ensembleBMA" package, and another 
14 were in the "BMA" package. 


      4.  The best way to "compare data from different experiments" 
depends on your evaluation of "3" above.  The "fda" package includes an 
"fRegress" function that might be useful. 


      Hope this helps. 
      Spencer Graves
trias wrote: