disciplinas:ce223:comandos2008

Não foi possível enviar o arquivo. Será algum problema com as permissões?

Você está aqui: start » disciplinas » ce223 » comandos2008

Diferenças

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

--- disciplinas:ce223:comandos2008 [2008/03/04 17:09]
paulojus
+++ disciplinas:ce223:comandos2008 [2008/03/19 15:08]
ehlers
@@ Linha 111: / Linha 111: @@
 </code>
-==== Semana 1 ====
+==== Semana 2 ====
 === 03/03/2008 ===
@@ Linha 162: / Linha 162: @@
 </code>
+Argumentos de funções: ordem dos argumentos, nomes dos argumentos, casamento parcial de nomes, argumentos com //default//
 <code R>
+args(dbinom)
+dbinom(x, 10, 0.3)
+dbinom(x=x, size=10, prob=0.3)
+dbinom(prob=0.3, x=x, size=10)
+dbinom(size=10, prob=0.3, x=x)
+dbinom(x, s=10, p=0.3)
+dbinom(x, p=0.3, s=10)
+dbinom(x, p=0.3, s=10)
+dbinom(x, 10, 0.3, log=F)
+dbinom(x, 10, 0.3, log=T)
 </code>
+=== 05/03/2008 ===
+Verificando máximos, mínimos e suas posições em um vetor.<code R>
+pesos <- c(67, 83, 56, 91, 58, 47, 82, 75)
+max(pesos)
+min(pesos)
+which(pesos == max(pesos))
+which.max(pesos)
+which(pesos == min(pesos))
+which.min(pesos)
+</code>
+Testando pela ocorrência de ''NA''. Note o uso do caracter de negação ''!'' <code R>
+dat <- c(43, 56, NA, 23, 48, 33, NA, 29, 33, 39)
+is.na(dat)
+!is.na(dat)
+which(is.na(dat))
+which(!is.na(dat))
+dat1 <- dat[!is.na(dat)]
+</code>
+==== Semana 3 ====
+=== 10/03/2008 ===
+Outras funções de probabilidades<code R>
+args(dchisq)
+args(dt)
+</code>
+Algumas contantes e operadores<code R>
+exp(1)
+exp(3)
+print(pi, dig=12)
+print(pi, dig=12)
+pi
+options(digits=12)
+pi
+exp(1)
+</code>
+Examinar opções de ''options()''<code R>
+options()
+</code>
+=== 12/03/2008 ===
+Funções para verificar exitência de ''NA'''s ''NaN'''s e valores infinitos (''Inf'') <code R>
+x <- c(5, 0, -2)
+x/0
+is.nan(x/0)
+is.finite(x/0)
+!is.finite(x/0)
+</code>
+Explicações sobre como o R armazana objetos (RAM e/ou dispositivos como por exemplo o HD)<code R>
+save.image()
+q()
+</code>
+Listando e apagando objetos
 <code R>
+ls()
+objects()
+rm(x)
+rm(list=ls())
+args(ls)
+ls()
+ls(all=T)
+</code>
+Tipos de objetos no R: matrizes
+<code R>
+x <- 1:24
+x
+m <- matrix(x, nr=6)
+m
+m <- matrix(x, nr=6, byrow=T)
+m
+x <- 1:25
+m <- matrix(x, nr=6)
+m
+attributes(x)
+attributes(m)
+</code>
+Tipos de objetos no R: arrays
+<code R>
+x <- 1:24
+a <- array(x, dim=c(4,3,2))
+a
+</code>
+Digitação e conversão de uma tabela de tripla entrada (dada no quadro durante a aula) em um objeto do tipo ''array''
+|            |  PR                   ||   SC                 ||   RS                ||
+|            | Masculino | Feminino  | Masculino | Feminino  | Masculino | Feminino |
+|Não Fuma    ^  45       ^  16       ^   21      ^  33       ^  40       ^   45     ^
+|Fuma pouco  ^  28       ^  22       ^   34      ^  21       ^  50       ^   37     ^
+|Fuma muito  ^  37       ^  15       ^   56      ^  30       ^  85       ^   29     ^
+Comentários sobre ordem de entrada dos dados, cliclagem das variáveis e definição das dimensões do array<code R>
+freqs <- scan()
+: 45
+: 28
+: 37
+: 16
+: 22
+: 15
+: 21
+: 34
+: 56
+: 33
+: 21
+: 30
+: 40
+: 50
+: 85
+: 45
+: 37
+: 29
+:
+freqs
+Af <- array(freqs, dim=c(3,2,3))
+Af
+</code>
+==== Semana 4 ====
+=== 17/03/2008 ===
+<code R>
+freqs = scan(file='http://leg.ufpr.br/~ehlers/CE223/fumo.dat')
+array(freqs, dim=c(3,2,3))
+nomes = list(c('PR','SC','RS'), c('M','F'), c('nao fuma','fuma pouco','fuma muito'))
+hf = array(freqs, dim=c(3,2,3), dimnames=nomes)
+hf
+m1 <- matrix(1:12, ncol = 3)
+m1
+dimnames(m1)
+dimnames(m1) <- list(c("L1", "L2", "L3", "L4"), c("C1", "C2", "C3"))
+m1
+m2 <- cbind(1:5, 6:10)
+m2
+m3 <- cbind(1:5, 6)
+m3
+margin.table(m1, margin = 1)
+apply(m1, 1, sum)
+rowSums(m1)
+margin.table(m1, margin = 2)
+apply(m1, 2, sum)
+colSums(m1)
+#operacoes com matrizes
+m4 <- matrix(1:6, nc = 3)
+m5 <- matrix(10 * (1:6), nc = 3)
+m4
+m5
+m4 + m5
+m4 * m5
+m5 - m4
+m5/m4
+m4 %*% m5
+t(m4)
+m6 = t(m4)%*% m5
+solve(m6)
+m6[3,3]=20
+solve(m6)
+mat <- matrix(c(1, 5, 2, 3, -2, 1, -1, 1, -1), nc = 3)
+vec <- c(10, 15, 7)
+solve(mat, vec)
+</code>
+<code R>
+x <- c(0, 1, 2, 3, 4, 5 ,6, 7, 8, 9, 10)
+x
+x <- 0:10
+x
+x <- seq(0,10, by=1)
+x <- scan()
+: 0
+: 1
+...
+: 9
+: 10
+:
+</code>
+Extendendo as possibilidades
+<code R>
+seq(0,1, by=0.1)
+(0:10)/10
+*(0:10)
+seq(0,20,by=2)
+:0
+seq(10,0, by=-1)
+</code>
+=== 19/03/2008 ===
+Data frames
+<code R>
+d1 = data.frame(x=1:10,y=c(51,54,61,67,68,75,77,75,80,82))
+d1
+      x  y
+   1 51
+   2 54
+   3 61
+   4 67
+   5 68
+   6 75
+   7 77
+   8 75
+   9 80
+10 82
+names(d1)
+d1$x
+d1$y
+d1[,1]
+d1[,2]
+plot(d1)
+plot(d1$x,d1$y)
+d2 = data.frame(Y = c(rnorm(5,mean=10,sd=2), rnorm(5,16,2), rnorm(5,14,2)))
+gl(3, 5)
+ [1] 1 1 1 1 1 2 2 2 2 2 3 3 3 3 3
+Levels: 1 2 3
+d2$lev = gl(3, 5)
+d2
+is.factor(d2$lev)
+# Aplicando uma funcao a um Data Frame separando por fatores
+by(d2$Y, d2$lev, summary)
+# Criando um Data Frame a partir de todas as combinacoes de 2 fatores.
+d3 = expand.grid(1:3, 4:5)
+d3
+is.data.frame(d3)
+# Adicionando colunas
+d4 = data.frame(peso=rnorm(15,65,5),altura=rnorm(15,160,10))
+d4
+d4=cbind(d4,sexo=c(rep('M',10),rep('F',5)))
+d4
+# Toda coluna que não seja composta exclusivamente de números é definida como um fator.
+is.factor(d4$sexo)
+[1] TRUE
+letters
+LETTERS
+d4=cbind(d4,nome=letters[1:15])
+is.factor(d4$nome)
+[1] TRUE
+d4$nome= as.character(d4$nome)
+d4
+dim(d4)
+names(d4)
+dimnames(d4)
+rownames(d4)
+colnames(d4)
+d4[d4$sexo=='M',1:2]
+d4[d4$sexo=='F',4]
+by(d4[,1:2],d4$sexo,function(x)x)
+by(d4[,4],d4$sexo,function(x)as.character(x))
+</code>
+Listas
+<code R>
+#Listas sao estruturas genericas e flexiveis que permitem armazenar
+#diversos formatos em um unico objeto.
+lis1 <- list(A = 1:10, B = "CE 223", C = matrix(1:9,ncol = 3))
+lis1
+$A
+[1]  1  2  3  4  5  6  7  8  9 10
+$B
+[1] "CE 223"
+$C
+     [,1] [,2] [,3]
+[1,]    1    4    7
+[2,]    2    5    8
+[3,]    3    6    9
+#Varias funcoes do R retornam listas
+d1
+lis2 = lm (y ~ x, data=d1)
+lis2
+is.list(lis2)
+class(lis2)
+summary(lis2)
+anova(lis2)
+names(lis2)
+lis2$pred
+lis2$residuals
+par(mfrow=c(2,2))
+plot(lis2)
+#Selecionando elementos de uma lista
+lis1$A
+lis2$coeff
+lis1[3]
+lis1[[3]]
 </code>

Diferenças

Navegação

Busca

Ferramentas

QR Code