Displaying 1 result from an estimated 1 matches for "b559lp".
Did you mean:
b559hp
2008 Nov 06
1
nls: Fitting two models at once?
...cytb.v
# fit first spectrum to two reference spectra
asfe.fit <-
nls(
asfe ~ (cytf.v * 28) * CYTF + hp.v * B559HP,
start = list(cytf.v = 0.00005, hp.v = 0.000003) # arbitrary
);
# fit second spectrum to three reference spectra
dias.fit <-
nls(
dias ~ (cytb.v * 24.6 * 2) * CYTB6 + lp.v * B559LP + c550.v * C550,
start = list(cytb.v = 1, lp.v = 1, c550.v = 1) # arbitrary
);
# draw stuff
plot(
1, 2,
type="n",
xlim = c(540, 575),
ylim=c(-0.002, 0.008),
);
# first spectrum and fit
lines(wl, asfe, type="b", pch=19); # solid circles
lines(wl, fitted(asfe.fit), col = &...