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

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

Sys.time()
 
format(as.Date("4/02/2011 03:24:13", "%d/%m/%Y %H:%M:%S"), "%j")
format(as.Date(Sys.time(), "%d/%m/%Y %H:%M:%S"), "%j")
# a A b B c d H I j m M p S U w W x X y Y z Z
 
x <- c("1jan1960", "2jan1960", "31mar1960", "30jul1960")
z <- strptime(x, "%d%b%Y")
z
 
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")
 
## fracoes 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 da área de transferência

s1 <- scan()
 
12  13  14  
21  21  14  
12  31  25
15  32  29
 
s1 <- scan("clipboard")
 
d1 <- read.table("clipboard")
d1
 
 
##
## exemplos de leitura e manipulacao para ajustes
##
## site MCT --> indicadores 
 
## Producao cientifica, tabela 5.1 (pesquisadores)
 
pcPesq <- read.table("clipboard", sep="\t", dec=",")
head(pcPesq)
 
pcPesq <- read.table("clipboard", sep="\t", dec=",", strip.white=T, as.is=T)
head(pcPesq)
 
str(pcPesq)
 
ptMil <- function(x){ 
	if(!is.numeric(x) & !is.integer(x)) return(as.numeric(gsub("\\.", "", x)))
        else return(x)
}
pcPesq <- as.data.frame(lapply(pcPesq, ptMil))
 
## Dados socio economicos --> 9.1
ScEc <- read.table("clipboard", dec=",", sep="\t", na.strings="\\u2026 ")
head(ScEc)
 
ScEc <- as.data.frame(lapply(pcPesq, ptMil))


QR Code
QR Code cursos:mct:exemplo5 (generated for current page)