Statystyka, prognozowanie, ekonometria, data mining Strona Główna
Reklama pqstat.pl
Statystyka, prognozowanie, ekonometria, data mining
Forum miłośników statystyki - Portal Statystyczny

FAQFAQ  SzukajSzukaj  UżytkownicyUżytkownicy  GrupyGrupy  StatystykiStatystyki
RejestracjaRejestracja  ZalogujZaloguj  Chat   Regulamin  Kadra forum
PORTAL STATYSTYCZNY
 Ogłoszenie 
FORUM STATYSTYCZNE MA JUŻ 10 LAT

Znasz statystykę lub ekonometrię, metody prognozowania, data mining i chcesz pomóc w rozwoju forum statystycznego ?
Pisz na: administrator(małpa)statystycy.pl

Rozpoczął swoją działalność portal statystyczny - masz pomysł na jego rozwój ?

Drogi forumowiczu! Zanim napiszesz posta zapoznaj się z regulaminem forum i przedstaw się
The International Year of Statistics (Statistics2013) Smarter Poland Portal statystyczny
Znalezionych wyników: 29
Statystyka, prognozowanie, ekonometria, data mining Strona Główna
Autor Wiadomość
  Temat: grafika w R
statnowy

Odpowiedzi: 4
Wyświetleń: 1297

PostForum: Biblioteki R, Pakiety R   Wysłany: 2017-07-21, 13:21   Temat: grafika w R
Kod:

data <- data.frame(klm=c(1L,2L,3L),
               year=c(rep(2000,3),rep(2016,3)),
               Category1=c(c(1313,1057,981),c(1456,1200,1124)),
               Category2=c(c(1661,1257,1301),c(1827,1444,1600)))


data$klm <- factor(data$klm, levels=1:3, labels=c("Belgium", "Netherlands","Germany"))

data$year <- factor(data$year,levels=c(2000,2016))

Średnia=
  ggplot(data=data, aes(x=klm, y=Category2, fill=year)) +
  geom_bar(stat="identity", position=position_dodge(),colour="black")+
  scale_fill_manual(values = c("blue", "red"))+
  geom_text(aes(y = Category2,label = format(Category2,big.mark = " ")),position = position_dodge(width = 1),vjust = -0.5, hjust = 0.5, cex = 3)


Średnia+
  geom_bar(data=data, aes(x=klm, y=Category1, fill=year),stat="identity", width = 0.5,position=position_dodge(), colour="black")+
  geom_text(aes(y = Category1,label = format(Category1,big.mark = " ")),position = position_dodge(width = 1),vjust = -0.5, hjust = 0.5, cex = 3)


Jak zmienić kolory warstw i dorobić legendę taką jak w poście 1

Z góry dzięki za pomoc
  Temat: grafika w R
statnowy

Odpowiedzi: 4
Wyświetleń: 1297

PostForum: Biblioteki R, Pakiety R   Wysłany: 2017-07-04, 08:28   Temat: grafika w R
Dzięki wielkie ale kodu nie widać, możesz jeszcze raz przesłać

co do uwagi, że nie jest to wykres do naśladowania - no może i nie ale taki chce uzyskać

Wiesz może jak taki uzyskać


Miłego dnia
  Temat: grafika w R
statnowy

Odpowiedzi: 4
Wyświetleń: 1297

PostForum: Biblioteki R, Pakiety R   Wysłany: 2017-07-03, 15:31   Temat: grafika w R


jak wykonać taki wykres w R, Razem z legendą

Próbuję z pakietem ggplot2 ale marne wyniki
  Temat: [R] zmiana nazwy wybranej kolumny
statnowy

Odpowiedzi: 2
Wyświetleń: 2272

PostForum: Biblioteki R, Pakiety R   Wysłany: 2017-03-09, 07:59   Temat: [R] zmiana nazwy wybranej kolumny
Ogólnie możesz zmienić kilka wybranych kolumn

Kod:
colnames(dane)[c(1,3,7)]=cbind("kielich","lodyga","lisc")


podmienisz nazwy pierwszej trzeciej i siódmej kolumny
  Temat: Problem w wynikiem - niedokładność
statnowy

Odpowiedzi: 0
Wyświetleń: 826

PostForum: Biblioteki R, Pakiety R   Wysłany: 2017-03-08, 10:20   Temat: Problem w wynikiem - niedokładność
Witam

Mam taki problem pytanie

licząc
Kod:
sin(pi)/2
dostaję wynik 6.123032e-17
oczywiście wiem, że przy tym poziomie to jest bliskie zeru,
zastanawiam się czy przy poważniejszych operacjach nie zaburzy to
wyniku?

Jak się z tym uporać?
  Temat: Wypełnienie tabeli pod warunkiem
statnowy

Odpowiedzi: 0
Wyświetleń: 1126

PostForum: Biblioteki R, Pakiety R   Wysłany: 2017-01-16, 09:20   Temat: Wypełnienie tabeli pod warunkiem
Cześć

Mam długą ankietę i zrobiłem bazę każde pytanie osobna kolumna, w wierszach kolejne ankiety
w pytaniu 1 znajduje się pytanie czy respondent miał wycieczkę z organizatorem (przejazd) (1=tak,2=nie)
w pytaniu 2 znajduje się pytanie jaki cell (1=osobisty 2=zawodowy)
w pytaniu 3 znajduje się pytanie :sposób podroży (1=pieszo, 2 coś na benzynę, 3 kolej, samolot, 4 inne) - zawsze tylko jedna odpowiedz ma być - nie zwracajcie uwagi na małą realność tego pytania

Część respondentów wpisała tylko sumę ogólną wydatków pole 4, a część rozbiła na 4a-4d
Jeśli w pyt1 była odpowiedz 2 to 4a(wydatki na organizatora przejazd)=0
ale w innych podkategoriach może coś być (tj. od 4b-4d)

wymyśliłem, że dla każdej kombinacji (pyt 1, pyt2, pyt3) zrobię strukturę wydatków a następnie uzupełnię dane przemnażając wydatki zadeklarowane ogółem przez tą strukturę (o ile były same zera w 4a-4d)

Kod:
tabela dane ma takie kolumny pyt1,pyt2,pyt3, wydatki ogółem, 4a,4b,4c,4d

dane10P=subset(dane, 4a+4b+4c+4d>0, ) #zawężenie do tych choć trochę wypełnionych

Tab_przest1=aggregate(list(dane10P[,5],dane10P[,6],dane10P[,7],dane10P[,8],dane10P[,4]),list(dane10P[,1],dane10P[,2],dane10P[,3]), FUN=sum, na.rm=TRUE) # tabela przestawna sumy wydatków w zależności od odpowiedzi

Tab_przestP=cbind(Tab_przest1,0,0,0,0) # dodanie miejsce na strukturę
Tab_przestP[,9:12]=Tab_przestP[,4:7]/Tab_przestP[,8] #obliczenie struktury wydatków dla każdej z kategorii

for(i in 1:2)# pyt 1
{
  for(j in 1:2)#pyt 2
  {
    for(k in 1:4)#pyt 3
    {baza_idx=which(dane[,1]==i & dane[,2]==j & dane[,3]==k & dane$4a+dane$d4b+dane$4c+dane$4d==0)
   
     dane[baza_idx,5]=dane[baza_idx,4]*Tab_przestP[Tab_przestP[,1]==i & Tab_przestP[,2]==j & Tab_przestP[,3]==k,9]
     baza1[baza_idx,6]=dane[baza_idx,222]*Tab_przestP[Tab_przestP[,1]==i & Tab_przestP[,2]==j & Tab_przestP[,3]==k,10]
     baza1[baza_idx,7]=dane[baza_idx,222]*Tab_przestP[Tab_przestP[,1]==i & Tab_przestP[,2]==j & Tab_przestP[,3]==k,11]
     baza1[baza_idx,8]=dane[baza_idx,222]*Tab_przestP[Tab_przestP[,1]==i & Tab_przestP[,2]==j & Tab_przestP[,3]==k,12]
     }
  }
}


z tym, że coś nie gra ;/
Kod:
Error in `[<-.data.frame`(`*tmp*`, baza_idx, 5, value = numeric(0)) :
  replacement has length zero


Nie wiem co może być nie tak
--------------------------------------------

Inne rozwiązanie wykonuje mi się kilkanaście godzin choć wiem, że działa
Kod:

dane<-read.csv2("testimp.csv")

dane10P=subset(dane, dane[,5]+dane[,6]+dane[,7]+dane[,8]>0, )

Tab_przest1=aggregate(list(dane10P[,5],dane10P[,6],dane10P[,7],dane10P[,8],dane10P[,4]), list(dane10P[,1],dane10P[,2],dane10P[,3]), FUN=sum, na.rm=TRUE)
colnames(Tab_przest1)=c("pyt1","pyt2","pyt3","4a","4b","4c","4d","ogolem")


Tab_przestP=cbind(Tab_przest1,0,0,0,0)
Tab_przestP[,9:12]=Tab_przestP[,4:7]/Tab_przestP[,8]
colnames(Tab_przestP)=c("pyt1","pyt2","pyt3","4a","4b","4c","4d","ogolem","p4a","p4b","p4c","p4d")

daneniewyp=subset(dane, dane[,5]+dane[,6]+dane[,7]+dane[,8]==0 & dane[,4]>0, )

for(i in 1:2)# pyt 1
{
  for(j in 1:2)# pyt2
  {
    for(k in 1:4)#pyt3
    {
      for(m in 1:length(daneniewyp[,1]))
      {if(daneniewyp[m,1]==i& daneniewyp[m,2]==j& daneniewyp[m,3]==k )
        {
        for(g in 1:length(Tab_przestP[,1]))
          {if(Tab_przestP[g,1]==i& Tab_przestP[g,2]==j&Tab_przestP[g,3]==k )
            {daneniewyp[m,5]<-Tab_przestP[g,9]*daneniewyp[m,4]
             daneniewyp[m,6]<-Tab_przestP[g,10]*daneniewyp[m,4]
             daneniewyp[m,7]<-Tab_przestP[g,11]*daneniewyp[m,4]
             daneniewyp[m,8]<-Tab_przestP[g,12]*daneniewyp[m,4]
            }
          }
        }
      }
    }
  }
}

for(o in 1:length(dane[,1]))
{for(p in 1:length(daneniewyp[,1]))
  if(dane[o,5]+dane[o,6]+dane[o,7]+dane[o,8]==0 & dane[o,4]==daneniewyp[p,4]& dane[o,1]==daneniewyp[p,1]& dane[o,2]==daneniewyp[p,2]&dane[o,3]==daneniewyp[p,3])
  {
    dane[o,5]=daneniewyp[p,5]
    dane[o,6]=daneniewyp[p,6]
    dane[o,7]=daneniewyp[p,7]
    dane[o,8]=daneniewyp[p,8]
  }
}
  Temat: Rozkład zmiennej Z=min(x+y)
statnowy

Odpowiedzi: 0
Wyświetleń: 1081

PostForum: Teoria i rachunek prawdopodobieństwa   Wysłany: 2016-09-08, 15:33   Temat: Rozkład zmiennej Z=min(x+y)
Mam dwie zmiene losowe X, Y mają one wartość oczekiwaną 0, oraz różne wariancje

Jak wyznaczyć wartość oczekiwaną zmiennej Z=min(X,Y) oraz jej wariancję.


Jak zacząć? Proszę o dużo wskazówek i opisów
  Temat: Dwa wykresy na jednym
statnowy

Odpowiedzi: 1
Wyświetleń: 1540

PostForum: Biblioteki R, Pakiety R   Wysłany: 2016-08-08, 14:12   Temat: Dwa wykresy na jednym
Kod:

rok    C00M
1977    193
1978    161
1979    166
1980    170
1981    186
1982    164
1983    161
1984    156
1985    158
1986    139
1987    141
1988    122
1989    111
1990    127
1991    81
1992    105
1993    89
1994    89
1995    75
1996    65


Chciałbym stworzyć wykres gdzie wartości dla lat 1977 były by czerwone, a następne (do 1996) niebieskie.

Próbowałem zrobić dwa oddzielne wykresy i je połączyć ale jakoś nie wychodzi

Proszę o wskazówki
  Temat: Kalibracja - jak zrobić
statnowy

Odpowiedzi: 0
Wyświetleń: 1144

PostForum: Biblioteki R, Pakiety R   Wysłany: 2016-07-12, 12:34   Temat: Kalibracja - jak zrobić
Witam

Wczytuję dwie kolumny dane oraz wagi

Dane razy wagi dają mi pewną wartość (suma iloczynów) powiedzmy 10000

Ale teraz muszę zmienić jak najmniej wagi by ta suma iloczynów dawała 11000

Jak najmniej rozumiem jako suma kwadratów zmian ma byś jak najmniejsza.

Znalazłem pakiet Survey ale nie potrafię użyć podanej tu funkcji.

Wiecie jak to można zrobić?

Pozdrawiam
  Temat: zapętlone generowanie wykresów i zapis do wielu plików
statnowy

Odpowiedzi: 5
Wyświetleń: 2587

PostForum: Biblioteki R, Pakiety R   Wysłany: 2016-06-15, 07:55   Temat: zapętlone generowanie wykresów i zapis do wielu plików
Trochę obok tematu więc z góry przepraszam, ale problem dotyczy podobnego problemu - jak zapisywać z pętli wyniki poszczególnych iteracji:

W jaki sposób zapisać wygenerowane w pętli macierze

W jednej pętli napisałem polecenie do generowania (obliczenia) pewnej macierzy i chciałbym ją zapisać (może być w csv)

Problem w tym, że takich macierzy muszę obliczyć i zapisać dużo a na chwilę obecną macierz się nadpisuje i mam tylko jedną - tą ostatnią

Ewentualnie proszę o podanie kodu aby w R po wykonaniu jednej iteracji tworzyła się zmienna mac1 (pierwsza macierz) a po wykonaniu drugiej iteracji zmienna mac2 (druga macierz)
Chyba, że jest inny sposób.

Pozdrawiam
  Temat: Własny pakiet - jak go zrobić, udostepnić i korzystać
statnowy

Odpowiedzi: 1
Wyświetleń: 1426

PostForum: Biblioteki R, Pakiety R   Wysłany: 2016-05-10, 08:40   Temat: Własny pakiet - jak go zrobić, udostepnić i korzystać
Witam

Mam do was pytanie:

Załóżmy, że napisałem pewną funkcję (dla przykładu coś łatwego)

Kod:
sumaliczb<-function(a,b){c=a+b^2; return(c)}



W jaki sposób mogę to zapisać w pakiet i przekazać komuś by on wpisując sumaliczb(4,5) otrzymał wynik 29

proszę o łopatologiczną instrukcję (mam RStudio) lub jakiś link do szczegółowego opisu

Z góry dzięki
  Temat: Wyskakujące "okienka" w R
statnowy

Odpowiedzi: 0
Wyświetleń: 1024

PostForum: Biblioteki R, Pakiety R   Wysłany: 2016-05-09, 14:25   Temat: Wyskakujące "okienka" w R
Witam

Piszę pewien program w R i w pewnym momencie program musi odwołać się do wartości w zmiennej k. Puki co zmienne ta jest wpisana na stałe na początku programu a chciałbym by program w trakcie zapytał się ile ona ma wynosić bo często ją zmieniam co jest uwarunkowane testowaniem różnych ustawień.

Wiem, że taka opcja istnieje w Java i podejrzewam, że w R też da się to zrobić. Zastanawiam się jednak jak

Im prostszy sposób tym lepszy to tylko dla mnie by szybciej testować rożne ustawienia parametrów w programie.
Może jest jakiś pakiet który to umożliwia lub jakaś literatura na ten temat

Pozdrawiam

----

Napisałem coś takiego z zabezpieczeniem
Kod:

readinteger <- function()
{
  n <- readline(prompt="Wpisz całkowitą: ")
  if(!grepl("^[0-9]+$",n))
  {
    return(readinteger())
  }
 
  return(as.integer(n))
}

print(readinteger())


ale ktos mi poradził coś o wiele prostszego
Kod:
x <- as.numeric(readline("Podaj wartość x"))


Nie wiem jednak jak przypisać odpowiedz do zmiennej bo będę chciał z tej wartości skorzystać
  Temat: Osie OX i OY na wykresie
statnowy

Odpowiedzi: 1
Wyświetleń: 1240

PostForum: Biblioteki R, Pakiety R   Wysłany: 2016-04-15, 15:37   Temat: Osie OX i OY na wykresie
Rozwiązanie

Pakiet MASS

Kod:
n <- 1000
X <- mvrnorm(n, mu=c(0,0), Sigma=matrix(c(1,.6,.6,1), ncol=2))
plot(X, xlab="X label", ylab="Y label", pch=19, cex=.4)
abline(h=0, v=0, lwd=1)


Generalnie

abline(h=0, v=0, lwd=1)

załatwia sprawę
  Temat: Wykres Webba
statnowy

Odpowiedzi: 0
Wyświetleń: 1543

PostForum: Biblioteki R, Pakiety R   Wysłany: 2016-04-15, 13:42   Temat: Wykres Webba
Witam,

Mam taki problem a właściwie pytanie: Czy wiecie jak zrobić wykres (diagram) Webba w R

Przykład wykresy Webba

Generalnie mam problem z narysowaniem osi OX i OY w R, ale może jakoś to obejdę

Poważniejszy problem jak połączyć w pary odpowiednie punkty. Wiecie może jak to zrobić? lub znacie odpowiedni pakiet jeżeli takowy istnieje.

Pozdrawiam
  Temat: Osie OX i OY na wykresie
statnowy

Odpowiedzi: 1
Wyświetleń: 1240

PostForum: Biblioteki R, Pakiety R   Wysłany: 2016-04-15, 11:55   Temat: Osie OX i OY na wykresie
Mam takie pytanie : Czy jest możliwość dodania osi OX, OY na wykresie wykonanym w R. Zakładam, że jest jakiś pakiet który to umożliwia, ale może da się inaczej.

np.
plot(sin, -pi, 2*pi) jak dodać osie???

Pozdrawiam
 
Strona 1 z 2
Skocz do:  
Ideą przyświecającą istnieniu forum statystycznego jest stworzenie możliwości wymiany informacji, poglądów i doświadczeń osób związanych ze statystyką, mierzenie się z różnego rodzaju problemami statystycznymi i aktuarialnymi. Poruszane problemy: Statystyka w badaniach sondażowych rynku, metody reprezentacyjne, Teoria i rachunek prawdopodobieństwa, statystyka opisowa, teoria estymacji, testowanie hipotez statystycznych, ekonometria, prognozowanie, metody data mining.
Copyright (C) 2006-2015 Statystycy.pl
Powered by phpBB modified by Przemo © 2003 phpBB Group
Strona wygenerowana w 0,2 sekundy. Zapytań do SQL: 15