To jest tylko wersja do druku, aby zobaczyć pełną wersję tematu, kliknij TUTAJ
Statystyka, prognozowanie, ekonometria, data mining
Forum miłośników statystyki - Portal Statystyczny

Biblioteki R, Pakiety R - [R] Zadania, zadania z uczelni, ćwiczeń, ogólnie

pRoblematyczna - 2013-12-03, 18:49

Prawie, ze z bolem mojego humanistycznego serca i rozumu doszlam do tego...tzn. teraz nie wiem jak podstawic moj zbior danych.
szw1710 - 2013-12-03, 19:30

Proponuję dokładnie opisać problem. Wtedy można konkretnie pomóc. Opisz co masz zrobić, podaj swoje dane itp.

Szkoda, że ludzie przywykli dzielić się na ścisłowców i humanistów. Humanistą jest ten, kto nie umie (nie cierpi) matematyki. A to nieprawda. Humanistyczną jest każda nauka mająca jakikolwiek związek z człowiekiem. Tak więc i matematyka jest nauką humanistyczną, a ja - zawodowy matematyk - za humanistę się uważam. Także z moją miłością do literatury.

MateuszEr - 2014-11-30, 07:04

Witam wszystkich,
moglby mi ktos pomoc z zadaniem?
Mam podane wektory o roznych nazwach (42 szt) z czego 40 ma dlugosc 100 a 2 wektory 152. Wszystkie wektory maja wartosci (0,1). Musze stworzyc matryce z tych wektorow z iloscia wierszy =100. Macie jakies pomysly?
pozdrawiam

ania28 - 2015-08-12, 12:58

Witam,

Mam pewnie problem z wykonaniem testu Tukeya. Dwa lata temu pracowałam w ówczesnej wersji programu wg schematu:
Załadowałam zbiór danych, stworzyłam model liniowy (do jednoczynnikowej analizy wariancji, np A~B) -> Anova (II) i wpisywałam komendę HSD.test(NazwaModelu, trt="B") oczywiście przy włączonym pakiecie agricolae i analiza była wykonana.

Obecnie wyskakuje tylko komunitat: niestosowalna metoda dla 'TukeyHSD' zastosowana do obiektu klasy "lm"

Próbowałam obecnie także komendy znalezionej w internecie: TukeyHSD(LinearModel.1, trt="B", ordered = FALSE, conf.level=0.95) ale też wyskakuje to samo.

Czy ktoś wie jak w obecnej wersji R (3.2.1.) wykonać test Tukeya?

Z góry dziękuję za odpowiedź, pozdrawiam Ania

Crunchy - 2015-08-12, 13:50

Prawdopodobnie nie zaktualizowałaś pakietu.
Kod:
> HSD.test(mod,"virus", group=T, console=T)

Study: mod ~ "virus"

HSD Test for yield

Mean Square Error:  22.48917

virus,  means

      yield      std r  Min  Max
cc 24.40000 3.609709 3 21.7 28.5
fc 12.86667 2.159475 3 10.6 14.9
ff 36.33333 7.333030 3 28.0 41.8
oo 36.90000 4.300000 3 32.1 40.4

alpha: 0.05 ; Df Error: 8
Critical Value of Studentized Range: 4.52881

Honestly Significant Difference: 12.39967

Means with the same letter are not significantly different.

Groups, Treatments and means
a        oo      36.9
ab       ff      36.33
bc       cc      24.4
c        fc      12.87
> installed.packages(priority='NA')
             Package        LibPath                              Version  Priority
agricolae    "agricolae"    "C:/Program Files/R/R-3.2.1/library" "1.2-1"  NA     
coda         "coda"         "C:/Program Files/R/R-3.2.1/library" "0.17-1" NA     
combinat     "combinat"     "C:/Program Files/R/R-3.2.1/library" "0.0-8"  NA     
deldir       "deldir"       "C:/Program Files/R/R-3.2.1/library" "0.1-9"  NA     
klaR         "klaR"         "C:/Program Files/R/R-3.2.1/library" "0.6-12" NA     
LearnBayes   "LearnBayes"   "C:/Program Files/R/R-3.2.1/library" "2.15"   NA     
sp           "sp"           "C:/Program Files/R/R-3.2.1/library" "1.1-1"  NA     
spdep        "spdep"        "C:/Program Files/R/R-3.2.1/library" "0.5-88" NA     
translations "translations" "C:/Program Files/R/R-3.2.1/library" "3.2.1"  NA     
             Depends                                                 
agricolae    NA                                                     
coda         "R (>= 2.14.0)"                                         
combinat     NA                                                     
deldir       "R (>= 0.99)"                                           
klaR         "R (>= 2.10.0), MASS"                                   
LearnBayes   NA                                                     
sp           "R (>= 3.0.0)"                                         
spdep        "R (>= 3.0.0), methods, sp (>= 1.0), Matrix (>= 1.0.12)"
translations NA                                                     
             Imports                                                          LinkingTo
agricolae    "klaR, MASS, nlme, cluster, spdep"                               NA       
coda         "lattice"                                                        NA       
combinat     NA                                                               NA       
deldir       NA                                                               NA       
klaR         "combinat"                                                       NA       
LearnBayes   NA                                                               NA       
sp           "methods, graphics, utils, lattice, grid"                        NA       
spdep        "LearnBayes, deldir, boot (>= 1.3-1), splines, coda, nlme, MASS" NA       
translations NA                                                               NA       
             Suggests                                                                                                     
agricolae    NA                                                                                                           
coda         NA                                                                                                           
combinat     NA                                                                                                           
deldir       "polyclip"                                                                                                   
klaR         "scatterplot3d (>= 0.3-22), som, mlbench, rpart, e1071"                                                       
LearnBayes   NA                                                                                                           
sp           "RColorBrewer, rgdal (>= 0.8-7), rgeos (>= 0.2-20), gstat"                                                   
spdep        "parallel, spam(>= 0.13-1), RANN, RColorBrewer, lattice,\nxtable, maptools (>= 0.5-4), foreign, igraph, knitr"
translations NA                                                                                                           
             Enhances License           License_is_FOSS License_restricts_use OS_type MD5sum
agricolae    NA       "GPL"             NA              NA                    NA      NA   
coda         NA       "GPL (>= 2)"      NA              NA                    NA      NA   
combinat     NA       "GPL-2"           NA              NA                    NA      NA   
deldir       NA       "GPL (>= 2)"      NA              NA                    NA      NA   
klaR         NA       "GPL-2"           NA              NA                    NA      NA   
LearnBayes   NA       "GPL (>= 2)"      NA              NA                    NA      NA   
sp           NA       "GPL (>= 2)"      NA              NA                    NA      NA   
spdep        NA       "GPL (>= 2)"      NA              NA                    NA      NA   
translations NA       "Part of R 3.2.1" NA              NA                    NA      NA   
             NeedsCompilation Built 
agricolae    "no"             "3.2.1"
coda         "no"             "3.2.1"
combinat     NA               "3.2.1"
deldir       "yes"            "3.2.1"
klaR         "no"             "3.2.1"
LearnBayes   "no"             "3.2.1"
sp           "yes"            "3.2.1"
spdep        "yes"            "3.2.1"
translations NA               "3.2.1"

ania28 - 2015-08-12, 14:03

Dziękuję bardzo za pomoc! Działa :-)
macko - 2016-04-19, 14:13

Cześć, mam jedno pytanie.

Jak powinien wyglądać kod, który wyświetli mi skorygowane R-kwadrat w prostym modelu regresji liniowej? Wyświetla strukturę modelu za pomocą funkcji str(), ale nie widzę tak odpowiedniego odwołania.

ania28 - 2018-02-23, 19:03

Witam,

Mam problem z wykonaniem testu Dunn'a w pakiecie R. Zainstalowałam dzisiaj nową wersję (3.4.3) i cały czas wyskakuje mi błąd, chociaż załadowałam pakiet. Czy ktoś mógłby podać mi jaka jest komenda do obecnej wersji?

Z góry dziękuję za pomoc,
Ania

mikos - 2018-02-23, 19:37

* z jakiego pakietu korzystasz ?

* np. pod linuksem pakiet: DescTools -- funkcja: DunnTest działa bez problemów.

ania28 - 2018-02-23, 20:43

Pakiety zaktualizowałam, pakiet DescTools załadowałam, wyświetla się błąd dla formuły jaką stworzyłam (na podstawie jakiś tam danych z internetu) "brakuje argumentu 'g', a nie ma określonej wartości domyślnej".
Ania

mikos - 2018-02-23, 20:57

U mnie działa taki zapis:

Kod:
y <- rnorm(90)
f <- factor(rep(letters[1:3],each=30))


Kod:
DunnTest(y,f)

 Dunn's test of multiple comparisons using rank sums : holm 

    mean.rank.diff   pval   
b-a     -10.366667 0.3730   
c-a      -4.533333 0.7743   
c-b       5.833333 0.7743   
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

igrek77 - 2018-03-29, 17:01

Witam. Potrzebuje pomocy w rozwiazaniu zadania w programie R.Dla mnie jest to czarna magia, ale dla osoby ktora sie na tym zna podejrzewam ze to kwestia pol godziny lub godziny. Zadanie polega na napisaniu kodu do wykresu ktory dostalem. Zadanie jest w jezyku angielskim.
pozdrawiam



Powered by phpBB modified by Przemo © 2003 phpBB Group