Problème fonction interp()

[Louis Burguet]

Bonjour à tous,

J'ai déjà utilisé plusieurs fois cette fonction dans le cadre de bathymétrie, mais j'ai un problème avec mon jeu de données, je ne vois pas ce qui cloche...

Code : Tout sélectionner

library(akima)
bathy<-structure(list(X = c(503.2303, 503.2303, 501.7727, 502.7517,
503.5129, 503.6458, 504.314, 504.7277, 505.1559, 505.8524), Y = c(500.6256,
500.6256, 500.3125, 500.5406, 500.6899, 500.7457, 500.8613, 500.7345,
500.6057, 500.579), Z = c(49.3431, 49.3431, 49.4225, 49.4655,
49.3434, 48.487, 48.5458, 48.4216, 48.9574, 49.7252)), .Names = c("X",
"Y", "Z"), row.names = c(NA, 10L), class = "data.frame")
interp_bathy<-interp(bathy[,1],bathy[,2],bathy[,3])


Je reçois ce message d'erreur :
Error in interp.old(x, y, z, xo = xo, yo = yo, ncp = 0, extrap = extrap, :
duplicate data points: need to set 'duplicate = ..'

J'ai testé en changeant mon jeu de données par un autre, ça fonctionne....

Merci d’avance pour votre aide.

[Serge Rapenne]

Bonjour,

les 2 1er enregistrements sont identiques :

Code : Tout sélectionner

bathy[1:2,]
          X        Y       Z
1  503.2303 500.6256 49.3431
2  503.2303 500.6256 49.3431

et comme par défaut interp utilise l'argument duplicate="error", le code sort un erreur.

Soit tu supprimes l'une des 2 lignes, soit tu indiques duplicate="strip" dans les arguments de interp.

Code : Tout sélectionner

?interp

pour les détails

Serge

[Louis Burguet]

Merci beaucoup!