Quanto è utile/interessante questa discussione:
| Autore |
Discussione |
|
|
gianpaolo
Nuovo Arrivato
15 Messaggi |
 Inserito il - 10 maggio 2010 : 17:54:26
|
Salve a tutti, ho un bel problema da porvi.
digitando i comandi:
> x=runif(20000,min=0,max=1)
> hist(x,breaks=20)
così facendo mi genero un hist con i 20000 numeri casuali, ed ogni classe contiene un valore medio di circa 1000 eventi. A questo punto a me servirebbe fare un istogramma di queste frequenze. In teoria dovrebbe uscir fuori una gaussiana centrata sul valore 1000, con tutte le altre classi di freq (per es.) 990, 980, 970....alla sinistra di 1000, e 1010, 1020, 1030 ...ecc alla destra di 1000. Il mio problema è proprio questo. Come faccio ad ottenere un istogramma del genere?????
Spero di essere stato un po chiaro.
Grazie...
|
|
|
|
|
dallolio_gm
Moderatore
Prov.: Bo!
Città: Barcelona/Bologna
2445 Messaggi |
Inserito il - 10 maggio 2010 : 18:26:50
|
non ho capito :-)
x=runif(20000,min=0,max=1) questo ti crea un vettore di numeri casuali compresi tra 0 e 1
Citazione:
mi genero un hist con i 20000 numeri casuali, ed ogni classe contiene un valore medio di circa 1000 eventi.
il valore medio é di 1000 perché hai diviso l'istogramma in 20 blocchi... se fossero stati 10, la media sarebbe 2000.
Citazione:
In teoria dovrebbe uscir fuori una gaussiana centrata sul valore 1000,
é questo che non capisco.. il 1000 a che cosa corrisponderebbe? che significato avrebbe? |
Il mio blog di bioinformatics (inglese): BioinfoBlog
Sono un po' lento a rispondere, posso tardare anche qualche giorno... ma abbiate fede! :-) |
 |
|
|
gianpaolo
Nuovo Arrivato
15 Messaggi |
Inserito il - 10 maggio 2010 : 19:03:15
|
Ciao, capisco di non essere stato molto chiaro. Il discorso del 1000 è appunto come dici tu, 20000/20=1000, cioè ho una media di 1000 numeri per ogni classe. Sulla l'asse delle x come si vede dall'istogramma ho 20 colonne con una frequenza media di 1000 conteggi. Ad es. avrò delle colonne di 980, 1100, 960, 990, 1080, 1060, ecc....
Adesso io vorrei un istogramma in cui lungo l'asse x ci siano le frequenze del primo istogramma ottenuto attraverso i comandi che ho scritto sopra. Quindi essendo che tutte quelle colonne oscillavano intorno al valore 1000, il mio nuovo istogramma dovrebbe essere centrato intorno al numero 1000, ottenendo così una sorta di gaussiana. Alla destra della colonna centrale (1000) dovrebbero starci le colonne 980, 1100, 960, 990, 1080, 1060, ecc..
Io non capisco come potrei ottenere ques'ultimo tipo di istogramma.
spero il discorso sia un tantino più chiaro.
grazie |
|
|
|
gianpaolo
Nuovo Arrivato
15 Messaggi |
Inserito il - 10 maggio 2010 : 19:05:53
|
SCUSATE, volevo scrivere :
Alla DESTRA e SINISTRA della colonna centrale (1000) dovrebbero starci le colonne 980, 1100, 960, 990, 1080, 1060, ecc.. |
|
|
|
chick80
Moderatore
Città: Edinburgh
11491 Messaggi |
Inserito il - 10 maggio 2010 : 19:44:57
|
Se vuoi vedere qualcosa che assomigli di più ad una gaussiana devi fare più di 20 breaks nel primo istogramma.
Ad es
x <- runif(50000, min=0, max=1)
par(mfrow=c(1,2))
h <- hist(x, 1000)
hist(h$counts, 100)
|
Sei un nuovo arrivato?
Leggi il regolamento del forum e presentati qui
My photo portfolio (now on G+!) |
|
|
|
gianpaolo
Nuovo Arrivato
15 Messaggi |
Inserito il - 10 maggio 2010 : 20:01:11
|
BINGO!!!!!!!!!!
grazie 1000 chick80, mi hai risolto il problema. Sto da poco usando R ed ancora ho molte difficoltà.
Sono circa 3 gg che cercavo di risolvere questo punto.
Grazie!!!! |
|
|
|
TMax
Utente Junior
Prov.: BG
Città: Capriate
270 Messaggi |
Inserito il - 11 maggio 2010 : 16:45:33
|
cosi?
x<-rnorm(10000, 1000, 100)
hist(x)
|
|
|
|
chick80
Moderatore
Città: Edinburgh
11491 Messaggi |
|
| |
Discussione |
|
|
|
Quanto è utile/interessante questa discussione:
| MolecularLab.it |
© 2003-18 MolecularLab.it |
|
|
|
Istogrammi delle frequenze con R
Guide tools online
