2017-03-14 35 views
1

私はRで初心者ですMomocsパッケージでアウトライン解析を行うために、楕円フーリエ係数のシーケンスをRにロードする方法は?

私は楕円フーリエ係数(高調波)のあらかじめ計算されたシーケンスを持つテキストファイルを持っています。最初の列はオブジェクト名を保持し、他の高調波はスペースで分けられます:A0 B0 C0 D0 ... An Bn Cn Dn。このファイルの一部です:

Lceph-sp1-1 0 0,27 -1 0 -0,29 0,2 0,1 -0,05 0,15 0,1 0,07 0,02 
Lceph-sp1-2 0 -0,36 1 0 -0,27 0,25 0,12 -0,09 -0,26 -0,03 -0,17 0 
Lceph-sp1-3 0 0,25 -1 0 -0,29 0,19 0,09 -0,05 0,15 0,1 0,06 0,01 
Lceph-sp1-4 0 -0,37 1 0 -0,27 0,26 0,09 -0,07 -0,19 -0,07 -0,12 -0,02 

私はread.delim2コマンドを使用して、Rにファイルをロードします。

ef <- read.delim2("filename", header=FALSE, sep="") 

その後、私はようMomocs packageから関数を用いて得られたデータを、分析したいです較正高調波電力,プロット,PCAおよびその他。

これを行うには、ロードされたデータをCoe(?)または他のオブジェクトに変換する必要があります。 (私には正確には分かりません)

Momocs Rパッケージの分析のためにロードされたデータをどのように準備するのですか?あなたがR/Momocs外の楕円フーリエ係数を計算する場合は、次のよう

答えて

0

、あなたはOutCoeオブジェクトとしてそれらをインポートする必要があります。

> x 
An OutCoe object [ elliptical Fourier analysis ] 
-------------------- 
    - $coe: 4 outlines described, 3 harmonics 
- $fac: 3 classifiers: 
    'ind' (factor 1): Lceph. 
'sp' (factor 1): sp1. 
'id' (factor 4): 1, 2, 3, 4. 

# we load Momocs and import your table. 
library(Momocs) 

#Here I use your data (using `text` but a path to your file will work by your side 
data <- read.table(dec=",", text=" 
Lceph-sp1-1 0 0,27 -1 0 -0,29 0,2 0,1 -0,05 0,15 0,1 0,07 0,02 
Lceph-sp1-2 0 -0,36 1 0 -0,27 0,25 0,12 -0,09 -0,26 -0,03 -0,17 0 
Lceph-sp1-3 0 0,25 -1 0 -0,29 0,19 0,09 -0,05 0,15 0,1 0,06 0,01 
Lceph-sp1-4 0 -0,37 1 0 -0,27 0,26 0,09 -0,07 -0,19 -0,07 -0,12 -0,02") 

# we can extract the `$fac` cofactor/variates table with `lf_structure` 
fac <- lf_structure(data[, 1], names=c("ind", "sp", "id"), split="-") 

# then we simply pass all columns but the first, the `fac` 
# and two others arguments `Momocs` cannot deduce `method` and `norm` 
# see ?Outcoe 
x <- OutCoe(coe=data[, -1], fac=fac, method="efourier", norm=TRUE) 

現在 '作業' Outcoeオブジェクトを持っています

分析を続けるために使用できるもの:

x %>% PCA %>% plot 

それはあなたが探していたものですか?

+0

遅れて申し訳ありませんが、以前これを見たことがありません。 –

関連する問題