Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

--- cursos:mct:exemplo5 [2011/05/09 10:42] – paulojus
+++ cursos:mct:exemplo5 [2011/05/27 19:00] (atual) – paulojus
@@ Linha 1: / Linha 1: @@
 ====== Miscelânia de tópicos ======
-<fs medium>Um breve exemplo de operações matriciais - coeficientes de regressao(modelo) linear
+<fs medium>Ordenação</fs>
-</fs><code R>
+<code R>
-require(MASS)
+x <- rpois(10, lam=10)
-data(hills)
+x
-names(hills)
+sort(x)
-X <- cbind(1, hills$dist)
+order(x)
-y <- hills$time
+x[order(x)]
+x[rev(order(x))]
-## Há formas e formas de fazer operações...
+args(order)
+order(x, decreasing=T)
-## uma forma não muito adequada:
-solve(t(X) %*% X) %*% t(X) %*% y
-## outra melhor:
-solve(crossprod(X), crossprod(X, y))
 </code>
@@ Linha 73: / Linha 68: @@
 </code>
-<fs medium>Lendo dados de teclado e da área de transferência</fs>
+<fs medium>Precisão numérica/ponto flutuante</fs>
+Encontrando dados replicados, e não replicados
 <code R>
-## digitar dados com <ENTER> após cada registro. Digitar um <ENTER> a mais para encerrar.
+x <- rpois(25, lam=5)
-s1 <- scan()
+x
+unique(x)
+duplicated(x)
+x[!duplicated(x)]
+x[duplicated(x)]
+</code>
+Representação de números e ponto flutuante.\\
+''isTRUE'' deve se usado para testar condições dentro de funções.
+<code R>
+.1 == .3/3
+all.equal(.1, .3/3)
+isTRUE(all.equal(0.1, .3/3))
+identical(.1, .3/3)
+</code>
+Mais coisas que deveriam ser iguais... mas não exatamente são ...
+<code R>
+seq(0, 1, by=0.1) == 0.3
+x <- c(.3, 0.4 - 0.1, 0.5 - 0.2, 0.6 - 0.3, 0.7 - 0.4)
+x
+unique(x)
+</code>
+Um exemplo com uma função para resolver equação do 2o grau:
+<code R>
+eq2g <- function(a, b, c){
+	delta <- b^2 - 4 * a * c
+	if(is.complex(delta) || all(delta > 0))
+	   delta <- sqrt(delta)
+	else
+	   delta <- sqrt(as.complex(delta))
+	return(cbind(-b + c(-1,1) * delta)/(2*a))
+}
+print(eq2g(1, -5, 6), dig=16)
-## recortar (CTRL-C) os seguintes
+eq2g(1, c(-5, 1), 6)
-  13  14
-  21  14
-  31  25
-  32  29
-## lendo como vetor...
+eq2g(1, -5, 6) - c(2, 3)
-s1 <- scan("clipboard")
-## .. e transformando para matriz
-## le por linhas, portanto para montar matriz de ser por linha
-matrix(s1, nc=3, byrow=T)
-## agora lendo diretamente como data-frame
+eq2g(1/3, -5/3, 6/3)
-d1 <- read.table("clipboard")
-d1
-## mais opcoes
+print(eq2g(1/3, -5/3, 6/3), dig=16)
-## recortar (CTRL-C) os seguintes
-X;      Y;   Z;    S
-,3;  13;  14,2;  A
-,2;  21;  14,0;  A
-,1;  31;  25,1;  B
-,1;  32;  29,9;  B
-d2 <- read.table("clipboard", sep=";", dec=",", header=T)
+eq2g(1/3, -5/3, 6/3) - c(2, 3)
-d2
 </code>
-<fs medium>Exemplo lendo tabelas da página do MCT/Indicadores</fs>
-Acessar tabelas na web por www.mct.gov.br/indicadores
+<fs medium>Objetos crescentes, substituição der valores e criação direta de vetores</fs>
+Os tempos de execução podem ser substancialmetne diferentes!!
 <code R>
+n<- 50000
+system.time({
+vec <- numeric(0)
+for(i in 1:n) vec <- c(vec, i)
+})
+system.time({
+vec <- numeric(n)
+for(i in 1:n) vec[i] <- i
+})
+system.time({
+vec <- 1:n
+})
 </code>
+Onderm de operações
+<code R>
+x <- rnorm(10000, m=100, sd=10)
+system.time({
+ lsum <- 0
+ for(i in 1:length(x)){
+   lsum <- lsum + log(x[i])
+  }
+ })
+system.time(lsum <- sum(log(x)))
+c(log(12.3), log(16.7), log(18.3))
+log(c(12.3, 16.7, 18.3))
+</code>

Diferenças

Navegação

Busca

Ferramentas

QR Code