我用ropls包畫PCA, PLS DA, OPLS DA圖,但是數(shù)據(jù)導(dǎo)入出了些問(wèn)題。
> data
sample calss X.2.methoxyethoxy.propanoic.acid X.gamma.Glu.Leu.Ile
1? asss1? ? 1? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 6.8? ? ? ? ? ? ? ? 5.9
2? asss2? ? 1? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 3.5? ? ? ? ? ? ? ? 3.6
3 assss3? ? 1? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 3.6? ? ? ? ? ? ? ? 3.8
4? asss4? ? 1? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 2.8? ? ? ? ? ? ? ? 3.9
5 assdd4? ? 2? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 4.7? ? ? ? ? ? ? ? 4.1
6? sdv6? ? 2? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 4.6? ? ? ? ? ? ? ? 4.2
7? ghj7? ? 2? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 4.5? ? ? ? ? ? ? ? 4.3
8? fdg9? ? 2? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 4.4? ? ? ? ? ? ? ? 5.1
我的數(shù)據(jù)導(dǎo)入后是這樣的册倒,但是運(yùn)行畫圖命令后出來(lái)了個(gè)error:
> pca <- opls(data)
Error: 'x' data frame must contain columns of 'numeric' vectors only
然后我就把first and second column刪掉(因?yàn)榈诙袛?shù)據(jù)沒(méi)用),沒(méi)想到出現(xiàn)一個(gè)
warning:
> data.exp <- data[,c(-1, -2)]
> pca <- opls(data.exp)
Warning message:
Single component model: only 'overview' and 'permutation' (in case of
single response (O)PLS(-DA)) plots available
也順帶出來(lái)了一些數(shù)據(jù)和一個(gè)圖抑党,結(jié)果這個(gè)圖是空白匙铡,只有橫縱坐標(biāo)
> pca
PCA
8 samples x 3 variables
standard scaling of predictors
R2X(cum) pre ort
Total? ? 0.704? 1? 0
我問(wèn)作者是不是數(shù)據(jù)有什么格式要求,就按照上面的格式給作者郵件伯顶,作者回復(fù):
The model tells you that one component is sufficient (ie the second one is
not bringing useful information).
If, however, you want to compute additional components (ie 2 components to
display the score plot), you can set the predI argument (eg to 2 or more
instead of the NA default value).
他說(shuō)一個(gè)成分足夠了囚灼,但是我沒(méi)有得到圖怎么辦呢?加主成分這個(gè)建議在SIMCA軟件里面也
有祭衩,不過(guò)我用別的包(muma包)不需要加主成分做出的PCA圖有3個(gè)主成分灶体,只是只能出圖,
沒(méi)有上面那些R2X數(shù)據(jù)掐暮。
我該怎么修改我的數(shù)據(jù)格式好呢蝎抽?如有描述不足或者不當(dāng)請(qǐng)多多指教