Michael,
I'm also playing with the nls function trying to get it to work
with a Gaussian. My lines that I have at the moment and I hope will help you
are
class(fo <- (x ~ (A/(sig*sqrt(2*pi)))* exp(-1*((bin-mu)^2/(2* sig^2)))))
nls.AB <- nls(fo,data=freq.tab, start= list(A=0.1*len, mu=0.01, sig=0.5),
trace=TRUE)
So first I create the eq to put into the nls function, my var
are the freq.tab which is a list of some binned frequency into 0.1 units and
the len is amount of features in the data. I'll say quickly that at the
moment I am getting a 'singular gradient error'. I hope that this helps.
Cheers.
Paul Benton
Research Technician
Mass Spectrometry
o The
/
o Scripps
\
o Research
/
o Institute
[[alternative HTML version deleted]]