Skip to content

Inconsistant results with adegenet::spca and adespatial::multispati #355

@arthurdegrandpre

Description

@arthurdegrandpre

Hello,

I have been looking into solutions for spatial multivariate analysis and tried both spca and multispati, but their results are vastly different.
Isn't spca basically a wrapper to multispati? Maybe I'm mistaken in their use and/or call synthax.

Following is a reproducible example.
Merci!
Arthur

`
library(ade4)
library(adespatial)
library(adegraphics)
library(adegenet)
library(sp)
library(spdep)

data("mafragh")

analysis with multispati

mxy = as.matrix(mafragh$xy)
nb = chooseCN(coordinates(mxy), type = 1, plot.nb = FALSE)
lw = nb2listw(nb, style = 'W', zero.policy = TRUE)
pca = dudi.pca(mafragh$env, scale=T, center=T, scannf=F, nf=2)
ms = multispati(pca, lw, scannf=F)

ms_loadings = ms$c1[,1]^2

analysis with spca

spca = spca(mafragh$flo, xy=mafragh$xy, ask=F, type=1, scannf=F, scale=T, center=T, nfposi=2, nfnega=0)
spca_loadings = spca$c1[,1]^2

comparison

loadingplot(ms_loadings)
loadingplot(spca_loadings)
ms
spca
`

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions