Maio 2016 - R-br - DINF UFPR

Verificar a normalidade em todo o data frame
by Elias Carvalho 15 Mai '16

15 Mai '16

Bom dia Pessoal Eu tenho o seguinte data frame: 'data.frame': 1999 obs. of 14 variables: $ AGE: int 21 31 44 46 49 50 52 64 23 33 ... $ PHA: Factor w/ 2 levels "NO","YES": 1 1 1 1 1 1 1 1 2 2 ... $ SMK: Factor w/ 2 levels "NO","YES": 1 1 1 1 1 1 1 1 1 1 ... $ ALC: Factor w/ 2 levels "NO","YES": 1 1 1 1 1 NA 2 NA 1 NA ... E preciso rodar uma rotina para verificar a normalidade de cada variável de uma só vez, então fiz assim (ps: me perdoem o 'for' eu sou programador das antigas e ainda estou aprendendo a user os ...apply): for (x in 1:length(names(data)))* # vai de 1 a 4 passando por todas as variáveis * { * # Insere o nome da variável em uma matriz* commandP <- paste("nmatrix[",x,",1] <- names(data)[x]", sep="") eval(parse(text=commandP)) * # aplicar apenas para variaveis continuas* * if (!is.factor(data$names(data)[x])) * { *# na coluna 2 da matriz insere valores referente a skewness* commandP <- paste("nmatrix[",x,",2] <- round(skewness(data$", names(data)[x],"),2)", sep="") eval(parse(text=commandP)) *# na coluna3 da matriz insere valores referente a kurtosis* commandP <- paste("nmatrix[",x,",3] <- round(kurtosis(data$", names(data)[x],"),2)", sep="") eval(parse(text=commandP)) * # gera um histograma para cada variável* já com par(mfrow...) calculado de acordo com o # numero de variáveis commandP <- paste("hist(data$", names(data)[x],")", sep="") eval(parse(text=commandP)) *# gera um gráfico de quartil para cada variável já com par(mfrow...) calculado de * * # acordo **com a quantidade de variáveis* commandP <- paste("qqnorm(data$", names(data)[x],")", sep="") eval(parse(text=commandP)) * # Insere a linha no gráfico * commandP <- paste("abline(0,1)", sep="") eval(parse(text=commandP)) } O Resultado que quero parecido com o abaixo, mas para todas as variáveis: nmatrix [,1] [,2] [,3] [1,] "AGE" "-0.13" "2.28" O meu problema é na linha if (!is.factor(data$names(data)[x])) que identifica todas as variáveis como caracter, o que não está errado, pois pega por exemplo "AGE" e então sempre retorna falso e na verdade eu queria algo como "!is.factor(data$AGE), !is.factor(data$SMK).... ou se alguem conhecer um pacote que faça isso de maneira mais fácil agradeço a sugestão Obrigado - Elias <http://lattes.cnpq.br/4248328961021251>

3 4

Retirar preposições, acentos e cedilha dos nomes
by Wagner Tassinari 15 Mai '16

15 Mai '16

Olá pessoal, Como estou trabalhando com dados não estruturados, gostaria de ter uma ideia para fazer uma "limpeza" nesses dados retirando: preposições, acentos, cedilhas e colocar os nomes em caixa alta. Por exemplo, eu tenho o seguinte nome "Maria de Fátima Açudes da Silva", gostaria que ficasse dessa forma: MARIA FATIMA ACUDES SILVA Agradeço a todos ----- Wagner S. Tassinari Departamento de Matemática Universidade Federal Rural do Rio de Janeiro. BR-465, Km 7 - Seropedica, RJ - Brasil CEP: 23890-000 Cel: (21) 96488-5982 (WhatsApp) Skype: wagner.tassinari wtassinari(a)gmail.com ----------------------------------------------------------------------

2 1

Contagem de sequencias baseado no tamanho da cadeia
by Fernando Gama 15 Mai '16

15 Mai '16

Eu tenho uma lista com os seguintes numeros: individuo[[1]]= [8 9 6 0 8 7 8 0 9 2 6 6 0 7 8 2 8 1 3 6 3 2 3 2 0 6 3 3 8 7] Gostaria de detectar (realizar a contagem de sequencias dada um tamanho de cadeia). Baseado nos seguintes filtros: 1) tamanho da cadeia = 2 2) pegar somente os valores iguais ou abaixo de 4 Ai ficaria assim: [8 9 6* 0* 8 7 8 *0* 9 *2* 6 6 *0* 7 8 *2 *8 *1 3* 6 *3 2* *3* *2* *0* 6 *3 3* 8 7] Tentei fazer começar com o which() para pegar as posicoes em sequencia: which(individuo[[1]]<=4) e com rle(individuo[[1]]), este ultimo só funciona para sequencias de numeros iguais. Porém, ainda nao encontrei uma maneira de detectar sequencias baseados no tamanho da cadeia. Como setei lá em cima(tamanho da cadeia = 2), então o output seria essas 4 sequencias: [*1 3* *3 2 * *3 2* *3 3*]. Teria portanto, 8 numeros com 4 sequencias. Alguem pode ajudar? :/ -- Att, Fernando.

3 4

Modelos Mistos - Erro função lme
by Fernando Antonio de souza 14 Mai '16

14 Mai '16

Caros colegas. Estou com dados de um experimento realizado em esquema de parcela subdividida em delineamento inteiramente casualizado. Há dois fatores experimentais A ( 2 níveis) e B (dois níveis) e a resposta Y foi medida repetidamente em tempos diferentes ( 23 tempos), sendo o tempo a subparcela. Eu pretendo analisá-lo através de Modelo Misto, utilizando a função lme{nlme}. Ao analisá-los utilizando a função lme{ nlme}, a função converge, porém recebo mensagens de warnings , avisando que NA's foram produzidos modelo0 <\- lme(GAS~HIDRAT*DILU*TEMP,random=~1|TEMP,weights=varIdent(form=~ 1|TEMP),data=dados) ℝ> anova(modelo0) numDF denDF F-value p-value (Intercept) 1 368 5594.089 <.0001 HIDRAT 1 368 246.749 <.0001 DILU 1 368 2.667 0.1033 TEMP 22 0 206.379 **NaN** HIDRAT:DILU 1 368 6.346 0.0122 HIDRAT:TEMP 22 368 4.199 <.0001 DILU:TEMP 22 368 0.374 0.9961 HIDRAT:DILU:TEMP 22 368 0.038 1.0000 Warning message: In pf(Fval[i], nDF[i], dDF[i]) : NaNs produzidos O mesmo **não** ocorre quando eu utilizo a função lmer{lmerTest} modelo <-lmer(GAS~HIDRAT*DILU*TEMP + (1|TEMP),REML=FALSE,data=dados) ℝ> anova(modelo) Analysis of Variance Table of type III with Satterthwaite approximation for degrees of freedom Sum Sq Mean Sq NumDF DenDF F.value Pr(>F) HIDRAT 2163 2163.24 1 150285 330.58 < 2.2e-16 *** DILU 52 52.32 1 150285 8.00 0.004689 ** TEMP 32798 1490.80 22 150285 227.82 < 2.2e-16 *** HIDRAT:DILU 21 20.66 1 150285 3.16 0.075610 . HIDRAT:TEMP 215 9.77 22 150285 1.49 0.064140 . DILU:TEMP 78 3.53 22 150285 0.54 0.960247 HIDRAT:DILU:TEMP 3 0.13 22 150285 0.02 1.000000 \--- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 A questão é que meus dados são heterocedasticos e gostaria de utilizar a função lme{nlme} devido as facilidades implementadas para lidar com dados assim. Procurei na literatura e estas facilidades ainda não se encontram implementadas na função lmer. Segue CMR dados<-structure(list(HIDRAT = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("1", "12"), class = "factor"), DILU = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L), .Label = c("1", "3"), class = "factor"), TEMP = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L), .Label = c("1", "2", "3", "4", "5", "6", "8", "10"), class = "factor"), BLOC = structure(c(1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L, 1L, 2L, 3L, 4L, 5L), .Label = c("1", "2", "3", "4", "5"), class = "factor"), GAS = c(4.61, 1.91, 0.16, 2.33, 2.64, 0.51, 2.02, 0.48, 1.44, 2.74, 7.24, 5.03, 0.8, 3.61, 7.02, 1.15, 4.7, 1.82, 2.75, 6.8, 8.13, 7.54, 1.49, 3.77, 9.61, 1.48, 7.2, 2.94, 3.59, 9.18, 8.64, 9.84, 4.03, 4.35, 11.67, 1.75, 9.23, 5.81, 3.96, 10.68, 9.27, 13.17, 6.77, 5.72, 13.18, 2.12, 11.49, 7.9, 4.92, 11.84, 9.99, 15.55, 9.31, 7.53, 15.36, 2.42, 13.21, 9.67, 6.44, 13.17, 13.48, 18.07, 14.69, 10.65, 19.52, 4.92, 18.46, 15.12, 9.48, 16.46, 16.96, 23.01, 20.95, 13.55, 24.01, 8.67, 23.56, 19.49, 12.43, 5.7, 5.11, 3.4, 4.73, 4.72, 4.05, 6.88, 4.78, 5.05, 4.5, 9.25, 9.92, 6.04, 7.96, 8.62, 7.12, 11.49, 8.4, 9.63, 8.89, 10.88, 13.5, 8.78, 8.6, 10.38, 8.66, 15.73, 11.27, 10.73, 11.02, 12.32, 16.06, 13.49, 9.42, 11.79, 9.82, 18.54, 16.25, 11.85, 12.44, 14.08, 19.72, 17.72, 10.98, 13.25, 11.11, 21.55, 19.76, 13.42, 13.37, 15.35, 22.71, 20.94, 12.7, 15.21, 12.21, 24.21, 22.36, 15.45, 14.53, 21.57, 26.13, 28.41, 16.8, 19.97, 17.38, 29.57, 29.18, 19.7, 17.78, 28.34, 31.73, 36.12, 22.04, 24.64, 24.89, 35.99, 35.4, 24.09, 21.46)), .Names = c("HIDRAT", "DILU", "TEMP", "BLOC", "GAS"), row.names = c(NA, -159L), class = "data.frame") install.packages(c("nlme","lmerTest")) #instale caso não tenha instalado em seu sistema library(nlme) library (lmerTest) modelo0 <\- lme(GAS~HIDRAT*DILU*TEMP,random=~1|TEMP,weights=varIdent(form=~ 1|TEMP),data=dados) modelo <-lmer(GAS~HIDRAT*DILU*TEMP + (1|TEMP),REML=FALSE,data=dados) ======================================================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal celular: (31)99796-8781 (Vivo) / (31)97358-4685 (Tim) e-mail:nandodesouza@gmail.com Lattes: http://lattes.cnpq.br/6519538815038307 blog: https://producaoanimalcomr.wordpress.com/ ========================================================================

1 0

Uso de if() sentro do while.
by Michelle Bau Graczyk 13 Mai '16

13 Mai '16

Caros, boa tarde, Eu estou tentando montar um programa de troca de calor com um reservatório. Dentro de uma janela de tempo de 63 minutos, se uma dada combinação é maior ou igual a 3, então o reservatório é aberto e recebe uma quantidade de calor. Ele é fechado se essa combinação atinge 0.5, ou se ela atinge um valor que é igual a de abertura mais 2, ou se atinge o limite de tempo, 63 minutos. Como é possível ver abaixo, a operação fecha duas vezes, abre e fecha no mesmo instante( então ele não esta mudando o j para que feche com o valor de abertura +2) e ainda está usando o tempo limite sem precisar. O resultado era para ser uma abertura no instante 17, um encerramento no 44, uma abertura no instante 45 e um encerramento por tempo no instante 63 já que ele não atinge 6.0176641+2 ou 0.5). Eu sei que que o erro está nas posições dos laços mas já mexi de várias formas e nada. Alguém saberia me dizer onde estou errando? Muito obrigada, Michelle vF18 [1] -0.3121193 -0.3121193 -0.3121193 -0.3121193 -2.0091919 -0.3121193 -0.3121193 -0.3121193 [9] -0.3121193 -0.3121193 -0.3121193 -0.3121193 -0.3121193 -0.3121193 -2.0091919 1.0156736 [17] 3.6570085 1.0599739 1.0599739 1.0599739 3.7013087 3.7013087 4.1291396 1.5321050 [25] 2.4320670 1.5321050 1.5321050 2.9041981 2.4763673 3.8041601 3.3763293 3.3763293 [33] 3.8484604 5.5455330 2.9484984 3.8484604 3.8484604 3.8484604 4.3205915 2.6235189 [41] 3.8927606 1.2957260 3.4206295 8.6146987 6.0176641 6.4454950 6.4454950 9.0425296 [49] 6.4454950 6.4454950 4.3205915 4.3205915 4.3205915 3.4649298 3.8927606 3.8927606 [57] 3.8927606 3.8927606 3.8927606 3.8927606 3.8927606 3.8927606 3.8927606 > matrizzF18<-matrix(NA,n,2) > colnames(matrizzF18)<-c("Abre F18","Fecha F18") > j <- 1 > > while (j <= n) { + k<-1 + + if(abs(vF18[j])>=3 && j<63){ + abre_F18<-5*s.F18[j]-round(c18_17)*5*s.F17[j]-round(c18_19)*5*s.F19[j] + matrizzF18[j,1]<-abre_F18 + while (abs(vF18[j + k])>0.5 && j+k<63) { + + if(abs(vF18[j + k])>abs(vF18[j])+2){ + fecha_F18<--5*s.F18[j+k]+round(c18_17)*5*s.F17[j+k]+round(c18_19)*5*s.F19[j+k] + matrizzF18[j + k,2]<-fecha_F18 + + k<-k+1 + break + } + else k <- k + 1 } + + fecha_F18<--5*s.F18[j+k]+round(c18_17)*5*s.F17[j+k]+round(c18_19)*5*s.F19[j+k] + matrizzF18[j + k,2]<-fecha_F18 + } + + j=j+k + } > print(matrizzF18) Abre F18 Fecha F18 [1,] NA NA [2,] NA NA [3,] NA NA [4,] NA NA [5,] NA NA [6,] NA NA [7,] NA NA [8,] NA NA [9,] NA NA [10,] NA NA [11,] NA NA [12,] NA NA [13,] NA NA [14,] NA NA [15,] NA NA [16,] NA NA [17,] 0.4 NA [18,] NA NA [19,] NA NA [20,] NA NA [21,] NA NA [22,] NA NA [23,] NA NA [24,] NA NA [25,] NA NA [26,] NA NA [27,] NA NA [28,] NA NA [29,] NA NA [30,] NA NA [31,] NA NA [32,] NA NA [33,] NA NA [34,] NA NA [35,] NA NA [36,] NA NA [37,] NA NA [38,] NA NA [39,] NA NA [40,] NA NA [41,] NA NA [42,] NA NA [43,] NA NA [44,] NA -0.45 [45,] 0.4 -0.40 [46,] NA NA [47,] NA NA [48,] NA -0.45 [49,] 0.4 -0.40 [50,] NA NA [51,] NA NA [52,] NA NA [53,] NA NA [54,] NA NA [55,] NA NA [56,] NA NA [57,] NA NA [58,] NA NA [59,] NA NA [60,] NA NA [61,] NA NA [62,] NA NA [63,] NA -0.35

1 0

Re: [R-br] Digest R-br, volume 65, assunto 23
by marcos paulo 13 Mai '16

13 Mai '16

Bom dia pessoal, Agradeço as dicas, consegui com a sugestão do Thiago. Com a função: par(mfrow=c(2,2)) não funciona Edson. Obrigado a todos que tiraram um pouco de seu tempo para essa questão. Cordialmente, Marcos Paulo Técnico Agrícola em Zootecnia Bacharel em Agronomia Mestrando em Produção Vegetal - Universidade Federal de GoiásEstagiário de Pós-Graduação da Embrapa Arroz e Feijão/CNPAF-Santo Antônio de Goiás-GO Contato: (62) 85075783 http://lattes.cnpq.br/4322347592884852 > From: r-br-request(a)listas.c3sl.ufpr.br > Subject: Digest R-br, volume 65, assunto 23 > To: r-br(a)listas.c3sl.ufpr.br > Date: Fri, 13 May 2016 12:00:02 -0300 > > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Multiplots usando lattice (marcos paulo) > 2. Re: Multiplots usando lattice (salah) > 3. Re: Multiplots usando lattice (Thiago V. dos Santos) > 4. Re: Colocar legenda sem borda (Éder Comunello) > 5. Re: Multiplots usando lattice (Éder Comunello) > 6. unsubscribe (Nelio Machado) > 7. Re: Multiplots usando lattice (Edson Lira) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 13 May 2016 06:14:26 +0300 > From: marcos paulo <marcospaulo_agronomo(a)hotmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Multiplots usando lattice > Message-ID: <BLU179-W177B236A6AB230A15BD101EE740(a)phx.gbl> > Content-Type: text/plain; charset="iso-8859-1" > > Pessoal estou querendo plotar 4 gráficos na mesma janela gráfica,usando o pacote lattice pra confecção desses plots. Porém não estou conseguindo colocar mais que dois na mesma janela, usando a função print. > Segue exemplo: > g1<- xyplot(NDVI~DAE|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="DAE", ylab="NDVI", main='Vitrine 2014/2016', strip=strip.custom(bg="gray90")) > g2<-xyplot(NDVI~IAF|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="IAF", ylab="NDVI", strip=strip.custom(bg="gray90")) > > g3<-xyplot(Chl~DAE|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="DAE", ylab="Clorofila", strip=strip.custom(bg="gray90")) > g4<-xyplot(Chl~IAF|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="IAF", ylab="Clorofila", strip=strip.custom(bg="gray90")) print(g1, split=c(1,1,1,4), more=TRUE)print(g2, split=c(1,2,1,4))print(g3, split=c(1,3,1,4))print(g4, split=c(1,4,1,4)) > Cordialmente, Marcos Paulo. > Técnico Agrícola em Zootecnia > Bacharel em Agronomia > Mestrando em Produção Vegetal - Universidade Federal de GoiásEstagiário de Pós-Graduação da Embrapa Arroz e Feijão/CNPAF-Santo Antônio de Goiás-GO > Contato: (62) 85075783 > http://lattes.cnpq.br/4322347592884852 > > > Livre de vírus. www.avast.com. > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160513/440c7009/atta…> > > ------------------------------ > > Message: 2 > Date: Fri, 13 May 2016 00:26:20 -0300 > From: salah <salah3.1416(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Multiplots usando lattice > Message-ID: <5735495C.4090807(a)gmail.com> > Content-Type: text/plain; charset="utf-8"; Format="flowed" > > Olá Marcos > > Talvez este exemplo lhe ajude > > library(lattice) > > # Data > w <- as.matrix(dist(Loblolly)) > x <- as.matrix(dist(HairEyeColor)) > y <- as.matrix(dist(rock)) > z <- as.matrix(dist(women)) > > # Plot assignments > pw <- levelplot(w, scales = list(draw = FALSE)) # "scales..." removes axes > px <- levelplot(x, scales = list(draw = FALSE)) > py <- levelplot(y, scales = list(draw = FALSE)) > pz <- levelplot(z, scales = list(draw = FALSE)) > > # Plot prints > print(pw, split = c(1, 1, 2, 2), more = TRUE) > print(px, split = c(2, 1, 2, 2), more = TRUE) > print(py, split = c(1, 2, 2, 2), more = TRUE) > print(pz, split = c(2, 2, 2, 2), more = FALSE) # more = FALSE is redundant > > retirado de > http://stackoverflow.com/questions/2540129/lattice-multiple-plots-in-one-wi… > > saudações > > Em 13/05/2016 00:14, marcos paulo escreveu: > > Pessoal estou querendo plotar 4 gráficos na mesma janela gráfica, > > usando o pacote lattice pra confecção desses plots. Porém não estou > > conseguindo colocar mais que dois na mesma janela, usando a função print. > > > > Segue exemplo: > > > > g1<- xyplot(NDVI~DAE|Cultivar, data=dados, col=1, layout=c(5,1), > > type=c("p","r"), > > xlab="DAE", ylab="NDVI", > > main='Vitrine 2014/2016', > > strip=strip.custom(bg="gray90")) > > > > g2<-xyplot(NDVI~IAF|Cultivar, data=dados, col=1, layout=c(5,1), > > type=c("p","r"), > > xlab="IAF", ylab="NDVI", > > strip=strip.custom(bg="gray90")) > > > > > > g3<-xyplot(Chl~DAE|Cultivar, data=dados, col=1, layout=c(5,1), > > type=c("p","r"), > > xlab="DAE", ylab="Clorofila", > > strip=strip.custom(bg="gray90")) > > > > g4<-xyplot(Chl~IAF|Cultivar, data=dados, col=1, layout=c(5,1), > > type=c("p","r"), > > xlab="IAF", ylab="Clorofila", > > strip=strip.custom(bg="gray90")) > > print(g1, split=c(1,1,1,4), more=TRUE) > > print(g2, split=c(1,2,1,4)) > > print(g3, split=c(1,3,1,4)) > > print(g4, split=c(1,4,1,4)) > > > > Cordialmente, Marcos Paulo. > > > > *Técnico Agrícola em Zootecnia > > Bacharel em Agronomia * > > * > > * > > *Mestrando em Produção Vegetal - Universidade Federal de Goiás* > > *Estagiário de Pós-Graduação da Embrapa Arroz e Feijão/CNPAF-Santo > > Antônio de Goiás-GO* > > * > > * > > *Contato: (62) 85075783* > > > > http://lattes.cnpq.br/4322347592884852 > > > > <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campai…> > > Livre de vírus. www.avast.com > > <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campai…>. > > > > > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel. > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160513/e3b2a537/atta…> > > ------------------------------ > > Message: 3 > Date: Fri, 13 May 2016 11:11:46 +0000 (UTC) > From: "Thiago V. dos Santos" <thi_veloso(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Multiplots usando lattice > Message-ID: > <1996625406.1461259.1463137906486.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset=UTF-8 > > library(gridExtra) > > grid.arrange(g1,g2,g3,g4, ncol=2, nrow=2) > > Greetings, > -- Thiago V. dos Santos > > PhD student > Land and Atmospheric Science > University of Minnesota > > > On Thursday, May 12, 2016 10:26 PM, salah <salah3.1416(a)gmail.com> wrote: > > > > Olá Marcos > > Talvez este exemplo lhe ajude > > library(lattice) > > # Data > w <- as.matrix(dist(Loblolly)) > x <- as.matrix(dist(HairEyeColor)) > y <- as.matrix(dist(rock)) > z <- as.matrix(dist(women)) > > # Plot assignments > pw <- levelplot(w, scales = list(draw = FALSE)) # "scales..." > removes axes > px <- levelplot(x, scales = list(draw = FALSE)) > py <- levelplot(y, scales = list(draw = FALSE)) > pz <- levelplot(z, scales = list(draw = FALSE)) > > # Plot prints > print(pw, split = c(1, 1, 2, 2), more = TRUE) > print(px, split = c(2, 1, 2, 2), more = TRUE) > print(py, split = c(1, 2, 2, 2), more = TRUE) > print(pz, split = c(2, 2, 2, 2), more = FALSE) # more = FALSE is > redundant > > retirado de http://stackoverflow.com/questions/2540129/lattice-multiple-plots-in-one-wi… > > saudações > > > Em 13/05/2016 00:14, marcos paulo escreveu: > > > >Pessoal estou querendo plotar 4 gráficos na mesma janela gráfica, > >usando o pacote lattice pra confecção desses plots. Porém não estou conseguindo colocar mais que dois na mesma janela, usando a função print. > > > > > >Segue exemplo: > > > > > >g1<- xyplot(NDVI~DAE|Cultivar, data=dados, col=1, layout=c(5,1), > > type=c("p","r"), > > xlab="DAE", ylab="NDVI", > > main='Vitrine 2014/2016', > > strip=strip.custom(bg="gray90")) > > > > > >g2<-xyplot(NDVI~IAF|Cultivar, data=dados, col=1, layout=c(5,1), > > type=c("p","r"), > > xlab="IAF", ylab="NDVI", > > strip=strip.custom(bg="gray90")) > > > > > > > > > >g3<-xyplot(Chl~DAE|Cultivar, data=dados, col=1, layout=c(5,1), > > type=c("p","r"), > > xlab="DAE", ylab="Clorofila", > > strip=strip.custom(bg="gray90")) > > > > > >g4<-xyplot(Chl~IAF|Cultivar, data=dados, col=1, layout=c(5,1), > > type=c("p","r"), > > xlab="IAF", ylab="Clorofila", > > strip=strip.custom(bg="gray90")) > > > >print(g1, split=c(1,1,1,4), more=TRUE) > >print(g2, split=c(1,2,1,4)) > >print(g3, split=c(1,3,1,4)) > >print(g4, split=c(1,4,1,4)) > > > > > >Cordialmente, Marcos Paulo. > > > >Técnico Agrícola em Zootecnia > >Bacharel em Agronomia > > > > > >Mestrando em Produção Vegetal - Universidade Federal de Goiás > >Estagiário de Pós-Graduação da Embrapa Arroz e Feijão/CNPAF-Santo Antônio de Goiás-GO > > > > > >Contato: (62) 85075783 > > > > > >http://lattes.cnpq.br/4322347592884852 > > > > > > Livre de vírus. www.avast.com. > > > > > >_______________________________________________ > R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel. > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne� c�igo m�imo reproduz�el. > > > ------------------------------ > > Message: 4 > Date: Fri, 13 May 2016 09:06:14 -0400 > From: Éder Comunello <comunello.eder(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Colocar legenda sem borda > Message-ID: > <CABmC8gnhOuvm2MnK0Uh2Y1TMCsA9++4skj3XhDATz+-Y9Ur4eg(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Bom dia, > > Outra solução... > > B <- matrix(c(578,1289, 439,1061), ncol=2) > colnames(B) <- c("Feminino","Masculino") > rownames(B) <- c("População","IBGE") > barplot(B, beside=T, legend.text=TRUE, col=c(4,"limegreen"), > main="", ylim=c(0,2000), args.legend = c(bty="n")) > > [image: Imagem inline 1] > > > > ================================================ > Éder Comunello > Agronomist (UEM), MSc in Environ. Sciences (UEM) > DSc in Agricultural Systems Engineering (USP/Esalq) > Brazilian Agricultural Research Corporation (Embrapa) > Dourados, MS, Brazil |<O>| > ================================================ > GEO, -22.2752, -54.8182, 408m > UTC-04:00 / DST: UTC-03:00 > > > > > Em 12 de maio de 2016 23:08, Maurício Lordêlo <mslordelo(a)gmail.com> > escreveu: > > > Segue uma sugestão: > > B<-matrix(c(578,1289, 439,1061),ncol=2) > > colnames(B)=c("Feminino","Masculino") > > rownames(B)=c("População","IBGE") > > barplot(B, beside=T,col=c("blue","limegreen"), main="",ylim=c(0,2000)) > > legend(x=5, y=1900,xpd=FALSE, > > legend=c("População","IBGE"),fill=c("blue","limegreen"),bty="n") > > > > Maurício > > > > Em 12 de maio de 2016 17:21, Andre Oliveira <andreolsouza(a)yahoo.com.br> > > escreveu: > > > >> Boa tarde, > >> como posso tirar borda desta legenda? Tentei bty="n" mas não deu certo! > >> B<-matrix(c(578,1289, 439,1061),ncol=2) > >> colnames(B)=c("Feminino","Masculino") > >> rownames(B)=c("População","IBGE") > >> barplot(B, beside=T,legend.text=TRUE,col=c(4,"limegreen"), > >> main="",ylim=c(0,2000)) > >> Obrigado > >> > >> André Oliveira Souza. Graduação em Matemática, mestrado em estatística > >> aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito > >> Santo. IFES > >> > >> _______________________________________________ > >> R-br mailing list > >> R-br(a)listas.c3sl.ufpr.br > >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > >> código mínimo reproduzível. > >> > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > > código mínimo reproduzível. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160513/7b7106fc/atta…> > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: image.png > Tipo: image/png > Tamanho: 9527 bytes > Descrição: não disponível > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160513/7b7106fc/atta…> > > ------------------------------ > > Message: 5 > Date: Fri, 13 May 2016 09:32:32 -0400 > From: Éder Comunello <comunello.eder(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Multiplots usando lattice > Message-ID: > <CABmC8gkJEK_kMbjDzAq=-MiUHK+E_6wfj+tXkGtKpQU-ssB0dw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Marcos, > > Tentou adicionar o argumento more=T nas instruções print de de g2 e g3? > > > ================================================ > Éder Comunello > Agronomist (UEM), MSc in Environ. Sciences (UEM) > DSc in Agricultural Systems Engineering (USP/Esalq) > Brazilian Agricultural Research Corporation (Embrapa) > Dourados, MS, Brazil |<O>| > ================================================ > GEO, -22.2752, -54.8182, 408m > UTC-04:00 / DST: UTC-03:00 > > > > > Em 12 de maio de 2016 23:14, marcos paulo <marcospaulo_agronomo(a)hotmail.com> > escreveu: > > > Pessoal estou querendo plotar 4 gráficos na mesma janela gráfica, > > usando o pacote lattice pra confecção desses plots. Porém não estou > > conseguindo colocar mais que dois na mesma janela, usando a função print. > > > > Segue exemplo: > > > > g1<- xyplot(NDVI~DAE|Cultivar, data=dados, col=1, layout=c(5,1), > > type=c("p","r"), > > xlab="DAE", ylab="NDVI", > > main='Vitrine 2014/2016', > > strip=strip.custom(bg="gray90")) > > > > g2<-xyplot(NDVI~IAF|Cultivar, data=dados, col=1, layout=c(5,1), > > type=c("p","r"), > > xlab="IAF", ylab="NDVI", > > strip=strip.custom(bg="gray90")) > > > > > > g3<-xyplot(Chl~DAE|Cultivar, data=dados, col=1, layout=c(5,1), > > type=c("p","r"), > > xlab="DAE", ylab="Clorofila", > > strip=strip.custom(bg="gray90")) > > > > g4<-xyplot(Chl~IAF|Cultivar, data=dados, col=1, layout=c(5,1), > > type=c("p","r"), > > xlab="IAF", ylab="Clorofila", > > strip=strip.custom(bg="gray90")) > > > > print(g1, split=c(1,1,1,4), more=TRUE) > > print(g2, split=c(1,2,1,4)) > > print(g3, split=c(1,3,1,4)) > > print(g4, split=c(1,4,1,4)) > > > > Cordialmente, Marcos Paulo. > > > > > > *Técnico Agrícola em ZootecniaBacharel em Agronomia * > > > > *Mestrando em Produção Vegetal - Universidade Federal de Goiás* > > *Estagiário de Pós-Graduação da Embrapa Arroz e Feijão/CNPAF-Santo Antônio > > de Goiás-GO* > > > > *Contato: (62) 85075783* > > > > http://lattes.cnpq.br/4322347592884852 > > > > > > <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campai…> Livre > > de vírus. www.avast.com > > <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campai…>. > > <#m_2740539342992165564_DDB4FAA8-2DD7-40BB-A1B8-4E2AA1F9FDF2> > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > > código mínimo reproduzível. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160513/255d26ed/atta…> > > ------------------------------ > > Message: 6 > Date: Fri, 13 May 2016 10:35:16 -0300 > From: Nelio Machado <neliomachado(a)gmail.com> > To: R-br(a)listas.c3sl.ufpr.br > Subject: [R-br] unsubscribe > Message-ID: > <CAKKvw8PyRtPsCy2V-W5=UWihRucwngY-krGCQnCUfMzfthC=Cw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Please, I'd like to unsubscribe this list. > > Thanks > > > Nelio Machado > >  Antes de imprimir pense em seu compromisso com o Meio Ambiente e no > comprometimento com os Custos. > > Se for encaminhar esta mensagem, por favor: > 1. Apague o meu e-mail e o meu nome; > 2. Encaminhe como cópia oculta (Cco ou Bcc) aos seus destinatários; > 3. Apague também, se houver, os endereços dos amigos antes de reenviar; > 4. Agindo sempre assim dificultaremos a disseminação de vírus, spams e > banners. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160513/033f4bc9/atta…> > > ------------------------------ > > Message: 7 > Date: Fri, 13 May 2016 13:47:27 +0000 (UTC) > From: Edson Lira <edinhoestat(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>, > "comunello.eder(a)gmail.com" <comunello.eder(a)gmail.com> > Subject: Re: [R-br] Multiplots usando lattice > Message-ID: > <1766220594.1777376.1463147247246.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Não sei se vai funcionar no lattice.Use: > par(mfrow=c(2,2)) > > > Edson Lira > Estatístico > Manaus-Amazonas > > Em Sexta-feira, 13 de Maio de 2016 9:33, Éder Comunello <comunello.eder(a)gmail.com> escreveu: > > > Marcos, > > Tentou adicionar o argumento more=T nas instruções print de de g2 e g3? > > ================================================ > Éder ComunelloAgronomist (UEM), MSc in Environ. Sciences (UEM) > DSc in Agricultural Systems Engineering (USP/Esalq) > Brazilian Agricultural Research Corporation (Embrapa)Dourados, MS, Brazil |<O>|================================================ > GEO, -22.2752, -54.8182, 408m > UTC-04:00 / DST: UTC-03:00 > > > > Em 12 de maio de 2016 23:14, marcos paulo <marcospaulo_agronomo(a)hotmail.com> escreveu: > > Pessoal estou querendo plotar 4 gráficos na mesma janela gráfica,usando o pacote lattice pra confecção desses plots. Porém não estou conseguindo colocar mais que dois na mesma janela, usando a função print. > Segue exemplo: > g1<- xyplot(NDVI~DAE|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="DAE", ylab="NDVI", main='Vitrine 2014/2016', strip=strip.custom(bg="gray90")) > g2<-xyplot(NDVI~IAF|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="IAF", ylab="NDVI", strip=strip.custom(bg="gray90")) > > g3<-xyplot(Chl~DAE|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="DAE", ylab="Clorofila", strip=strip.custom(bg="gray90")) > g4<-xyplot(Chl~IAF|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="IAF", ylab="Clorofila", strip=strip.custom(bg="gray90")) print(g1, split=c(1,1,1,4), more=TRUE)print(g2, split=c(1,2,1,4))print(g3, split=c(1,3,1,4))print(g4, split=c(1,4,1,4)) > Cordialmente, Marcos Paulo. > Técnico Agrícola em Zootecnia > Bacharel em Agronomia > Mestrando em Produção Vegetal - Universidade Federal de GoiásEstagiário de Pós-Graduação da Embrapa Arroz e Feijão/CNPAF-Santo Antônio de Goiás-GO > Contato: (62) 85075783 > http://lattes.cnpq.br/4322347592884852 > > | | Livre de vírus. www.avast.com. | > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. > > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne� c�igo m�imo reproduz�el. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160513/91ca0564/atta…> > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > > Fim da Digest R-br, volume 65, assunto 23 > *****************************************

1 0

Multiplots usando lattice
by marcos paulo 13 Mai '16

13 Mai '16

Pessoal estou querendo plotar 4 gráficos na mesma janela gráfica,usando o pacote lattice pra confecção desses plots. Porém não estou conseguindo colocar mais que dois na mesma janela, usando a função print. Segue exemplo: g1<- xyplot(NDVI~DAE|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="DAE", ylab="NDVI", main='Vitrine 2014/2016', strip=strip.custom(bg="gray90")) g2<-xyplot(NDVI~IAF|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="IAF", ylab="NDVI", strip=strip.custom(bg="gray90")) g3<-xyplot(Chl~DAE|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="DAE", ylab="Clorofila", strip=strip.custom(bg="gray90")) g4<-xyplot(Chl~IAF|Cultivar, data=dados, col=1, layout=c(5,1), type=c("p","r"), xlab="IAF", ylab="Clorofila", strip=strip.custom(bg="gray90")) print(g1, split=c(1,1,1,4), more=TRUE)print(g2, split=c(1,2,1,4))print(g3, split=c(1,3,1,4))print(g4, split=c(1,4,1,4)) Cordialmente, Marcos Paulo. Técnico Agrícola em Zootecnia Bacharel em Agronomia Mestrando em Produção Vegetal - Universidade Federal de GoiásEstagiário de Pós-Graduação da Embrapa Arroz e Feijão/CNPAF-Santo Antônio de Goiás-GO Contato: (62) 85075783 http://lattes.cnpq.br/4322347592884852 Livre de vírus. www.avast.com.

5 4

unsubscribe
by Nelio Machado 13 Mai '16

13 Mai '16

Please, I'd like to unsubscribe this list. Thanks Nelio Machado  Antes de imprimir pense em seu compromisso com o Meio Ambiente e no comprometimento com os Custos. Se for encaminhar esta mensagem, por favor: 1. Apague o meu e-mail e o meu nome; 2. Encaminhe como cópia oculta (Cco ou Bcc) aos seus destinatários; 3. Apague também, se houver, os endereços dos amigos antes de reenviar; 4. Agindo sempre assim dificultaremos a disseminação de vírus, spams e banners.

1 0

Modelo Misto - Parcela Subdividida em DIC
by Fernando Antonio de souza 12 Mai '16

12 Mai '16

Caros colegas, estou com a seguinte dúvida e gostaria de discutir com vocês. Estou com dados de um experimento realizado em esquema de parcela subdividida em delineamento inteiramente casualizado. Há dois fatores experimentais A ( 2 níveis) e B (dois níveis) e a resposta Y foi medida repetidamente em tempos diferentes ( 23 tempos), sendo o tempo a subparcela. Eu pretendo analisá-lo através de Modelo Misto, utilizando a função lme{nlme}. A minha dúvida é sobre a declaração do modelo aleatório a ser analisado. Neste tipo de experimento a parcela será cada unidade experimental K que recebeu o nível i do tratamento A e nível j do tratamento B e a subparcela o tempo. Estou certo? O modelo aleatório (descreve o erro entre parcelas) então será a interação entre os fatores A e B e repetição (K)? O meu entendimento está correto? O codigo abaixo não é reproduzível, serve apenas para ilustrar o que estou pensando. parcela<\- with(dados, interaction(A,B,K)) modelo<-lme(fixed=Y~A*B*Tempo, random=1|Parcela, data=dados) Atenciosamente ======================================================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal celular: (31)99796-8781 (Vivo) / (31)97358-4685 (Tim) e-mail:nandodesouza@gmail.com Lattes: http://lattes.cnpq.br/6519538815038307 blog: https://producaoanimalcomr.wordpress.com/ ========================================================================

3 3

Acessar microdados RAIS Estabelecimentos
by Carlos Andrade 12 Mai '16

12 Mai '16

Prezados amigos da lista, Solicito ajuda para uma tarefa. Estou tentando acessar os dados do arquivo de microdados da RAIS Estabelecimentos. Tentei com o comando read.fwf, pois o arquivo não possui delimitadores, informando o tamanho dos campos, segundo o arquivo de layout. Sem sucesso. Abri o arquivo txt no editor vim e verifiquei que há 172 colunas, mas quando faço a soma dos campos informados pelo layout só dá 87. Observei também colunas vazias. Repositório dos microdados: ftp://ftp.mtps.gov.br/pdet/microdados/RAIS/2014/ Baixei o arquivo ESTB com UF e CEP.7z Agradeço antecipadamente. -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia

3 4