PMF coding

From Jimenez Group Wiki
Revision as of 07:26, 6 May 2008 by DonnaS (talk | contribs) (PMF Analysis)
Jump to: navigation, search

PMF Analysis

To Do Priority

  • Allow user to find solutions ranging over many seed values instead of fpeak values

To Do (Less Urgent)

  • Allow users to run PMF on partitioned drives

PMF Evaluation Panel

  • When more than 10 traces are on a plot, make sure the colors get looped
  • Colors in the diurnal plots should be applied to the median waves (when the factor is set to -1)
  • Typo in the popped current t series graph - "Curret" instead of "Current"
  • On the diurnal plots, have checkboxes for "Display Mean" and Display Median" so that you could display one or the other or both
  • Double check the calculation of the residual errors - from Ingrid on 3/31/08 "Messing around with the panel, I'm wondering if there's a calculation error somewhere in the residuals. For the dataset I'm working on now, I get negative values for the Sum(abs(Resid)) / Sum(Signal)"
  • List from Ingrid on 3/20/08:

1. Are mass fraction bars normalized to total measured mass instead of for total fit mass? We should add a wave for "Residual" (it could probably be black).

2. Can we (doesn't have to be this week) make a tseries plot of % of each tseries (so stacked factors as %) including the Residual (so it addes up to 100%) (like Figure 1d in paper http://www.atmos-chem-phys.net/5/3289/2005/acp-5-3289-2005.pdf )

3. Is is possible to make the axes on the RR plot have the lower value be the min(0, wave minimum)?

4. (Probably also lower priority) It would be nice if the "factors to plot" box was enabled. Most useful plots to enable this way are residuals and bars (though bars isn't very reliable).

5. I think it would be helpful to people to have a button somewhere for "save this solution". I do this with a function in my "notUsed" ipf called pmf_makeStaticFactorWvs (and you could remove the part where it asks for a datafolder list; I built that for some synthetic data runs). It has a naming convention and saves the waves to root:StaticFactors: in a subfolder for the number of factors in the solution (e.g., w3factors).

6. It would be nice if the status of the XaxisTo100 box for the factors plot was examined when remaking plots after moving the sliders.

7. People have asked me before if more of the axes can have labels, like Q/Qexp for the Q plots, Mass on the total Tseries plot, and some sort of # for the Current Species Histogram.

8. Another longer-term thing might be to make a table of string waves where the user can name their factors, and then make bar plots with the same factors in the same color/order. Ken has done a version of that in the attached plot.

PMF Scatter Panel

  • List From Ingrid on 3/20/08

1. The user can't modify the R vs Fpeak plot in the panel (change wave colors, axes, etc.).

2. When you change the slider in the scatter panel, the R vs Fpeak plot doesn't update.

3. In the R vs Fpeak plot, I think the row index (ExtSpecies) for RcorrMx4d_foo[ExtSpecies][FactorNum][*][pdex] is off; it looks like you're starting from 1 instead of from 0.

4. Scatter plots should be the same color as their factor label. Factor 1 profile or time series should be black (like its factor label).

5. (possibly also lower priority) It would be nice to be able to compare more things (that I probably didn't tell you about before)! I'd divide them into two categories:

  a)  All of the factor mass spectra (or time series) in one solution vs each other (a la a Scatter Plot Matrix, which can be found under Windows/New/Scatter Plot Matrix).  This can be used for assessing rotation and I think Manjula does it.
  b) Comparing the various types of Residual (time series or mass spectrum) the the factors (this how we see, e.g., that the residual in Pittsburgh looked like the OOA-II time series).