GoldenLine
Zaloguj się
Justyna Góra

Justyna Góra Student, UG

Temat: Zadanie z R

dostałam do zrobienia zadanie w programie R nie znam sie na tym i bardzo proszę o pomoc. Mam pewien kod do tego zadania jednak czy mógłby mi ktoś sprawdzić czy ten kod jest zgodny z treścią zadania i ewentualnie go poprawić? Bardzo proszę o pomoc, dziękuje. J.

Treść zadania: Sprawdzić symulacyjnie czy przedział ufnoœci dla œ Średniej wartoœci cechy mierzalnej
o rozkładzie normalnym w (1-alfa)100% przypadków obejmuje rzeczywistą œsrednią,
jeżeli próba liczy 20 jednostek.

Kod:

#liczba powtórzeń
nrepl=10000

#długość próbki
n=20
alfa=0.1 #90rzedział ufności
MW=matrix(NA, nrepl,1)

#początek pentli monte-carlo
for (i in 1:nrepl) #petelka od 1 do wartosci nrepl
{ # początek i
y=rnorm(n,175,10) #y jest wektorem posiadajacym losowe liczby o rozkladzie normalnym
sry=mean(y) #srednia z ww.
sdy=sd(y) #standard deviation
talfa=qt(0.05, (n-1)) #rozklad t studenta
talfa=abs(talfa) #to chyba jest wartosc bezwzgledna?
dolna=sry-((talfa*sdy)/sqrt(n)) #dolny przedział ufności
gorna=sry+((talfa*sdy)/sqrt(n)) #górny przedział ufności
indeks1=which(y<:Dolna || y>gorna) #"||" oznacza lub (or)
wiarygodnosc=length(indeks1)/n
MW(i,1)=wiarygodnosc
}# koniec i
mean(MW)
Link do wypowiedziLink
Kamil Bęczyński

Kamil Bęczyński R, SAS, analizy

Temat: Zadanie z R

>talfa=abs(talfa) #to chyba jest wartosc bezwzgledna?
tak, to wartość bezwzględna - wpisz w konsolę "?abs" bez cudzysłowu

ja bym zrobił tak, zamiast :
>indeks1=which(y<:Dolna || y>gorna) #"||" oznacza lub (or)
>wiarygodnosc=length(indeks1)/n
>MW(i,1)=wiarygodnosc

napisał :
MW(i,1)=1-(175<dolna)||(175>gorna)

w ten sposób sprawdzasz czy rzeczywista średnia równa 175 (bo napisałaś y=rnorm(n,175,10)) mieści się wewnątrz wyliczonego przez Ciebie przedziału ufności o postaci : [dolna,gorna]

Jeżeli się nie mieści to wynikiem ciągu działań: 1-(175<dolna)||(175>gorna) jest 0, jeśli się mieści to wynikiem obliczeń jest 1.

Pozdrawiam :)))
Link do wypowiedziLink
Justyna Góra

Justyna Góra Student, UG

Temat: Zadanie z R

Dziękuję za pomoc :)
Link do wypowiedziLink
Konrad J.

Konrad J. Student, UMK

Temat: Zadanie z R

Witam i proszę o pomoc muszę zrobić jedno zadanie w programie R a mianowicie:
Mamy dwie funkcje y=x^2-2, y=x+4. Interesuje nas która funkcja przyjmuje większe wartości w przedziale od -10 do wybranego przez nas „z”, ma to wyglądać tak że mamy skorzystać z pętli for oraz z instrukcji if i utworzyć macierz w pierwszej kolumnie ma być numer wiersza, w drugiej wartości pierwszej funkcji, w trzeciej wartości drugiej funkcji a w ostatniej 4 kolumnie 1 jeżeli wartość funkcji pierwszej jest większa od wartości funkcji drugiej i 0 gdy jest na odwrót. Do tej pozy udało mi zrobić to, ale nie działa to zbyt dobrze, ponieważ ta wewnętrzna pętla for stoi na jednej wartości i nie mam pojęci dlaczego:

funkcja=function (z)
{
A=matrix(0,(z+11),4)
for(i in 1:(z+11))
{for (x in (-10):z)
{y1=((x)^2)-2
y2=(x)+4
A[i,1]=i
A[i,2]=y1
A[i,3]=y2

if (y1>=y2)
{ A[i,4]=1}
else
{ A[i,4]=0}
}
}
A<<-A
}
funkcja(10)
fix(A)
Link do wypowiedziLink
Olivia Mario

Olivia Mario Student, Szkoła
Główna Gospodarstwa
Wiejskiego w
Warszawie

Temat: Zadanie z R

.Ten post został edytowany przez Autora dnia 13.05.14 o godzinie 17:36
Link do wypowiedziLink
Olivia Mario

Olivia Mario Student, Szkoła
Główna Gospodarstwa
Wiejskiego w
Warszawie

Temat: Zadanie z R

Proszę o pomoc, gdyby ktoś umiał to zadanie rozwiązać z R.

Studium przypadku:
-funkcja logistyczna
-funkcja produkcji Cobba - Douglasa

Zbadaj jak kształtują się zmiany średnich cen grunt ów ornych w Polsce (w oparciu o dane GUS). Czy można do
opisu tej zależności wykorzystać regresję logistyczną? W oparciu o regresję logistyczną odpowiedz ile wynosi
wysycenie dla średnich cen gruntów ornych? Czy w oparciu o regresję logistyczną można stwierdzić w jakim
punkcie średnie ceny ziemi zaczęły maleć? Czy regresja logistyczna będzie najlepszym modelem opisującym tą
zależność?
(Sugestia1: Jak wyglądają zmiany cen ziemi na świecie?).
Link do wypowiedziLink
« Wróć do tematów
Odpowiedz

Podobne tematy


Następna dyskusja:

prognozowanie - zadanie




Wyślij zaproszenie do
Anuluj