Jak obliczyć krytyczną wartość t za pomocą R?

13

Przepraszam, jeśli to nowe pytanie; Po raz pierwszy próbuję nauczyć się statystyki. Myślę, że mam podstawową procedurę, ale staram się ją wykonać z R.

Próbuję więc ocenić znaczenie współczynników regresji w wielokrotnej regresji liniowej formy

y^=Xβ^

Myślę, że statystyki t do testowania podajeH0:β^j=0,Ha:β^j0

t0=β^j0se(β^j)=β^jσ^2Cjj=β^jCjjSSRes/(np)
gdzie to przekątna pozycja w . j t h ( X X ) - 1Cjjjth(XX)1

Na razie w porządku. Wiem, jak obliczyć wszystkie te wartości za pomocą operacji macierzowych w R. Ale żeby odrzucić wartość zerową, książka mówi, że potrzebuję

|t0|>tα/2,np

Jak obliczyć tę wartość krytyczną przy użyciu R? tα/2,np W tej chwili jedynym sposobem, w jaki wiem, jak znaleźć te wartości, jest zajrzenie do tabeli na końcu książki. Musi być lepszy sposób.

Jurajski
źródło
Działa to funkcja qt ().
gość

Odpowiedzi:

22

Witamy w Stats Stackexchange. Możesz obliczyć wartość krytyczną w R przez wykonanie .tα/2,npqt(1-alpha/2, n-p)

W poniższym przykładzie proszę R o podanie mi wartości krytycznej dla . Wynikiem jest lista pierwszych dziesięciu wartości krytycznych dla rozkładu t na danym poziomie ufności:d f = 1 , 2 , , 1095%df=1,2,,10
> qt(.975, 1:10)
[1] 12.706205 4.302653 3.182446 2.776445 2.570582 2.446912 2.364624 2.306004 2.262157 2.228139

Christopher Aden
źródło