7 mar 2012

NIT: Fatty acids study in R - Part 004

We already know that sample 66 has a MD of 11.6. Let´s see the values for the six constituents for this sample:
> fattyac_msc[66,1:6]
   C16_0  C16_1  C18_0  C18_1  C18_2  C18_3
66  15.8     2     6    62.3   10.2    0.6

Let´s compare with the summary


> summary(fattyac_msc)
        C16_0           C16_1                 
 Min.   : 0.00   Min.   :1.500     
 1st Qu.:20.10   1st Qu.:2.000     
 Median :21.00   Median :2.200     
 Mean   :21.34   Mean   :2.267     
 3rd Qu.:22.90   3rd Qu.:2.500     
 Max.   :26.00   Max.   :3.500    
   
  C18_0            C18_1     
Min.   : 5.800   Min.   :43.80 
1st Qu.: 8.600   1st Qu.:51.95 
Median : 9.400   Median :54.50 
Mean   : 9.711   Mean   :53.93 
3rd Qu.:10.500   3rd Qu.:56.15 
Max.   :14.000   Max.   :62.30 

 C18_2            C18_3      
 Min.   : 5.500   Min.   :0.3000 
 1st Qu.: 7.600   1st Qu.:0.5000 
 Median : 8.500   Median :0.6000 
 Mean   : 8.503   Mean   :0.6032 
 3rd Qu.: 9.100   3rd Qu.:0.7000 
 Max.   :14.700   Max.   :1.3000  


Sample 66 has the higher value for C18:1 (oleic acid), but it is not isolated in the histogram. For some reasons this sample differs from the others especially from 100 to 1050 nm. We will wait forward to take a decision about this sample.
Until now we have been managing with the X matrix.
Now we start to study the Y matrix. First thing to do is to have a look to the summary, and of course to the histograms.
If you want to follow this tutorial, please send me an e_mail. I´ll send you the “txt” file attached.

> hist(C16_0,col="red")
> hist(C16_1,col="blue")
> hist(C18_0,col="green")
> hist(C18_1,col="brown")
> hist(C18_2,col="violet")
> hist(C18_3,col="orange")


No hay comentarios:

Publicar un comentario