1 / 8

a lkfos PCA

a lkfos PCA. SPH 247 Statistical Analysis of Laboratory Data. > library( ISwR ) > data( alkfos ) > names( alkfos ) [1] " grp " "c0" "c3" "c6" "c9" "c12" "c18" " c24“ > summary( alkfos ) grp c0 c3 c6

autumn-burke
Télécharger la présentation

a lkfos PCA

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. alkfos PCA SPH 247 Statistical Analysis of Laboratory Data

  2. > library(ISwR) > data(alkfos) > names(alkfos) [1] "grp" "c0" "c3" "c6" "c9" "c12" "c18" "c24“ > summary(alkfos) grp c0 c3 c6 Min. :1.000 Min. : 81.0 Min. : 79.0 Min. : 79.0 1st Qu.:1.000 1st Qu.:124.0 1st Qu.:110.0 1st Qu.:112.5 Median :1.000 Median :161.0 Median :134.0 Median :146.0 Mean :1.465 Mean :159.9 Mean :143.3 Mean :149.0 3rd Qu.:2.000 3rd Qu.:190.5 3rd Qu.:170.5 3rd Qu.:167.5 Max. :2.000 Max. :237.0 Max. :277.0 Max. :298.0 c9 c12 c18 c24 Min. : 74.0 Min. : 73.0 Min. : 69.0 Min. : 68.0 1st Qu.:107.5 1st Qu.:103.0 1st Qu.:114.0 1st Qu.:111.8 Median :140.5 Median :134.0 Median :140.0 Median :146.5 Mean :140.3 Mean :143.3 Mean :149.1 Mean :150.1 3rd Qu.:163.5 3rd Qu.:168.0 3rd Qu.:185.5 3rd Qu.:189.8 Max. :240.0 Max. :289.0 Max. :288.0 Max. :260.0 NA's : 1.0 NA's : 2.0 NA's : 4.0 NA's : 5.0 SPH 247 Statistical Analysis of Laboratory Data

  3. > af1 <- alkfos[complete.cases(alkfos),] > afmat <- af1[,-1] > grp <- af1$grp > pca.af <- prcomp(afmat,scale=T) > pca.tam <- prcomp(afmat[grp==2,],scale=T) > pca.placebo <- prcomp(afmat[grp==1,],scale=T) > plot(pca.af$x[,1:2],col=grp) > plot(pca.af) SPH 247 Statistical Analysis of Laboratory Data

  4. SPH 247 Statistical Analysis of Laboratory Data

  5. SPH 247 Statistical Analysis of Laboratory Data

  6. > library(MASS) > lda.af <- lda(afmat,grp) > lda.af Prior probabilities of groups: 1 2 0.6 0.4 Group means: c0 c3 c6 c9 c12 c18 c24 1 156.7143 161.8571 173.9048 158.4286 163.8571 164.3333 163.2857 2 164.2143 125.1429 123.7143 118.7143 117.2857 130.7857 134.8571 Coefficients of linear discriminants: LD1 c0 0.0618073455 c3 -0.0329471378 c6 0.0004421163 c9 -0.0232320119 c12 -0.0248954902 c18 0.0113410946 c24 0.0003473940 SPH 247 Statistical Analysis of Laboratory Data

  7. > af.pred<- predict(lda.af) > table(af.pred$class,grp) grp 1 2 1 20 0 2 1 14 > plot(grp,af.pred$x) SPH 247 Statistical Analysis of Laboratory Data

  8. SPH 247 Statistical Analysis of Laboratory Data

More Related