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) Free statistics help forum. Discuss statistical research, statistical consulting Smarter Poland Portal statystyczny

Poprzedni temat «» Następny temat

Tagi tematu: Brak tagów.

R_zmienna_x_numeryczna?
Autor Wiadomość
Lady Tilly 
Chorąży


Pomogła: 9 razy
Wiek: 37
Posty: 129
Skąd: Poznań
Wysłany: 2010-07-03, 20:30   R_zmienna_x_numeryczna?

Dane są pliki tak jak w zadaniu tym http://www.statystycy.pl/t2836_r-gdzie_blad%3F.php
kod:

> p<-function(dane, proba) {
+ d<-dane[,4]
+ D<-max(d)
+ dane_s<-dane[proba,]
+ xs<-dane_s[,2]
+ ys<-dane_s[,3]
+ ds<-dane_s[,4]
+ d_ciag<-levels(factor(d))
+ dane_r<-dane[(- proba),]
+ dr<-dane_r[,4]
+ sumxrd_ciag<-rowsum((dane_r[,2]),dr, reorder=FALSE)
+
+ nd_ciag<-rep(NA,D)
+ ssxsd_ciag<-rep(NA,D)
+ sumysd_ciag<- rep(NA,D)
+ sumxysd_ciag<- rep(NA,D)
+ for (i in 1:D) {
+ dane_sd<-dane_s[(dane_s[,4]==i),]
+ nd_ciag[i]<-nrow(dane_sd)
+ ssxsd_ciag[i]<-sum((dane_sd[,2])^2)
+ sumysd_ciag[i]<-sum(dane_sd[,3])
+ sumxysd_ciag[i]<-sum(dane_sd[,2]*dane_sd[,3])
+ }
+
+ library(nlme)
+ vc<-VarCorr( lme( ys~ xs - 1, random= ~ xs -1 |ds) )
+ esig2_v<-as.numeric(vc[1,1])
+ esig2_e<-as.numeric(vc[2,1])
+
+ b_ciag<-esig2_e+esig2_v*ssxsd_ciag
+ bodw_ciag<-(1/b_ciag)
+ b_m2_ciag<-b_ciag^(-2)
+ b_m3_ciag<-b_ciag^(-3)
+
+ ebeta1<-(sum(ssxsd_ciag*bodw_ciag))^(-1)
+ ebeta2<-sum(sumxysd_ciag*bodw_ciag)
+ ebeta<-ebeta1*ebeta2
+ pred<-matrix(NA,ncol=D,nrow=1)
+ eMSE<-matrix(NA,ncol=D,nrow=1)
+
+ b1<-sum(((nd_ciag-1)/(esig2_e^2))+ b_m2_ciag)
+ b2<-sum(b_m2_ciag*(ssxsd_ciag^2))
+ b3<-sum(b_m2_ciag*ssxsd_ciag)
+ b<-b1*b2-b3^2
+
+ Ivv<-(2/b)*b2
+ Ive<-(2*b)*b3
+ Iee<-(2/b)*b1
+
+ nawias<-( (esig2_v^2)*Ivv - 2*esig2_e*esig2_v*Ive+(esig2_e^2)*Iee)
+
+ for (i in 1:D) {
+ d_w<-d_ciag[i]
+ Nd_w<-nrow(dane[(dane[,4]==d_w),])
+ nd_w<-nrow(dane_s[(dane_s[,4]==d_w),])
+ Nrd_w<-Nd_w - nd_w
+ sumysd_w<-sumysd_ciag[i]
+ sumxrd_w<-sumixrd_ciag[i]
+ bodw_w<-bodw_ciag[i]
+ sumxysd_w<-sumxysd_ciag[i]
+ ssxsd_w<-ssxsd_ciag[i]
+ b_m3_w<-b_m3_ciag[i]
+
+ pred[i]<-sumysd_w +sumxrd_w *ebeta +bodw_w * esig2_v * sumxrd_w * (sumxysd_w - ebeta* ssxsd_w)
+ eg1<-Nrd_w*esig2_e
+ +esig2_e*esig2_v*bodw_w*(sumxrd_w^2)
+ eg2<-ebeta1*(esig2_e*bodw_w*sumxrd_w)^2
+ eg3<-(sumxrd_w^2)*ssxsd_w* b_m3_w*nawias
+
+ eMSE[i]<-eg1+eg2+2*eg3
+ }
+ wyniki<-matrix(c(pred,eMSE), nrow=2, ncol=D, byrow =TRUE)
+ nazwa<-list(c("predyktor","ocena MSE"),c(1:D))
+ dimnames(wyniki)<-nazwa
+ wyniki
+ }
>
> s<-c(15,16,21,25,30,34,36)
> p(dane,s)
Error in rowsum.default((dane_r[, 2]), dr, reorder = FALSE) :
'x' must be numeric
>

o co chodziz xem?

[ Dodano: 2010-07-04, 11:22 ]
już wiem jak ma być
na początku zabrakło
dane<-as.data.frame(dane)
hahaha człowiek się uczy :mrgreen: na swoich błędach :lol:
_________________
Co to jest statystyka?
- Jeśli trzymasz głowę w lodzie a nogi w ogniu to średnio biorąc jest ci dobrze.
 
 
     
Google

Wysłany:    Reklama google.

 
 
Wyświetl posty z ostatnich:   
Odpowiedz do tematu
Nie możesz pisać nowych tematów
Nie możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach
Nie możesz załączać plików na tym forum
Możesz ściągać załączniki na tym forum
Dodaj temat do Ulubionych zakładek(IE)
Wersja do druku

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,05 sekundy. Zapytań do SQL: 13