[R] Prostej regresja liniowa - lm()

Witam,
Mam do stworzenia prostą regresję dla dwóch ciągów danych znajdujących się w dwóch plikach. W jednym pliku znajdują się prognozy, w drugim rzeczywiste pomiary. W obydwu plikach są to ciągi danych w dwóch kolumnach V1 (data pomiaru) i V2 (wartość pomiaru). Udało mi się wczytać dane z obydwu plików (za pomocą funkcji read.table), tak że po ich wywołaniu dane się wyświetlają, natomiast nie wiem jak z tych danych skorzystać do budowy funkcji regresji. Wiem, że trzeba użyć funkcji lm(), ale kombinuję i po prostu nie wiem jak korzystać z danych z dwóch plików itd....

[R] Symulacje Monte Carlo

Witam,
Chciałbym poczytać o metodach MC w rozwiązywaniu jakiś ciekawych problemów (w tym teoretycznych), możecie polecić jakieś strony bądź książki z porównaniem metod analitycznych i symulacji (przykłady kodu mile widziane)?

[R] Generowanie liczb pseudolosowych - rnorm (), rexp()

Hej,
Czy macie jakieś pomysły jak poprawić/ulepszyć/przyspieszyć poniższy kod, generujący liczby pseudolosowe. Na ogień wziąłem rozkład U(0,1) (czyli w zasadzie po prostu liczby losowe), rozkład wykładniczy i normalny. Zwłaszcza z tym ostatnim mam problemy, jak to napisać wektorowo?

::
#congruential random numbers generator:
Runif<-function(N, x0=unclass(Sys.time()), m=2^31-1, b=7^5)
{
random.seed<-x0
X<-rep(NA,N)
for (j in 1:N)
{
random.seed <- (b * random.seed) %% m
X[j] <- random.seed / m
}
return(X)
}

#check
Runif(N=10,x0=2, m=7,b=3)

#exponential distribution
# inverse transformation method
Rexp<-function(N)
{
U=Runif(N)
X=-log(U) #transforms of uniforms to exp
return(X)
}

#check
Rexp(10)

#random normal (rejection/acceptance method)
Rnorm<-function(N=100)
{
x<-rep(NA,N)
i=0
while (is.na(sum(x))==T)
{[br ...

[R] cor.test()

wyznaczono mase cząsteczkową kwasu bursztynowego. uzyskano następujące wyniki:
Stężenie roztworu: 0,120; 0,238; 0,357; 0,479; 0,558.
Przyrost temp: 0,063; 0,124; 0,186; 0,250; 0,291.
Masa cząsteczkowa: 121,04; 120,04; 117,04; 124,80; 125,32.
czy wyzaczona wartość masy cząsteczkowej zależy od stężenia roztworu. Oblicz średnią masę cząsteczkowa kwasu i porównać ją z wartością teoretyczną m=119,09. α=0,05.

Współczynnik korelacji wyszedł mi 0,5967 , czyli zależność masy cz. od stężenia roztworu jest wysoka(zmienne są skorelowane). Potem wyznaczyłam wartośc t- Studenta= 1,28810152, a t krytyczne= 3,182446305. Więc t<t krytyczne =>nie ma podstaw do odrzucenia hipotezy Ho (zmienne nieskorelowane). No i coś mi tu nie pasuje,ale nie wiem gdzie jest błąd, czy w ogole to jest dobrze policzone.

[R] LINUX

Szukam jakiegoś prostego, lekkiego (czyli nie Emacs :) ) edytora skryptów do R który będzie podświetlał składnię, obsługiwał copy/paste a jeśłi będzie można wysyłać zaznaczenie (jak Tinn-R w Win) do terminala to wogóle będzie pięknie :)

Moje Google-fu dało wyniki w postaci RKward (to jest chyba w zasadzie cały GUI) oraz Vim-a z odpowiednią wtyczką - czy może ktoś ma doświadczenie w pracy z któryms z tych edytorów? Może polecić coś innego?

[R] Wykres funkcja par()

Chciałem wygenerować kilka wykresów obok siebie. Robiłem to milion razy i zawsze działało poprawnie. Nie wiem co jest nie tak z nowym R.

Najpierw ustalam rozmiary

:: par(mfrow=c(2,2)) ::

potem dodaje wykres pierwszy jednak on widoczny jest podwójnie..nie wiem dlaczego..

Dzieki z góry za pomoc.

Pozdrawiam,
B.

[R]Wykorzystanie R w innych językach programowania (C/C++)

Witam,
dopiero poznaję R. Zaciekawił mnie jeden fragment książki p. Biecka, mianowicie:
:: R pozwala na wykonywanie funkcji z bibliotek przygotowywanych w innych
jezykach (C, C++, Fortran) oraz wykonywanie funkcji w R z poziomu innych
jezyków (Java, C++, C i wiele innych) ::

Czy to oznacza, że funkcje napisane w R mogę wykorzystywać przy pisaniu programów w m.in. C czy Fortranie?

[R] seq()

czy instnieje pakiet ktory ulatwi
Mam funkcje

f(a,b)

chce uzyskac tabele ktora zawiera wartosc fukcji w zaleznosci od parametrow a,b czyli cos takiego:

a od 2 od 2.4 krok 0.2

b od 1 do 3 krok 1

1 2 3
2.0 3 3 2
2.2 6 3 2
2.4 3 3 2

wraz z reprezentacja graficzna ;

Ewentualnie czy mozna jakos elegancko stworzyc jakis array i ewentuanie wypelnic ja wartosciami funkcji ?
Z gory dziekuje za odpowiedz