Statistiques inverse avec R

Question

Ce que je veux faire des sons simples. Je veux tracer une courbe de QI normal avec R avec une moyenne de 100 et un écart-type de 15. Ensuite, je voudrais pouvoir superposer un diagramme de dispersion des données sur le dessus de celui-ci.

Quelqu'un sait comment faire cela?

Answer 1

Je devine ce que vous voulez faire est la suivante: vous voulez tracer la modèle densité normale avec une moyenne de 100 et sd = 15, et que vous voulez superposer au-dessus de ce que le densité empirique de certains ensemble d'observations qui suivent la densité de soi-disant modèle normal, de sorte que vous pouvez visualiser dans quelle mesure la densité du modèle correspond à la densité empirique. Le code ci-dessous devrait le faire (ici, x serait le vecteur d'observations réelles, mais à des fins d'illustration je générer avec une distribution mixte normale N (100,15) + 15 * N (0,1), à savoir la prétendue N (100,15) la distribution plus le bruit).

require(ggplot2)
x <- round( rnorm( 1000, 100, 15 )) + rnorm(1000)*15
dens.x <- density(x)
empir.df <- data.frame( type = 'empir', x = dens.x$x, density = dens.x$y )
norm.df <-  data.frame( type = 'normal', x = 50:150, density = dnorm(50:150,100,15))
df <- rbind(empir.df, norm.df)
m <- ggplot(data = df, aes(x,density))
m + geom_line( aes(linetype = type, colour = type))

Answer 2

Eh bien, il est plus comme un histogramme, car je pense que vous attendez ces d'être plus comme un entier processus arrondi:

x<-round(rnorm(1000, 100, 15))
y<-table(x)
plot(y)
par(new=TRUE)
plot(density(x), yaxt="n", ylab="", xlab="", xaxt="n")

Si vous voulez que la valeur théorique de dnorm superposées, puis utilisez l'une de ces:

lines(sort(x), dnorm(sort(x), 100, 15), col="red")

-ou

points(x, dnorm(x, 100, 15))

Answer 3

Vous pouvez générer des scores QI PDF avec:

curve(dnorm(x, 100, 15), 50, 150)

Mais pourquoi voulez-vous la diffusion sur superposition courbe de densité? À mon humble avis, qui est très inhabituel ...

Answer 4

En plus des autres bonnes réponses, vous pourriez être intéressé à tracer un certain nombre de panneaux, chacun avec son propre graphique. Quelque chose comme ça .