1H NMR Analysis (SIMCA): Difference between revisions

From Powers Wiki
No edit summary
No edit summary
 
(8 intermediate revisions by 3 users not shown)
Line 1: Line 1:
[[Category:Protocols]]
[[Category:Metabolomics]]
=1H NMR Analysis (SIMCA)=


==Excel Processing==  
==Excel Processing==  
Secondary observation ID labels can be added by inserting a blank second row and filling in the desired labels. This does NOT affect the raw data. It merely makes “viewing” easier in the SIMCA output.
Secondary observation ID labels can be added by inserting a blank second row and filling in the desired labels. This does NOT affect the raw data. It merely makes “viewing” easier in the SIMCA output.
 
Note:


Use the standard method to remove instrument noise.
Use the standard method to remove instrument noise. Click here [http://bionmr.unl.edu/wiki/Noise_removal_for_PCA] for details.  


For PLS-DA analysis, the y-variable (e.g., the paralysis score for mouse urine metabolomics) values are placed in a row below the row containing the last NMR integral bucket value. When the spreadsheet is transposed in SIMCA, this last row becomes the last column.   
For PLS-DA analysis, the y-variable values are placed in a row below the row containing the last NMR integral bucket value. When the spreadsheet is transposed in SIMCA, this last row becomes the last column.   


After noise removal, the data should be autoscaled. By default, SIMCA-P will use "UV" autoscale the imported data set.     
After noise removal, the data should be autoscaled. By default, SIMCA-P will use "UV" autoscale the imported data set.     
Line 25: Line 19:
4. Click “Open”
4. Click “Open”


5. Click on the comma option. The data formatting should now look correct
5. Click on the comma option  The data formatting should now look correct.


6. Click “OK”
6. Click “OK”
Line 34: Line 28:


9. Respond “No” to the question regarding the one row that is “empty or contains only text”  
9. Respond “No” to the question regarding the one row that is “empty or contains only text”  
   It is probably the row of labels that was inserted using Excel.
   It is the row of labels that was inserted using Excel. You can delete the blank row to avoid this.
 
10. Use the “Commands” button (bottom left) to “transpose” the data set
The rows should now correspond to a given NMR spectrum.  Each row is an “observation”. The integral values contained in each row are referred to as “variables”.


10. Use the “Commands” button (bottom left) to “transpose” the data set. The rows should now correspond to a given NMR spectrum. Each row is an “observation”. The integral values contained in each row are referred to as “variables”.
11. Click on the button at the top of the column to set the first column as the “primary observation ID’s”  
It may already be labeled as primary.


11. After transposing, make sure that first column has been identified as the “primary observation ID’s”.  To do this, click on the button at the top of the column
12. Click on the “Observation IDs primary” button  
    (Note: it may already be labeled as primary) and then click on the “Observation IDs primary” button. The observation IDs should now be color coded with the observation ID primary color (dark yellow or mustard color). If secondary observation IDs were inserted using Excel, then click on the second column button and then click on “Observation IDs secondary”. Again, the column should be color coded to the correct color (light yellow).     
The observation IDs should now be color coded with the observation ID primary color. If secondary observation IDs were inserted using Excel, then click on the second column button and then click on “Observation IDs secondary”. Again, the column should be color coded to the correct color (light yellow).     


12. Click on the butoon for the first row (note: it may already be labeled as primary) and then click on the “Variable IDs primary” button to color code the first row (green).  Next, repeat for the second row containing the ppm ranges that define each bucket. These will be the “Variable IDs secondary” and are turquoise.  Click “Next” button in the lower-right.  
12. Click on the butoon for the first row
It may already be labeled as primary and then click on the “Variable IDs primary” button to color code the first row (green).  Next, repeat for the second row containing the ppm ranges that define each bucket. These will be the “Variable IDs secondary” and are turquoise.  Click “Next” button in the lower-right.  


13. Click the “Finish” button
13. Click the “Finish” button
Line 53: Line 52:
17. Click “Next” then click “Finish”
17. Click “Next” then click “Finish”


18. Using the menu bar, click on the “autofit” button. This should calculate the first and second primary components.  Additional components can be calculated using the “Calculate next component” button.   
18. Using the menu bar, click on the “autofit” button.  
    This should calculate the first and second primary components.  Additional components can be calculated using the “Calculate next component” button.   


19. View the results, click on the “Create four overview plots” button. This produces the “Score Scatter” plot in the upper-left hand corner. The lower-left hand corner contains the “Loading Scatter” plot.   
19. View the results, click on the “Create four overview plots” button.  
This produces the “Score Scatter” plot in the upper-left hand corner. The lower-left hand corner contains the “Loading Scatter” plot.   


20. Expand the score scatter plot for better viewing.   
20. Expand the score scatter plot for better viewing  


21. Click on data point
21. Click on data point
Line 64: Line 65:


23. Choose the color tab and choose coloring type by “identifiers”  
23. Choose the color tab and choose coloring type by “identifiers”  
    The default then is to color by secondary observation IDs.  This uses the labels inserted using Excel. If desired, change the default colors.  
The default then is to color by secondary observation IDs.  This uses the labels inserted using Excel. If desired, change the default colors.  


24. Click “Apply”.
24. Click “Apply”


25. Click “OK”.
25. Click “OK”


==PLS Analysis==
==PLS Analysis==


Note:
1. The data are prepared in an Excel file as described above.   
1. The data are prepared in an Excel file as described above.   


2. The EXCEL file is opened in SIMCA as described above for PCA analysis.
2. The Excel file is opened in SIMCA as described above for PCA analysis.


3. The opened file should then be transposed. Again, the commands button found in the lower left provides access to this command.   
3. The opened file should then be transposed. Again, the commands button found in the lower left provides access to this command.   


4. The label information for observations and variables should be processed as described for PCA analysis above.   
4. The label information for observations and variables should be processed as described for PCA analysis above.   


5. The difference between the PCA approach and the PLS approach occurs with labeling of the variables (i.e., the bucketed intensities and the  
5. The difference between the PCA approach and the PLS approach occurs with labeling of the variables (i.e., the bucketed intensities and the discriminator values).   
discriminator values).   


6. The region containing the bucketed intensities is highlighted. These values are labeled as “x-variables” by clicking on the VARIABLE button  
6. The region containing the bucketed intensities is highlighted.  
in the left hand frame and choosing x-variables.   
These values are labeled as “x-variables” by clicking on the VARIABLE button in the left hand frame and choosing x-variables.   


7. Next, the discriminator values (e.g., the single column of paralysis scores in mouse urine metabolomics) are highlighted.  These values are labeled as “y-variables” by clicking on the VARIABLES button and choosing “y-variables”. The column of y-variables will not have a variable ID number in Row 1 because that row was not created in the ACD processing. Incrementing the variable number and typing the value into the cell will save SIMCA from asking about that issue.  Also, SIMCA does not like a mix of text and numerical values in the “y-variables” column.  I discovered this when I used simply “0, 1, 2, 3, 4” instead of “EAE-0, EAE-1,…,EAE-4”.  One response is to write in the text value using the missing value box in the left hand panel (it is below the “exclude data button”).  This approach preserves the numerical content for PLS.  Another answer is to use all text or all numerical values, but not a mix.
7. Next, the discriminator values should be set as 1 or 0.  


8. Selected data may still be excluded prior to the statistical analysis.
8. Selected data may still be excluded prior to the statistical analysis


9. The dataset is “fit”, additional components can be added/substracted and the results are visualized using the same commands as for PCA analysis described above.
9. The dataset is “fit”, additional components can be added/substracted and the results are visualized using the same commands as for PCA analysis described above.
Line 95: Line 94:
==OPLS-DA Analysis==
==OPLS-DA Analysis==


1. This approach applies orthogonal signal correction prior to PLS analysis. I need to do more reading…  But, I believe that I have basic command sequence needed to explore the approach using SIMCA-P+ Version 12.0.
1. This approach applies orthogonal signal correction prior to PLS analysis. This manual is for SIMCA-P+ 12.0 or higher version.  
 
2. Start from an Excel file that has NOT been noise corrected and NOT been autoscaled.


3. The data file should be first analyzed using either PCA or PLS as described above.
2. Start from an Excel file that has NOT been autoscaled and any outliers should be removed prior to OPLS analysis.


4. The analysis is reportedly hindered by extreme outliers found in PCA 2D scores plot (where did I read this…?...the Umetrics manual maybe…). So, any outliers should be removed prior to OPLS analysis.
3. The data file should be first analyzed using either PCA or PLS as described above to find two unsupervised, separated groups for comparison. Wild-type or control groups should be assigned as "0", mutant or treatment group should be assigned as "1". This column of "0" and "1" should be taken as Y values.


5. Open the Excel file in SIMCA, go to the Dataset pull-down menu.  Choose Spectral Filters.  In the available column, scroll down and select OCS.  Press the button labeled as => to move OCS into the selected column.  Click OK.
4. Follow the same procedure as PCA to import the data.  


6. The OSC panel should appear.  Refer to your NMR spectral data to identify bins that contain strong peaks for metabolites (e.g., citrate peaks in mouse urine spectra).  Highlight several of these bins (maybe 5-10). Click on Y to change the state of these to Y values.  Click on the “Next >” button.  There may be a message regarding exclusion of variables with no variance.
5. Click on Y to change the state of these to Y values.  Click on the “Next >” button.  There may be a message regarding exclusion of variables with no variance. A new Workset window will appear.


7. The result, both in terms of the scores plot and the plots showing the difference between PC score points, depends which peaks are assigned as the Y values.
6. Click "Yes to All" and then another workset window will pop-up with "model number", "type of model", etc.  


8. A new OSC panel will appear.  There should be a table with columns labeled No, Angle in Degreees, Remaining SS in % and Eigenvalue.  Click on the “next component” button.  Generally, two components are recommended by UmetricsClick on the Next button.
7. Double click it and click "Workset" icon on that pop-up window.   


9. Check the destination folder and file name. Click on the “Finish” button.
8. Select all the variables, and use "Par" for scaling. Model type "OPLS/O2PLS".


10. Read and close the OSC message box.
9. Click Autofit button on the toolbar


11. The current model will probably say “PLS <unfitted>”  in the Type column.  Fit the data using the usual commands for autofit and plot visualization


12. Go the Analysis pull-down menu and select Change Model Type.  From the list, choose OPLS/O2PLS.  Choose the Analysis pull-down menu and select Autofit (…or just use the Autofit button on the toolbar).
[[category:Data_Processing_and_Analysis|H NMR Analysis (SIMCA)]]

Latest revision as of 05:14, 20 January 2022

Excel Processing

Secondary observation ID labels can be added by inserting a blank second row and filling in the desired labels. This does NOT affect the raw data. It merely makes “viewing” easier in the SIMCA output.

Use the standard method to remove instrument noise. Click here [1] for details.

For PLS-DA analysis, the y-variable values are placed in a row below the row containing the last NMR integral bucket value. When the spreadsheet is transposed in SIMCA, this last row becomes the last column.

After noise removal, the data should be autoscaled. By default, SIMCA-P will use "UV" autoscale the imported data set.

PCA Analysis

1. Start a “new” project by opening the desired integral table spreadsheet

2. Click on new project icon

3. Select the desired Excel file

4. Click “Open”

5. Click on the comma option The data formatting should now look correct.

6. Click “OK”

7. Project type should be SIMCA-P Project

8. Click “next” button

9. Respond “No” to the question regarding the one row that is “empty or contains only text”

  It is the row of labels that was inserted using Excel. You can delete the blank row to avoid this. 

10. Use the “Commands” button (bottom left) to “transpose” the data set The rows should now correspond to a given NMR spectrum. Each row is an “observation”. The integral values contained in each row are referred to as “variables”.

11. Click on the button at the top of the column to set the first column as the “primary observation ID’s” It may already be labeled as primary.

12. Click on the “Observation IDs primary” button The observation IDs should now be color coded with the observation ID primary color. If secondary observation IDs were inserted using Excel, then click on the second column button and then click on “Observation IDs secondary”. Again, the column should be color coded to the correct color (light yellow).

12. Click on the butoon for the first row It may already be labeled as primary and then click on the “Variable IDs primary” button to color code the first row (green). Next, repeat for the second row containing the ppm ranges that define each bucket. These will be the “Variable IDs secondary” and are turquoise. Click “Next” button in the lower-right.

13. Click the “Finish” button

14. Exclude the solvent region and “ends” of the spectra

15. Highlight the desired rows by dragging the cursor along the top set of buttons and then click the “exclude” button (along the left edge)

16. Repeat for each desired region

17. Click “Next” then click “Finish”

18. Using the menu bar, click on the “autofit” button.

   This should calculate the first and second primary components.  Additional components can be calculated using the “Calculate next component” button.  

19. View the results, click on the “Create four overview plots” button. This produces the “Score Scatter” plot in the upper-left hand corner. The lower-left hand corner contains the “Loading Scatter” plot.

20. Expand the score scatter plot for better viewing

21. Click on data point

22. Right-click and choose “Properties”

23. Choose the color tab and choose coloring type by “identifiers” The default then is to color by secondary observation IDs. This uses the labels inserted using Excel. If desired, change the default colors.

24. Click “Apply”

25. Click “OK”

PLS Analysis

1. The data are prepared in an Excel file as described above.

2. The Excel file is opened in SIMCA as described above for PCA analysis.

3. The opened file should then be transposed. Again, the commands button found in the lower left provides access to this command.

4. The label information for observations and variables should be processed as described for PCA analysis above.

5. The difference between the PCA approach and the PLS approach occurs with labeling of the variables (i.e., the bucketed intensities and the discriminator values).

6. The region containing the bucketed intensities is highlighted. These values are labeled as “x-variables” by clicking on the VARIABLE button in the left hand frame and choosing x-variables.

7. Next, the discriminator values should be set as 1 or 0.

8. Selected data may still be excluded prior to the statistical analysis

9. The dataset is “fit”, additional components can be added/substracted and the results are visualized using the same commands as for PCA analysis described above.

OPLS-DA Analysis

1. This approach applies orthogonal signal correction prior to PLS analysis. This manual is for SIMCA-P+ 12.0 or higher version.

2. Start from an Excel file that has NOT been autoscaled and any outliers should be removed prior to OPLS analysis.

3. The data file should be first analyzed using either PCA or PLS as described above to find two unsupervised, separated groups for comparison. Wild-type or control groups should be assigned as "0", mutant or treatment group should be assigned as "1". This column of "0" and "1" should be taken as Y values.

4. Follow the same procedure as PCA to import the data.

5. Click on Y to change the state of these to Y values. Click on the “Next >” button. There may be a message regarding exclusion of variables with no variance. A new Workset window will appear.

6. Click "Yes to All" and then another workset window will pop-up with "model number", "type of model", etc.

7. Double click it and click "Workset" icon on that pop-up window.

8. Select all the variables, and use "Par" for scaling. Model type "OPLS/O2PLS".

9. Click Autofit button on the toolbar