# Endy:F2620/Data Processing/Algorithm

### From OpenWetWare

(Difference between revisions)

Line 13: | Line 13: | ||

* evaluate the fitting equations for the total synthesis at many timepoints and select the hightest value from cog-AHL and store this time value. | * evaluate the fitting equations for the total synthesis at many timepoints and select the hightest value from cog-AHL and store this time value. | ||

* output to transfer.tsv all values of total function for all series at time calculated above and low/high repeat, st dev (Mathworld definition, I think it is what you call standard error), 95% confidence | * output to transfer.tsv all values of total function for all series at time calculated above and low/high repeat, st dev (Mathworld definition, I think it is what you call standard error), 95% confidence | ||

- | * plot using gnuplot 3D superimposed lines "ser*.pdf" files (we don't use them anymore) surfaces (those are the green ones we post on wiki) into files "sersurf*.pdf") and the transfer functions to "tranfer.pdf" with 95% confidence intervals, transferlh.pdf transfer functions with low/high errors, | + | * plot using gnuplot 3D superimposed lines "ser*.pdf" files (we don't use them anymore) surfaces (those are the green ones we post on wiki) into files "sersurf*.pdf") and the transfer functions to "tranfer.pdf" with 95% confidence intervals, transferlh.pdf transfer functions with low/high errors, where * is series number |

- | where * is series number | + | |

Notes: | Notes: | ||

1. Evaluation of the total synthesis function is done using arbitrary precision numbers. B/c they exceeded double range | 1. Evaluation of the total synthesis function is done using arbitrary precision numbers. B/c they exceeded double range |

## Revision as of 11:50, 5 September 2006

- Load the data from the excel file. Read headers of the column to know the colony number and AHL type used. Form a lookup table. Separate the medium column at this point. Count how many repeats there is for each type.
- Find the GFP background by finding the non-induced column
- GFPpre-fit by fitting all the GFP background columns
- Fit medium for OD (find mean)
- For each OD subtract this mean (media)
- For each GFP we subtract the fitted background calculated for non-induced cells of this type (e.g. for cog-AHL we subtract fitted results for non induced cog-AHL from the raw data)
- Calibrate GFP relative units to conc of GFP (from Barry's calibration run): 1.16*10e-6*GFP+9.95*10e-4
- Fit GFP
- Fit OD (don't fit media anymore)
- Plot GFP, GFP/OD, dOD/dt, gamma, (GFP/OD)/dT, total sythesis and output (this is for debugging and outputs a huge .pdf and .tex)
- to params.csv outputs all the fitting parameters and errors (the errors are not really important)
- to surfaces.tsv output points for superimposed transfer functions in 3D (so you can plot the 3D surfaces yourself)
- evaluate the fitting equations for the total synthesis at many timepoints and select the hightest value from cog-AHL and store this time value.
- output to transfer.tsv all values of total function for all series at time calculated above and low/high repeat, st dev (Mathworld definition, I think it is what you call standard error), 95% confidence
- plot using gnuplot 3D superimposed lines "ser*.pdf" files (we don't use them anymore) surfaces (those are the green ones we post on wiki) into files "sersurf*.pdf") and the transfer functions to "tranfer.pdf" with 95% confidence intervals, transferlh.pdf transfer functions with low/high errors, where * is series number

Notes: 1. Evaluation of the total synthesis function is done using arbitrary precision numbers. B/c they exceeded double range