Création d'une ellipse de distribution normale bivariée sur un tracé ayant des axes x et y log-transformés

Question

Bien que j'ai trouvé de nombreuses réponses sur la façon de dessiner une ellipse de distribution normale bivariée sur un tracé existant, j'ai une question concernant le tracé de cette ellipse sur un existant tracé pour lequel les axes x et y sont log-transformés.

Comme exemple, je l'ai données suivantes auxquelles j'ajouter une ellipse

library(mixtools) 
library(truncnorm) 
x<-rtruncnorm(n=100, a=0, b=20) 
y=1+.3*x+.3*rnorm(100) 
data<-cbind(x,y) 
mu <-c(mean(x), mean(y)) 
sigma <- var(data) 
plot(data) 
ellipse(mu, sigma, alpha=0.1, npoints = 200, newplot = FALSE) 

Cependant, mes données réelles exige que j'utiliser le log transformé x et y axe, comme si

plot(data,log="xy") 

Lors du traçage de la fonction "ellipse", je ne reçois une ellipse

ellipse(mu, sigma, alpha=0.1, npoints = 200, newplot = FALSE) 

Ajout "log" aux spécifications de la fonction ellipse est également aucune option

ellipse(mu, sigma, alpha=0.1, npoints = 200, newplot = FALSE,log="xy") 

message d'avertissement: En plot.xy (xy.coords (x, y), type = type, ...): "log" n'est pas un paramètre graphique.

Quelqu'un peut-il m'aider avec celui-ci?

Merci!

Answer 1

Le module car peut gérer des ellipses dans un espace log-transformé. Voici un avant-goût de ce qu'il peut faire. Une ligne est sans l'argument log alors que l'ellipse presque complète est avec l'argument log.

library(car) 
plot.default(data, log = "xy") 
car::ellipse(center = log(mu), shape = sigma, radius = 3, center.pch = NULL)