2012-07-17 38 views
11

I necesitan para determinar los valores críticos T para colas de un solo lado de 75% y 99%, para 40 grados de libertad.valores t críticos en R

El siguiente es el código para un de dos caras 99% los valores críticos t:

qt(0.01, 40) 

pero ¿cómo puedo determinar por un valor t crítico de un solo lado?

+7

qué te hace pensar que 'qt (0,01, 40)' es el valor crítico para la prueba de dos caras? Yo sugeriría leer '? Qt', y después de pensar un poco más acerca de lo uno y dos caras pruebas significa. –

+1

Realmente se trata de una cuestión en la comprensión de lo que esos valores críticos significan realmente. Esta es más una pregunta estadística y probablemente debería migrarse al sitio stackexchange de estadísticas. – Dason

Respuesta

20

El código que envió da el valor crítico para una prueba unilateral (ver here Por lo tanto la respuesta a usted pregunta es simplemente:..

abs(qt(0.25, 40)) # 75% confidence, 1 sided (same as qt(0.75, 40)) 
abs(qt(0.01, 40)) # 99% confidence, 1 sided (same as qt(0.99, 40)) 

Tenga en cuenta que la distribución t es simétrica Para un 2 prueba -sided (decir con confianza del 99%) se puede utilizar el valor crítico

abs(qt(0.01/2, 40)) # 99% confidence, 2 sided 
4

los comentarios de Josh son acertadas. Si no está muy familiarizado con los valores críticos me gustaría sugerir jugar con qt, leer el manual (?qt) junto con mirar una tabla de consulta (LINK). Cuando me mudé de SPSS a RI creé una función que hacía que el valor t crítico fuera bastante fácil (nunca usaría esto ahora, ya que lleva demasiado tiempo y con los valores p que generalmente se proporcionan en la salida es un punto discutible) Aquí está el código para eso:

critical.t <- function(){ 
    cat("\n","\bEnter Alpha Level","\n") 
    alpha<-scan(n=1,what = double(0),quiet=T) 
    cat("\n","\b1 Tailed or 2 Tailed:\nEnter either 1 or 2","\n") 
    tt <- scan(n=1,what = double(0),quiet=T) 
    cat("\n","\bEnter Number of Observations","\n") 
    n <- scan(n=1,what = double(0),quiet=T) 
    cat("\n\nCritical Value =",qt(1-(alpha/tt), n-2), "\n") 
} 

critical.t() 
+1

n-2 como sus grados de libertad con ninguna opción para cambiarlo? ¿Parece que solo lo usaste para una regresión lineal simple? – Dason

+0

del punto en que era Estadísticas I. Pero eso suena como que el cartel está en demasiado. –

3

Extendiendo @Ryogi respuesta anterior, se puede aprovechar el parámetro lower.tail así:

qt(0.25/2, 40, lower.tail = FALSE) # 75% de confianza

qt(0.01/2, 40, lower.tail = FALSE) # 99% de confianza

Cuestiones relacionadas