campionaria.theta2: Distribuzione campionaria

In masspastore/ADati: Analisi dei dati in Psicologia

Distribuzione campionaria

Description

Produce la distribuzione campionaria di una statistica theta basata su due campioni.

Usage

campionaria.theta2( Omega1, n1 = 2, n2 = n1, Omega2 = NULL, 
      theta = 'diff', replace = FALSE, exact = FALSE, 
      grafico = TRUE, B = 1000, parziali = FALSE )

Arguments

Omega1	vettore numerico, contiene i valori della popolazione 1.
n1	numerosita' dei campioni estratti da Omega1.
n2	numerosita' dei campioni estratti da Omega2. Se non specificato e' uguale a n1.
Omega2	vettore numerico, contiene i valori della popolazione 2. Se non specificato e' uguale ad Omega1.
theta	funzione statistica a due argomenti, di cui produrre la distribuzione campionaria, per default diff.
replace	logico, indica se i campioni sono con reinserimento (TRUE) oppure senza reinserimento (FALSE).
exact	logico, indica se creare tutti i campioni possibili (TRUE) oppure solo un numero B di campioni estratti casualmente (FALSE).
grafico	logico, indica se produrre il grafico (TRUE) oppure no (FALSE).
B	numero di campioni casuali da estrarre se exact=FALSE.
parziali	logico, indica se produrre i grafici parziali durante il ricampionamento.

Details

Richiede la funzione combinations del pacchetto gtools.

ATTENZIONE: dato che il numero di campioni possibili aumenta con l'aumento della dimensione di Omega, quando tale numero supera 50000 la distribuzione viene prodotta per approssimazione campionando per B volte.

Stampa una tabella con i valori della popolazione e le statistiche della distribuzione campionaria.

Value

Restituisce una lista con i seguenti elementi:

$Omega1	i valori della popolazione 1.
$Omega2	i valori della popolazione 2.
$Stheta	distribuzione campionaria della statistica theta.
$theta	funzione statistica passata in input.

Author(s)

Massimiliano Pastore

Examples

## dati uniformi
require(gtools)
campionaria.theta2(0:3,theta='max',exact=TRUE,replace=TRUE) # con reinserimento
campionaria.theta2(0:3,theta='max',exact=TRUE) # senza reinserimento

## dati normali
fx <- function(x1,x2){max(mean(x1),mean(x2))}
campionaria.theta2(rnorm(50),n1=8,theta='fx')
campionaria.theta2(rnorm(50),n1=8,Omega2=rnorm(50,5,3),theta='fx')

## dati non normali
campionaria.theta2(rchisq(50,1),n1=8,Omega2=rnorm(50,5,3),theta='fx')

masspastore/ADati documentation built on April 18, 2024, 5:55 p.m.