Não foi possível enviar o arquivo. Será algum problema com as permissões?
Essa é uma revisão anterior do documento!
Miscelânia de tópicos
Um breve exemplo de operações matriciais - coeficientes de regressao(modelo) linear
require(MASS) data(hills) names(hills) X <- cbind(1, hills$dist) y <- hills$time ## Há formas e formas de fazer operações... ## uma forma não muito adequada: solve(t(X) %*% X) %*% t(X) %*% y ## outra melhor: solve(crossprod(X), crossprod(X, y))
Representação e Manipulação de datas
## transformando caracteres em datas
dt1 <- as.Date("4/02/2011 03:24:13", "%d/%m/%Y %H:%M:%S")
dt1
## exibindo a data de otras formas
format(dt1, "%j") # dia do ano
format(dt1, "%A") # dia da semana
format(dt1, "%c") # referencia completa (de acordo com o locale do sistema)
format(dt1, "%m") # número do mes
## opcoes disponíveis para exibição:
# a A b B c d H I j m M p S U w W x X y Y z Z
Sys.time()
format(as.Date(Sys.time(), "%d/%m/%Y %H:%M:%S"), "%j")
## as datas podem ser fornecidas de outras maneiras...
x <- c("25set1967", "28set1969", "30jan2003", "05jul2004")
z <- strptime(x, "%d%b%Y")
z
##
## Ver detalhes aqui!!!
##
## e note a confusão com a omisão do século !!!
x <- c("1jan90", "2jan90", "31mar90", "30jul90")
z <- strptime(x, "%d%b%y")
z
x <- c("25set67", "28set69", "30jan03", "05jul04")
z <- strptime(x, "%d%b%y")
z
## concatenando datas e horarios para formar registros únicos
dates <- c("02/27/92", "02/27/92", "01/14/92", "02/28/92", "02/01/92")
times <- c("23:03:20", "22:29:56", "01:03:30", "18:21:03", "16:56:26")
x <- paste(dates, times)
strptime(x, "%m/%d/%y %H:%M:%S")
dt2 <- strptime(x, "%m/%d/%y %H:%M:%S")
format(dt2, "%c")
dt2 <- sort(dt2)
dt2
diff(dt2)
## frações de segundos
z <- strptime("20/2/06 11:16:16.683", "%d/%m/%y %H:%M:%OS")
z
op <- options(digits.secs=3)
z
Lendo dados de teclado e da área de transferência
## digitar dados com <ENTER> após cada registro. Digitar um <ENTER> a mais para encerrar.
s1 <- scan()
## recortar (CTRL-C) os seguintes
12 13 14
21 21 14
12 31 25
15 32 29
## lendo como vetor...
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
d1 <- read.table("clipboard")
d1
## mais opcoes
## recortar (CTRL-C) os seguintes
X; Y; Z; S
12,3; 13; 14,2; A
21,2; 21; 14,0; A
12,1; 31; 25,1; B
15,1; 32; 29,9; B
d2 <- read.table("clipboard", sep=";", dec=",", header=T)
d2
Exemplo lendo tabelas da página do MCT/Indicadores
Acessar tabelas na web por www.mct.gov.br/indicadores