Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
| Próxima revisão Ambos lados da revisão seguinte | |||
|
cursos:mct:exemplo5 [2011/05/06 00:44] paulojus criada |
cursos:mct:exemplo5 [2011/05/09 10:42] paulojus |
||
|---|---|---|---|
| Linha 1: | Linha 1: | ||
| ====== Miscelânia de tópicos ====== | ====== Miscelânia de tópicos ====== | ||
| - | <fs small>Um breve exemplo de operações matriciais - coeficientes de regressao(modelo) linear | + | <fs medium>Um breve exemplo de operações matriciais - coeficientes de regressao(modelo) linear |
| </fs><code R> | </fs><code R> | ||
| require(MASS) | require(MASS) | ||
| Linha 19: | Linha 19: | ||
| - | <fs small>Representação e Manipulação de datas</fs> | + | <fs medium>Representação e Manipulação de datas</fs> |
| <code R> | <code R> | ||
| - | Sys.time() | + | ## 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 | ||
| - | format(as.Date("4/02/2011 03:24:13", "%d/%m/%Y %H:%M:%S"), "%j") | + | ## opcoes disponíveis para exibição: |
| - | 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 | # 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") | + | 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 <- strptime(x, "%d%b%Y") | ||
| z | 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") | 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") | times <- c("23:03:20", "22:29:56", "01:03:30", "18:21:03", "16:56:26") | ||
| Linha 36: | Linha 58: | ||
| strptime(x, "%m/%d/%y %H:%M:%S") | strptime(x, "%m/%d/%y %H:%M:%S") | ||
| - | ## fracoes de segundos | + | 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 <- strptime("20/2/06 11:16:16.683", "%d/%m/%y %H:%M:%OS") | ||
| z | z | ||
| Linha 43: | Linha 73: | ||
| </code> | </code> | ||
| - | <fs small>Lendo dados da área de transferência</fs> | + | <fs medium>Lendo dados de teclado e da área de transferência</fs> |
| <code R> | <code R> | ||
| + | ## digitar dados com <ENTER> após cada registro. Digitar um <ENTER> a mais para encerrar. | ||
| s1 <- scan() | s1 <- scan() | ||
| + | ## recortar (CTRL-C) os seguintes | ||
| 12 13 14 | 12 13 14 | ||
| 21 21 14 | 21 21 14 | ||
| Linha 52: | Linha 84: | ||
| 15 32 29 | 15 32 29 | ||
| + | ## lendo como vetor... | ||
| s1 <- scan("clipboard") | 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 <- read.table("clipboard") | ||
| d1 | 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) |
| - | ## exemplos de leitura e manipulacao para ajustes | + | d2 |
| - | ## | + | </code> |
| - | ## site MCT --> indicadores | + | |
| - | ## Producao cientifica, tabela 5.1 (pesquisadores) | + | <fs medium>Exemplo lendo tabelas da página do MCT/Indicadores</fs> |
| - | pcPesq <- read.table("clipboard", sep="\t", dec=",") | + | Acessar tabelas na web por www.mct.gov.br/indicadores |
| - | head(pcPesq) | + | |
| - | pcPesq <- read.table("clipboard", sep="\t", dec=",", strip.white=T, as.is=T) | + | <code R> |
| - | head(pcPesq) | + | |
| - | str(pcPesq) | ||
| - | ptMil <- function(x){ | + | </code> |
| - | 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)) | ||
| - | |||
| - | </code> | ||
