6
C'è un modo per aggiungere intervalli di confidenza ad un qqplot?Aggiunta di intervalli di confidenza ad un grafico qq?
Ho un set di dati di valori di espressione genica, che ho osservata con un PCA:
pca1 = prcomp(data, scale. = TRUE)
sto ora alla ricerca di valori anomali controllando la distribuzione dei dati contro la distribuzione normale attraverso:
qqnorm(pca1$x,pch = 20, col = c(rep("red", 73), rep("blue", 33)))
qqline(pca1$x)
Questo è il mio dati:
data = [2.48 104 4.25 219 0,682 0,302 1,09 0,586 90,7 344 13,8 1,17 305 2,8 79,7 3,18 109 0,932 562 0,958 1.87 0.59 114 391 13.5 1.41 208 2.37 166 3.42]
avrei ora piace tracciare intervalli di confidenza al 95% per verificare quali sono i punti dati all'esterno. Qualche consiglio su come farlo?
Così si vuole sottrarre la tua distribuzione campione dalla distribuzione normale teorica? Sembra che quello che vuoi fare sia usare 'nls' per adattare i tuoi dati a una normale funzione dist e catturare i dati di confidenza dall'output di' nls'. –
È molto più probabile che tu riceva una risposta utile se fornisci un [set di dati minimo e riproducibile] (http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example/5963610 # 5963610) insieme al codice che hai provato. Grazie! – Henrik
Ho modificato il mio post iniziale con alcuni dati. È possibile acquisire i dati di confidenza dall'output di qqnorm? – user2846211