Novembro 2013 - R-br

Re: [R-br] gráfico linhas
by Ivan Bezerra Allaman 18 Nov '13

18 Nov '13

Ivan, obrigado pela atenção. Na sua dica, tentei e não consegui: 1. diminuir a distância entre os gráficos; Veja ?par 2. acrescentar uma legenda indicando os quatro modelos; Veja ?legend 3. acrescentar um título (único) para os três gráficos, mantendo os nomes que separam cada gráfico. Veja ?mtext (s,f,p) Allaman

2 1

[OFF-TOPIC] Shapefile do Brasil com estados divididos em slots
by ASANTOS 18 Nov '13

18 Nov '13

Boa tarde pessoal, Gostaria de saber se alguém sabe onde encontrar um shapelife com o mapa do Brasil com os estados divididos em slots pelo nome, pois gostaria de após ler o shapefile com readShapeLines () e posterior conversão em data frame, fazer seleções do tipo: contorno_line[as.character(contorno_line$ESTADO)=='SaoPaulo',] Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

2 2

Extrair p-valor de um objeto lmer
by Alisson Lucrecio 18 Nov '13

18 Nov '13

Caro Colegas r-br, Bom dia. Eu fiz uma analise com a função lmer "lmer1 <- lmer(Tiametoxam ~ Solo * Dose + Julian + Bloco + (1|Dia), data = dados)" do pacote lme4, mas a anova "anova(lmer1)" não retorna o p-valor. Alguém saberia como obter o p-valor para esse modelo? Obrigado. Alisson Lucrécio da Costa

3 2

Médias diárias em série de 10 anos
by Thiago V. dos Santos 17 Nov '13

17 Nov '13

Problema chato, pessoal: Tenho um data.frame com quatro colunas: variable, date, value e unit. O data.frame foi construído com medidas diárias de 10 anos de cada variável, portanto a data vai de 1 de janeiro de 1980 a aproximadamente 31 de dezembro de 1989 para cada variável. Aqui vai uma amostra do df: > tsvars[1:20,] variable date value unit 2 toa radiation 1980-01-01 258.8793 W/m^2 3 toa radiation 1980-01-02 289.1413 W/m^2 4 toa radiation 1980-01-03 264.9647 W/m^2 5 toa radiation 1980-01-04 291.8841 W/m^2 6 toa radiation 1980-01-05 269.4446 W/m^2 7 toa radiation 1980-01-06 286.0237 W/m^2 8 toa radiation 1980-01-07 274.6440 W/m^2 9 toa radiation 1980-01-08 277.9688 W/m^2 10 toa radiation 1980-01-09 272.6379 W/m^2 11 toa radiation 1980-01-10 273.0654 W/m^2 12 toa radiation 1980-01-11 278.9035 W/m^2 13 toa radiation 1980-01-12 267.1813 W/m^2 14 toa radiation 1980-01-13 270.5033 W/m^2 15 toa radiation 1980-01-14 279.6346 W/m^2 16 toa radiation 1980-01-15 262.5689 W/m^2 17 toa radiation 1980-01-16 267.4670 W/m^2 18 toa radiation 1980-01-17 273.8332 W/m^2 19 toa radiation 1980-01-18 271.3236 W/m^2 20 toa radiation 1980-01-19 271.2731 W/m^2 21 toa radiation 1980-01-20 269.9610 W/m^2 Ao todo eu tenho quatro variáveis e pra cada uma delas há 365 valores (um para cada dia do ano). O dataframe é muito longo e não vou postá-lo completo aqui, mas vejam a saída do str(): > str(tsvars) 'data.frame':18290 obs. of 4 variables: $ variable: Factor w/ 5 levels "latent","sensible",..: 4 4 4 4 4 4 4 4 4 4 ... $ date : Date, format: "1980-01-01" "1980-01-02" "1980-01-03" ... $ value : num 259 289 265 292 269 ... $ unit : chr "W/m^2" "W/m^2" "W/m^2" "W/m^2" ... O que eu preciso fazer é criar a média de cada dia individual dos 10 anos. Ou seja, preciso de uma média para o dia 01/jan, que vai ser a média de todos os 01/jan da série, 02/jan, 03/jan e assim por diante. O detalhe é que a coluna "date" está de fato formatada como data, o que talvez complique as coisas. E eu preciso que ela continue assim, pois facilita na hora de plotar a série temporal. Na sua opinião, qual é a melhor abordagem para tirar a média de cada dia da série, como descrito acima? Saudações, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/ind… Phone: (612) 323 9898

2 2

Carregar MDT no formato *.asc no R
by Helder Gramacho 16 Nov '13

16 Nov '13

Olá Pessoal, Preciso carregar um Modelo Digital do Terreno (raster) no formato *.asc no R gostaria de saber qual pacote é necessário e qual o procedimento. O arquivo quando aberto possui o cabeçalho abaixo: ncols 141 nrows 144 xllcorner 190611.999818 yllcorner 7536478.000059 cellsize 100.000000 nodata_value -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 -9999 Att, Hélder Gramacho *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

3 4

Acentos em nomes da variáveis.
by Andre Oliveira 15 Nov '13

15 Nov '13

Pessoal estou lendo um banco de dados. dados9 <- read.table("arquivos/tab1.txt", h = T) nomes X1 X2 X3 1 Walter 4 0.308 10 2 Jo\xe3o 3 0.235 11 3 S\xe9rgio 5 0.207 4 4 Rui NA 0.270 10 5 Ad\xe3o 7 0.283 5 Como devo ler para sair os acentos corretamente?

3 4

[Ajuda] Knitr - iniciante
by Andre Oliveira 15 Nov '13

15 Nov '13

Pessoal, sou iniciante Knitr já fiz uma busca de material na internet montei uma pasta de material aqui, mas, ainda estou com dúvidas básicas. Aceito indicação de literatura sobre opções de chunks. Estou aqui com um problema com estes códigos. ---------------------------------------------------------------------------------------------------------------------------------------------- <<echo=TRUE, results='asis', fig.width=4, fig.height=3, fig.cap="Histograma teste", echo=T>>= x=rnorm(3000) hist(x,main="",ylab="",xlab="") @ Este gráfico não esta no lugar correto do texto! ----------------------------------------------------------------------------------------------------------------------------------------------- <<????????>>= x=rnorm(66) y=x+rnorm(66) fit=lm(y~x) require(xtable) xtable(anova(fit)) @ Qual comando devo inserir para imprime saida LaTeX? --------------------------------------------------------------------------------------------------------------------------------------------------- Obrigado André Oliveira Souza

5 11

(sem assunto)
by Rodrigo Muniz 15 Nov '13

15 Nov '13

Caros membros da lista, Estou tentando fazer uma função discriminante quadrática para o seguinte conjunto de dados, sendo este um exemplo com apenas 6 linha, mas o conjunto inteiro possui 4000 observações. Os comandos que estou usando segue abaixo bem como o erro encontrado, ressalto que o mesmo dado para a função discriminante linear funciona perfeitamente. Alguém poderia me ajudar com esse erro? Desde já agradeço. LST Altitude MNDWIN NDVI SAVI USO 1 19.36 502 0.50 0.08 0.11 2 2 20.42 611 0.41 0.25 0.38 2 3 17.75 629 0.24 0.39 0.58 1 4 16.73 662 0.22 0.36 0.53 3 5 22.93 596 0.37 0.06 0.10 6 6 19.24 477 -0.04 -0.05 9.41 7 is(LST1985) [1] "data.frame" "list" "oldClass" "vector" qda(LST1985$LST ~ Altitude + MNDWIN + NDVI + SAVI + USO, data=LST1985) Error in qda.default(x, grouping, ...) : some group is too small for 'qda' -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

2 1

Re: [R-br] gráfico linhas
by Ivan Bezerra Allaman 15 Nov '13

15 Nov '13

É isso jovem? x <- seq(1,10) y1 <- rnorm(10) y2 <- rnorm(10) y3 <- rnorm(10) par(mfrow=c(3,1)) plot(y1 ~ x, axes = F, xlab = '', ylab = '', main = 'A', type = 'l') axis(2) plot(y2 ~ x, axes = F, xlab = '', ylab = '', main = 'B', type = 'l') axis(2) plot(y3 ~ x, axes = F, xlab = '', ylab = '', main = 'C', type = 'l') axis(1) axis(2) (s,f,P) Allaman

2 1

Círculo
by FHRB Toledo 15 Nov '13

15 Nov '13

Senhores, Segue um CMR extremamente simples... a pergunta o resultado está certo? require(plotrix) plot(c(-1.2, 1.2), c(-1.2, 1.2), type = 'n', axes = FALSE, xlab = '', ylab = '') rect(-1, -1, 1, 1, lwd = 1.5) draw.circle(0, 0, 1) Por que será que obtenho um gráfico em que o círculo ultrapassa o limite superior do quadrado? att, FH

2 4