zajecia2-uzupelnienie

zajecia1/.Rhistory

@@ -0,0 +1,126 @@
+start <- Sys.time()
+for (x in 1:100000000){
+pierw <- sqrt(x)
+}
+print(Sys.time() - start)
+start <- Sys.time()
+for (x in 1:100000000){
+pierw <- x^(1/2)
+}
+print(Sys.time() - start)
+library(schoolmath)
+library("schoolmath")
+install.packages("schoolmath")
+library(schoolmath)
+last <- 0
+for (x in schoolmath.primlist){
+if (x >= 1000){
+print(last)
+break
+}
+last <-x
+}
+install.packages("schoolmath")
+library(schoolmath)
+last <- 0
+for (x in schoolmath.primlist){
+if (x >= 1000){
+print(last)
+break
+}
+last <-x
+}
+install.packages("schoolmath")
+library(schoolmath)
+last <- 0
+for (x in schoolmath.primlist){
+if (x >= 1000){
+print(last)
+break
+}
+last <-x
+}
+library(schoolmath)
+last <- 0
+for (x in primlist){
+if (x >= 1000){
+print(last)
+break
+}
+last <-x
+}
+library(schoolmath)
+last <- 0
+for (x in schoolmath.primlist){
+if (x >= 1000){
+print(last)
+break
+}
+last <-x
+}
+install.packages("schoolmath")
+library(schoolmath)
+last <- 0
+for (x in primlist){
+if (x >= 1000){
+print(last)
+break
+}
+last <-x
+}
+last <- 0
+for (x in data("primlist")){
+if (x >= 1000){
+print(last)
+break
+}
+last <-x
+}
+print(data("primlist"))
+print(data(primlist))
+print(primlist)
+last <- 0
+for (x in primlist){
+if (x >= 1000){
+print(last)
+break
+}
+last <-x
+}
+last <- 0
+number <- 0
+for (x in primlist){
+if (x >= 1000){
+print(last)
+break
+}
+last <-x
+if (x > 100 && x < 500){
+number <- number + 1
+}
+}
+print(number)
+print(rep("a", 3))
+print( paste(rep("a", 3), rep("b", 3)))
+print( paste(rep("a", 3), rep("123", 3)))
+print( paste(rep("a", 3), 1))
+print( paste(rep("a", 3), 1, 2, 3))
+for (x in rep("a", 3)){
+print(paste(x, 1))
+}
+print(paste(rep("a", 3), 1:3))
+print(paste(rep("a", 3), 1:3, sep = ""))
+print(paste(paste(rep("a", 3), 1:3, sep = ""), rep("b", 3), 1:3, sep = ""))
+print(paste(paste(rep("a", 3), 1:3, sep = "")))
+print(paste(paste(rep("b", 3), 1:3, sep = "")))
+print(paste(paste(rep("a", 3), 1:3, sep = "")) + paste(paste(rep("b", 3), 1:3, sep = "")))
+print(c(paste(paste(rep("a", 3), 1:3, sep = "")), paste(paste(rep("b", 3), 1:3, sep = ""))))
+litery <- c("X", "Y", "Z")
+print(paste(1:30, rep(litery, 10)))
+litery <- c("X", "Y", "Z")
+print(paste(1:30, rep(litery, 10), sep="."))
+install.packages("car")
+library(car)
+odpowiedzi <- c("a", "b", "c", "d", "e")
+x <- sample(odpowiedzi, 100, replace = TRUE)
+print(recode(x, "c('a', 'b') = 1; c('c', 'd') = 2; 'e' = 3"))

zajecia2/.Rhistory

@@ -0,0 +1,8 @@
+moja_lista <- list(c("Jakub", "Adamski"), pi)
+moja_lista <- list(c("Jakub", "Adamski"), pi)
+print(moja_lista)
+pierwiastek <- function(x){
+return(sqrt(x))
+}
+moja_lista <- list(c("Jakub", "Adamski"), pi, pierwiastek)
+print(moja_lista)

zajecia2/NY_temp.RData

@@ -0,0 +1,13 @@
+"miesiace" "NY_FARENHEIT" "NY_CELCIUS"
+"1" "Styczeń" 32 0
+"2" "Luty" 33 0.555555555555556
+"3" "Marzec" 41 5
+"4" "Kwiecień" 52 11.1111111111111
+"5" "Maj" 62 16.6666666666667
+"6" "Czerwiec" 72 22.2222222222222
+"7" "Lipiec" 77 25
+"8" "Sierpień" 75 23.8888888888889
+"9" "Wrzesień" 68 20
+"10" "Październik" 58 14.4444444444444
+"11" "Listopad" 47 8.33333333333333
+"12" "Grudzień" 35 1.66666666666667

zajecia2/dane1.csv

@@ -0,0 +1,98 @@
+Wiek;Rozmiar.guza;Wezly.chlonne;Nowotwor;Receptory.estrogenowe;Receptory.progesteronowe;Niepowodzenia;Okres.bez.wznowy;VEGF;
+1;29;1;0;2;(-);(++);brak;22;914
+2;29;1;0;2;(++);(++);brak;53;1118
+3;30;1;1;2;(-);(+);brak;38;630
+4;32;1;0;3;(++);(++);brak;26;1793
+5;32;2;0;NA;(-);(++);brak;19;963
+6;33;1;1;3;(-);(++);wznowa;36;2776
+7;34;1;0;3;(-);(+);brak;33;1676
+8;35;2;1;2;(+);(++);brak;38;3827
+9;35;1;1;3;(-);(-);wznowa;38;22554
+10;36;1;1;2;(-);(++);brak;37;834
+11;36;1;0;2;(++);(++);brak;47;6938
+12;37;1;0;3;(-);(+);wznowa;40;3331
+13;38;2;0;2;(-);(+++);wznowa;30;3505
+14;38;1;1;3;(++);(+++);wznowa;16;2759
+15;38;1;1;3;(-);(-);brak;44;511
+16;40;1;1;2;(+);(++);brak;27;3038
+17;40;1;1;2;(-);(++);brak;36;1027
+18;41;1;1;2;NA;NA;brak;19;1981
+19;41;2;1;2;(++);(+++);brak;41;1489
+20;41;1;0;3;(+);(-);brak;18;989
+21;41;1;1;3;(+);(+);wznowa;41;6175
+22;42;1;1;2;(+++);(++);wznowa;29;5994
+23;42;1;1;2;(++);(++);brak;47;3673
+24;42;1;1;2;(+++);(+++);brak;47;532
+25;42;1;0;2;(+);(-);brak;53;419
+26;42;1;0;NA;(-);(+);brak;33;1531
+27;42;2;1;NA;(-);(+);brak;30;1339
+28;43;1;0;3;NA;NA;brak;39;1348
+29;44;1;0;1;(++);(++);brak;38;962
+30;44;1;0;1;(++);(+);brak;47;2402
+31;44;1;0;2;(-);(+);brak;44;14637
+32;44;2;0;2;(++);(+++);brak;48;483
+33;44;1;0;3;(-);(+);brak;33;1234
+34;44;1;1;3;(+);(++);brak;36;596
+35;44;2;0;3;(-);(-);wznowa;43;3842
+36;44;2;1;3;(-);(-);brak;53;164
+37;44;1;1;NA;(-);(-);brak;38;3836
+38;45;1;0;2;(+++);(++);brak;33;951
+39;45;1;0;2;(++);(++);brak;41;4344
+40;45;1;1;2;(++);(++);brak;54;1275
+41;45;2;0;2;(+++);(+++);brak;51;515
+42;46;1;0;2;(+);(+++);brak;23;2018
+43;46;1;1;2;(-);(-);brak;42;18201
+44;46;1;0;2;(++);(+++);brak;42;1197
+45;46;1;0;2;(-);(+++);brak;45;261
+46;46;1;0;2;(+);(++);brak;51;780
+47;46;2;0;2;(-);(+);brak;42;1524
+48;46;1;1;3;(-);(+);brak;36;2703
+49;46;1;1;3;(+);(++);brak;38;770
+50;46;1;0;NA;NA;NA;brak;28;1526
+51;47;1;0;2;(-);(-);brak;27;2208
+52;47;1;0;2;(+);(++);brak;31;286
+53;47;1;0;2;(+);NA;brak;47;488
+54;47;2;0;2;(+++);(+++);brak;33;2442
+55;47;2;1;2;(-);(-);brak;42;4188
+56;47;1;0;3;(+);(++);brak;38;326
+57;48;1;0;2;(++);(+);brak;34;5785
+58;48;1;1;3;(-);(-);wznowa;21;5194
+59;49;1;0;1;(+);(++);brak;34;4235
+60;49;1;0;1;(++);(++);brak;36;4355
+61;49;1;0;2;(++);(++);brak;26;3240
+62;49;1;1;2;(+);(++);brak;36;3101
+63;49;1;0;2;NA;NA;brak;39;1096
+64;49;2;1;3;(-);(-);brak;39;189
+65;50;1;1;1;(+);(+);brak;49;227
+66;50;2;1;1;(-);(-);brak;28;1485
+67;50;1;1;2;(-);(-);wznowa;19;2221
+68;50;1;0;2;(++);(+++);brak;29;118
+69;50;1;0;2;NA;NA;brak;46;916
+70;50;2;1;2;(++);(-);brak;33;1694
+71;50;1;0;3;(-);(-);brak;45;1652
+72;50;2;0;3;(-);(-);brak;39;1738
+73;50;2;0;3;(-);(-);brak;43;252
+74;50;2;1;3;(+++);(+++);brak;49;3946
+75;50;1;0;NA;(-);(-);brak;41;1945
+76;50;2;1;NA;(+++);(-);brak;27;7665
+77;51;1;0;2;(+);(+);brak;33;4236
+78;51;1;1;2;(++);(++);brak;33;629
+79;51;1;1;2;(+);(+);brak;36;2879
+80;51;1;0;2;(+);(+++);brak;50;223
+81;51;2;1;2;(+++);(++);brak;52;1098
+82;51;2;0;3;(-);(-);wznowa;10;13953
+83;51;1;0;NA;(+);(++);brak;49;175
+84;51;2;1;NA;(-);(-);brak;30;8064
+85;52;1;0;2;(-);(-);brak;30;547
+86;52;1;0;2;(++);(++);brak;42;357
+87;52;2;0;2;(++);(+++);brak;47;2170
+88;52;2;1;2;(+);(+);wznowa;48;1927
+89;52;2;0;3;(+++);(+++);brak;35;456
+90;52;1;0;NA;(-);(+);brak;48;3547
+91;53;1;0;2;(++);(+);brak;43;1067
+92;53;1;0;3;(-);(-);wznowa;50;590
+93;54;2;0;2;(++);(+++);brak;30;624
+94;55;1;0;1;(+);(++);brak;36;1354
+95;55;1;1;2;(++);(++);brak;29;373
+96;55;1;0;NA;NA;NA;brak;NA;1255
+97;57;1;1;2;(++);(++);brak;46;380

zajecia2/zajecia2.R

@@ -0,0 +1,169 @@
+
+# ZAD 1
+pierwiastek <- function(x){
+  return(sqrt(x))
+}
+moja_lista <- list(c("Jakub", "Adamski"), pi, pierwiastek, c(seq(0.02, 1, by=0.02)))
+print(moja_lista)
+
+moja_lista[1] <- NULL
+moja_lista[2] <- NULL
+print(moja_lista)
+
+print(sapply(moja_lista, gamma))
+
+# ZAD 2
+library(Matrix)
+
+m = matrix( rev(c(1,5,3,2,0,5,1,2,1)), nrow = 3, ncol = 3)
+print(rankMatrix(m))
+print(det(m))
+print(solve(m))
+print(eigen(m))
+
+print(sum(m[1,]))
+print(sum(m[2,]))
+print(sum(m[3,]))
+
+print(mean(m[1,]))
+print(mean(m[2,]))
+print(mean(m[3,]))
+
+print(m %*% solve(m))
+
+# ZAD 3
+ostatnia <- function(x){
+  return(x %% 10)
+}
+kwadraty <- c(1:100)^2
+jednosci <- sapply(kwadraty, ostatnia)
+print(table(jednosci))
+
+# ZAD 4
+tabliczka <- outer(paste(c(1:5), "*"), paste(c(1:5), "="), c(1:5)*c(1:5), FUN = "paste")
+print(tabliczka)
+
+# ZAD 5
+dane <- read.table("dane1.csv", header = TRUE, sep = ";")
+print(dane[c(FALSE, TRUE), ])
+filtrowanie1 <- dane[dane$Wiek > 50, ]
+filtrowanie2 <- filtrowanie1[filtrowanie1$Wezly.chlonne == 1, ]
+print(filtrowanie2)
+
+# ZAD 6
+tocelc <- function(x){
+  return(((x-32) * 5)/9)
+}
+miesiace <- c("Styczeń", "Luty", "Marzec", "Kwiecień", "Maj", "Czerwiec", "Lipiec", "Sierpień", "Wrzesień", "Październik", "Listopad", "Grudzień")
+ny_f <- c(32, 33, 41, 52, 62, 72, 77, 75, 68, 58, 47, 35)
+ny_c <- sapply(ny_f, tocelc)
+temperatury <- data.frame(miesiace, ny_f, ny_c)
+names(temperatury)[names(temperatury) == "ny_f"] <- "NY_FARENHEIT"
+names(temperatury)[names(temperatury) == "ny_c"] <- "NY_CELCIUS"
+print(temperatury)
+write.table(temperatury, "NY_temp.RData")
+
+# ZAD 1
+x <- 1:10
+iloczyn <- 1
+for (liczba in x){
+  iloczyn <- iloczyn * liczba
+}
+print(iloczyn)
+
+x <- 1:10
+iloczyn <- 1
+licznik <- 1
+while (licznik <= length(x)){
+  iloczyn <- iloczyn * x[licznik]
+  licznik <- licznik + 1
+}
+print(iloczyn)
+
+x <- 1:10
+iloczyn <- 1
+licznik <- 1
+repeat {
+  iloczyn <- iloczyn * x[licznik]
+  licznik <- licznik + 1
+  if (licznik > length(x)) {
+    break
+  }
+}
+print(iloczyn)
+
+# ZAD 2
+wspolczynnik <- function(n, k){
+  mianownik <- factorial(k) * factorial((n-k))
+  return(factorial(n)/mianownik)
+}
+licznik <- 0
+for (n in 1:100){
+  for (r in 1:n){
+    if (wspolczynnik(n, r) > 1000000){
+      licznik <- licznik + 1
+    }
+  }
+}
+print(licznik)
+
+# ZAD 3
+x <- c(3, 3, 2, 2, 3, 3)
+y <- c(1,2,3,4)
+
+palindrom <- function(x){
+  odwr <- rev(x)
+  licznik <- 1
+  while (licznik <= length(x)){
+    if (odwr[licznik] != x[licznik]){
+      return(FALSE)
+    }
+    licznik <- licznik + 1
+  }
+  return(TRUE)
+}
+
+print(palindrom(x))
+print(palindrom(y))
+
+# ZAD 4
+torad <- function(x){
+  return(x * pi / 180)
+}
+
+katy <- c(0, 30, 45, 60, 90)
+sinn <- sapply(sapply(katy, torad), sin)
+coss <- sapply(sapply(katy, torad), cos)
+tann <- sapply(sapply(katy, torad), tan)
+ctgg <- 1 / sapply(sapply(katy, torad), tan)
+tab <- data.frame(sinn, coss, tann, ctgg)
+print(tab)
+
+
+# ZAD 5
+x <- c(1,2,3,4,5,6,7)
+
+minmax <- function(x){
+  if (length(x) < 3){
+    return("za krótki argument")
+  }
+  liczby <- c()
+  pom <- x
+  for(licznik in 1:3){
+    mini <- min(pom)
+    pom <- c(pom[1:match(mini, pom)-1], pom[(match(mini, pom)+1):length(pom)])
+    liczby <- c(liczby, mini)
+  }
+  pom <- x
+  for(licznik in 1:3){
+    maxi <- max(pom)
+    pom <- c(pom[1:match(maxi, pom)-1], pom[(match(maxi, pom)+1):length(pom)+1])
+    pom <- pom[!is.na(pom)]
+    liczby <- c(liczby, maxi)
+  }
+  return(liczby)
+}
+
+print(minmax(x))
+
+

zajecia2/zajecia2.Rproj

@@ -0,0 +1,13 @@
+Version: 1.0
+
+RestoreWorkspace: Default
+SaveWorkspace: Default
+AlwaysSaveHistory: Default
+
+EnableCodeIndexing: Yes
+UseSpacesForTab: Yes
+NumSpacesForTab: 2
+Encoding: UTF-8
+
+RnwWeave: Sweave
+LaTeX: pdfLaTeX

zajecia3/zajecia3.Rproj

@@ -0,0 +1,13 @@
+Version: 1.0
+
+RestoreWorkspace: Default
+SaveWorkspace: Default
+AlwaysSaveHistory: Default
+
+EnableCodeIndexing: Yes
+UseSpacesForTab: Yes
+NumSpacesForTab: 2
+Encoding: UTF-8
+
+RnwWeave: Sweave
+LaTeX: pdfLaTeX