BIOL398-01/S10:Sample Microarray Analysis Vibrio cholerae

From OpenWetWare

Revision as of 00:01, 13 April 2010 by Kam D. Dahlquist (Talk | contribs)
Jump to: navigation, search
BIOL398-01: Bioinformatics Laboratory

Loyola Marymount University

Home       People        Molecular Genetics Explorer       HIV Evolution       HIV Structure       DNA Microarrays       Help  

This page has been written with the analysis of the Vibrio cholerae dataset in mind. However, these steps are similar to what needs to be performed with any microarray dataset, although the details will differ with the particular experimental design.

Before we begin...

  • The data from the Merrell et al. (2002) paper was accessed from this page at the Stanford Microarray Database.
  • The Log(base2) of R/G Normalized Ratio (Median) has been copied from the raw data files downloaded from the Stanford Microarray Database.
    • Patient A
      • Sample 1: 24047.xls (A1)
      • Sample 2: 24048.xls (A2)
      • Sample 3: 24213.xls (A3)
      • Sample 4: 24202.xls (A4)
    • Patient B
      • Sample 5: 24049.xls (B1)
      • Sample 6: 24050.xls (B2)
      • Sample 7: 24203.xls (B3)
      • Sample 8: 24204.xls (B4)
    • Patient C
      • Sample 9: 24053.xls (C1)
      • Sample 10: 24054.xls (C2)
      • Sample 11: 24205.xls (C3)
      • Sample 12: 24206.xls (C4)
    • Stationary Samples (We will not be using these, they are listed here for completeness, but do not appear in your compiled raw data file.)
      • Sample 13: 24059.xls (Stationary-1)
      • Sample 14: 24060.xls (Stationary-2)
      • Sample 15: 24211.xls (Stationary-3)
      • Sample 16: 24212.xls (Stationary-4)
  • Download the Merrell_Compiled_Raw_Data_Vibrio.xls file to your Desktop.
    • Save a copy of the file with a different filename that includes your initials and the date. For example, I would call mine "Merrell_Compiled_Raw_Data_Vibrio_KD_20091020.xls".

Normalize the log ratios for the set of slides in the experiment

To scale and center the data (between chip normalization) perform the following operations:

  • Insert a new Worksheet into your Excel file, and name it "scaled_centered".
  • Go back to the "compiled_raw_data" worksheet, Select All and Copy. Go to your new "scaled_centered" worksheet, click on the upper, left-hand cell (cell A1) and Paste.
  • Insert two rows in between the top row of headers and the first data row.
  • In cell A2, type "Average" and in cell A3, type "StdDev".
  • You will now compute the Average log ratio for each chip (each column of data). In cell B2, type the following equation:

and type "Enter". Excel is computing the average value of the cells specified in the range given inside the parentheses. Instead of typing the cell designations, you can click on the beginning cell, scroll down to the bottom of the worksheet, and shift-click on the ending cell.

  • You will now compute the Standard Deviation of the log ratios on each chip (each column of data). In cell B3, type the following equation:

and type "Enter".

  • Excel will now do some work for you. Copy these two equations (cells B2 and B3) and paste them into the empty cells in the rest of the columns. Excel will automatically change the equation to match the cell designations for those columns.
  • You have now computed the average and standard deviation of the log ratios for each chip. Now we will actually do the scaling and centering based on these values.
  • Insert a new column to the right of each data column and label the top of the column as follows: A1_scaled_centered, A2_scaled_centered, etc.
  • In cell C4, type the following equation:

In this case, we want the data in cell B4 to have the average subtracted from it (cell B2) and be divided by the standard deviation (cell B3). We use the dollar sign symbols surrounding the "B" to tell Excel to always reference that cell in the equation, even though we will paste it for the entire column. Why is this important?

  • Copy and paste this equation into the entire column.
  • Repeat the scaling and centering equation for each of the columns of data. Be sure that your equation is correct for the column you are calculating.

Perform statistical analysis on the ratios

We are going to perform this step on the scaled and centered data you produced in the previous step.

  • Insert a new worksheet and name it "statistics".
  • Go back to the "scaling_centering" worksheet and copy the first column ("ID").
  • Paste the data into the first column of your new "statistics" worksheet.
  • Go back to the "scaling_centering" worksheet and copy Column C ("A1_scaled_centered).
  • Go to your new worksheet and click on the B1 cell. Select "Paste Special" from the Edit menu. A window will open: click on the radio button for "Values" and click OK. This will paste the numerical result into your new worksheet instead of the equation which must make calculations on the fly.
  • Go to a new column on the right of your worksheet. Type the header "Avg_LogFC_A", "Avg_LogFC_B", and "Avg_LogFC_C" into the top cell of the next three columns.
  • Compute the average log fold change for the replicates for each patient by typing the equation:

into cell N2. Copy this equation and paste it into the rest of the column.

  • Create the equation for patients B and C and paste it into their respective columns.
  • Now you will compute the average of the averages. Type the header "Avg_LogFC_all" into the first cell in the next empty column. Create the equation that will compute the average of the three previous averages you calculated and paste it into this entire column.
  • Insert a new column next to each average fold change column that you computed in the previous step. Label the column "Tstat". This will compute a T statistic that tells us whether the scaled and centered average log ratio is significantly different than 0 (no change). Enter the equation:
=AVERAGE(N2:P2)/(STDEV(N2:P2)/SQRT(number of replicates))

(NOTE: in this case the number of replicates is 3. Be careful that you are using the correct number of parentheses.) Copy the equation and paste it into all rows in that column.

  • Label the top cell in the next column "Pvalue". In the cell below the label, enter the equation:
=TDIST(ABS(R2),degrees of freedom,2)

The number of degrees of freedom is the number of replicates minus one, so in our case there are 2 degrees of freedom. Copy the equation and paste it into all rows in that column.

  • Insert a new worksheet and name it "forGenMAPP".
  • Go back to the "foldchange_stats" worksheet and Select All and Copy.
  • Go to your new sheet and click on cell A1 and selct Paste Special, click on the Values radio button, and click OK. We will now format this worksheet for import into GenMAPP.
  • Select Columns B through Q (all the fold changes). Select the menu item Format > Cells. Under the number tab, select 2 decimal places. Click OK.
  • Select Columns R and S. Select the menu item Format > Cells. Under the number tab, select 4 decimal places. Click OK.
  • Select Columns N through S and Cut. Select Column B by left-clicking on the "B" at the top of the column. Then right-click on the Column B header and select "Insert Cut Cells". This will insert the data without writing over your existing columns.
  • Delete Rows 2 and 3 where it says "Average" and "StDev" so that your data rows with gene IDs are immediately below the header row 1.
  • Insert a column to the right of the "ID" column. Type the header "SystemCode" into the top cell of this column. Fill the entire column (each cell) with the letter "N".
  • Select the menu item File > Save As, and choose "Text (Tab-delimited) (*.txt)" from the file type drop-down menu. Excel will make you click through a couple of warnings because it doesn't like you going all independent and choosing a different file type than the native .xls. This is OK. Your new *.txt file is now ready for import into GenMAPP. But before we do that, we want to know a few things about our data as shown in the next section.
    • Upload both the .xls and .txt files that you have just created to your journal page in the class wiki. Make sure that your file name is distinct from your other classmates so that nobody overwrites anyone else's file.

Personal tools