Forum

Nome Utente:
Password:
Riconoscimi automaticamente
 Tutti i Forum
 MolecularLab
 Bioinformatica
 Aiuto per esercizio di statistica
 Nuova Discussione  Nuovo Sondaggio Nuovo Sondaggio
 Rispondi Aggiungi ai Preferiti Aggiungi ai Preferiti
Cerca nelle discussioni
I seguenti utenti stanno leggendo questo Forum Qui c'è:

Tag   esercizio statistica  Aggiungi Tag

Quanto è utile/interessante questa discussione:

Autore Discussione  

Pandorina
Nuovo Arrivato

Prov.: Na
Città: Napoli


4 Messaggi

Inserito il - 01 gennaio 2015 : 17:10:29  Mostra Profilo  Visita l'Homepage di Pandorina Invia a Pandorina un Messaggio Privato  Rispondi Quotando
Salve a tutti!
Ho un problema nella risoluzione di un esercizio di biostatistica, condivido con voi sperando possiate aiutarmi a risolverlo

L'esercizio pone un segmento di lunghezza 1 che viene rotto in maniera casuale due volte. Con l'espressione
(-1)^g+1 (matrice n+1 | g) (1-gu)^n


dove g rappresenta la lunghezza del segmento nell'intervallo (0,1).

Dopo la risoluzione l'esercizio chiede di scrivere un programma in R che simuli la procedura e di ricavarne l'istogramma.

Qualcuno sa darmi una manina?


Colgo l'occasione anche per fare un'altra domandina: è possibile salvare un'intera sessione di lavoro con R su ubuntu? Ossia posso salvare un file che contenga tutto ciò che scrivo su R da terminale?

Grazie ancora a tutti

chick80
Moderatore

DNA

Città: Edinburgh


11491 Messaggi

Inserito il - 01 gennaio 2015 : 18:11:14  Mostra Profilo  Visita l'Homepage di chick80 Invia a chick80 un Messaggio Privato  Rispondi Quotando
Citazione:
è possibile salvare un'intera sessione di lavoro con R su ubunt


Sì, puoi usare il comando savesession o, più semplicemente, consiglio di installare RStudio che è di gran lunga la miglior GUI per R disponibile. RStudio ti chiederà se salvare la sessione all'uscita.

In ogni caso le sessioni interattive si usano per dei test veloci, se vuoi avere delle analisi riproducibili salva degli script R, dove potrai anche aggiungere commenti etc.

Per il programma cosa non riesci a fare? Che istogramma vuoi ricavare? Cosa sarebbe la matrice nella formula? Cos'è gu?

Sei un nuovo arrivato?
Leggi il regolamento del forum e presentati qui

My photo portfolio (now on G+!)
Torna all'inizio della Pagina

Pandorina
Nuovo Arrivato

Prov.: Na
Città: Napoli


4 Messaggi

Inserito il - 01 gennaio 2015 : 20:31:25  Mostra Profilo  Visita l'Homepage di Pandorina Invia a Pandorina un Messaggio Privato  Rispondi Quotando
In primis: grazie per la risposta, proverò subito a scaricare sia RStudio che ad usare il comando!

Il programma che devo scrivere deve simulare la procedura di risoluzione dell'esercizio (ossia dividere il segmento di lunghezza 1 in 2 punti casuali e dimostrare che la probabilità che i frammenti non formino un triangolo è 3/4).

La formula che devo usare per questa dimostrazione è: http://i.imgur.com/Nu0SJfr.png ed infine elaborare (sempre tramite R) un istogramma della distribuzione del quoziente (sempre relativo all'esercizio)e produrre la stima della distribuzione.

I parametri della formula sono: u (viene assegnato dall'esercizio come 1/2) n ossia il numero di punti ed infine g che rappresenta la lunghezza del segmento.
Torna all'inizio della Pagina

chick80
Moderatore

DNA

Città: Edinburgh


11491 Messaggi

Inserito il - 03 gennaio 2015 : 12:20:19  Mostra Profilo  Visita l'Homepage di chick80 Invia a chick80 un Messaggio Privato  Rispondi Quotando
Hmmm... non mi è ancora del tutto chiaro il problema... (es. cos'è nella formula la sommatoria va da 0 a hu, cos'è hu?)

In ogni caso, l'implementazione della formula in R è triviale, ti faccio un esempio di come ottenere l'istogramma. Consiglio di usare la funzione replicate, che semplifica questo tipo di procedure.

Esempio banale: vogliamo calcolare l'istogramma dei risultati del lancio di n dadi


somma.dadi <- function(num.dadi)
  {
  # "Tiriamo il dado" con la funzione sample che, in questo caso sceglie
  # n numeri fra 1 e 6.  
  return (sum(sample(1:6, num.dadi, replace=T)))
  }

n.dadi <- 2
# Chiamiamo la funzione 10000 volte e mettiamo il risultato nel vettore res
res <- replicate(10000, somma.dadi(n.dadi))
# Creiamo l'istogramma
hist(res, 1:(6*n.dadi), col="black", las=1, freq=F)


Come atteso la distribuzione del lancio di 2 dadi sarà una gaussiana centrata sul 6

Sei un nuovo arrivato?
Leggi il regolamento del forum e presentati qui

My photo portfolio (now on G+!)
Torna all'inizio della Pagina

Pandorina
Nuovo Arrivato

Prov.: Na
Città: Napoli


4 Messaggi

Inserito il - 03 gennaio 2015 : 15:43:07  Mostra Profilo  Visita l'Homepage di Pandorina Invia a Pandorina un Messaggio Privato  Rispondi Quotando
Grazie mille dell'aiuto, proverò subito!

Se ti amano sei bravo...se ti odiano sei il migliore...
Torna all'inizio della Pagina
  Discussione  

Quanto è utile/interessante questa discussione:

 Nuova Discussione  Nuovo Sondaggio Nuovo Sondaggio
 Rispondi Aggiungi ai Preferiti Aggiungi ai Preferiti
Cerca nelle discussioni
Vai a:
MolecularLab.it © 2003-18 MolecularLab.it Torna all'inizio della Pagina