The Thanksgiving weekend provided time to ponder/reflect on the SRM stats that I’ve run thus far, what else needs to be done, and how to finish up. Received guidance from Dave-o. Note: next time, include a housekeeping protein in list of targets.

This is my “finalized” list of SRM & environmental stats to run. In the last few days I’ve completed much of this. In bold are those remaining. Tomorrow I’ll hopefully be done, and will post my scripts and results.

SRM Protein & Environmental Data Analysis Steps

Each Protein: (assume proteins are independent)

1. Test for normality
2. Lambda transformation
3. Test for normality post transformation
4. Assess outliers, remove if necessary
5. N-way ANOVA by: a) location b) habitat c) site d) region
6. Determine P-adjusted, correct for multiple comparisons (bonferroni method, P/13)
7. Post-hoc test to ID differences (is this really necessary?)
8. Ultimate goal: which proteins are different between locations?
9. Compare total abundance between sites (sum peptide abundance)

Each environmental variable:

1. Download tidal chart data for each site
2. Edit pH, DO & Salinity data:
   a. Remove data from exposed time points, as determined from tidal charts
   b. Identify and remove outliers from pH, DO & Salinity data
   c. Recombined outlier-scrubed data with Temp, Tide data.
3. Assess Normality of each env. variable (all time points)
   • Found to be non-parametric (pH is kinda, but let’s assume not). Dataset is large (>6000 for each parameter), so did not determine lambda via tukeytransform function. Instead, used Krusgal-Wallis non-parametric analysis in lieu of ANOVA
4. KW test for each env. variable by location, by region
5. Dunn Test post-hoc test to ID differences
6. Use bonferroni correction for P-adjusted in tests
7. Ultimate goal: which env. variables are different between locations?
   a. basically all of them.

Prep for regression model:

1. Calculate summary statistics: mean, variance, sd, min, max, median, %>1 sd from mean, %>2 sd from mean
2. Plot() all env. variables- are any linearly related, aka not independent? If so, need to include interaction parameter in regression model.
3. Plot() protein peptides against each other to confirm linear correlation; equation should be ~1:1.
4. If all correlated select 1 peptide to use in regression model; highest abundance is best.

Run regression models for each representative peptide:

1. Step-wise linear regression models with all env. variables; I would expect that only the variables that were found to be different via the ANOVA would significantly contribute to the model
2. General linear model with variables ID’d in step-wise lm
3. Figure out when to add a constant, and if I should do that in this scenario
4. Run anova on best fit model, find P-value of the env. variables to determine confidence in the influence of each env. variable on proteins.
5. Run model on the other peptides in the protein (not used as representative peptides); ID the R^2 and P values