search for: b559hp

Displaying 1 result from an estimated 1 matches for "b559hp".

2008 Nov 06
1
nls: Fitting two models at once?
...ref <- read.table("http://home.arcor.de/ballaschk/cytbf-help/reference.csv", sep="\t", dec=".", header=T); attach(ref); ### FITTING, problem: 2*cytf.v == 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, typ...