- R-br - DINF UFPR

out of topic
by gabriel henrique pires de mello ribeiro 21 Mai '16

21 Mai '16

Prezados estou tentando definir o perímetro de copas de árvores após a derrubada. Utilizando um true pulse, medi as distâncias e azimutes ao redor da copa, até o retorno ao ponto de partida. Ainda não consegui achar um pacote adequado, e.g. o pacote geosphere realiza algo similar (mas no caso a partir de coordenadas geográficas ele computa a distância entre linhas). Eu gostaria de fazer o delineamento do perímetro a partir da distância e azimuth entre os pontos (ao redor da copa de cada árvore derrubada). Qualquer sugestão será bem vinda. grato Gabriel Ribeiro

2 3

Wordcloud
by Daniel Spinelli 20 Mai '16

20 Mai '16

Prezados, boa tarde. Preciso fazer uma análise de um texto, contabilizado suas respectivas palavras. Alguma recomendação de pacote ou script para faze - lo? Tenho feito algumas pesquisas do uso da ferramenta Word cloud. Att.: Daniel Spinelli Graduando em Ciências sociais IFCS/UFRJ

3 3

[off-topic] RStudio download lento
by Jefferson Ferreira-Ferreira 20 Mai '16

20 Mai '16

Bom dia pessoal. Me desculpem esse offtopic, mas os servidores do rstudio.org ( https://download1.rstudio.org/rstudio-0.99.902-amd64.deb) parecem estar enfrentando algum problema. O download está terrivelmente lento de modo que não consigo baixar a ultima versão (ubuntu). Será que alguém teria esse .deb que pudesse me disponibilizar? -- *Jefferson Ferreira-Ferreira* Geógrafo – GEOPROCESSAMENTO IDSM | Coordenadoria de TI Jefferson.ferreira(a)mamiraua.org.br *Instituto de Desenvolvimento Sustentável Mamirauá* Ministério da Ciência, Tecnologia e Inovação Telefone: +55 97 3343-9710 *Google Maps* - Mapas deste e-mail: Exibir mapa ampliado <https://maps.google.com.br/maps?q=-3.355557,-64.731151&ll=-3.355471,-64.731…> *Contatos particulares:* *(55) 9615-0100*

1 0

Edição do eixo x
by suzana de souza rodrigues 19 Mai '16

19 Mai '16

Olá,Eu criei um gráfico com esses comandos no R studio:plot (Semana, RESPMe5, xlab = "Year", ylab = "Admissions", main = "Respiratory diseases under 5 years", axes = F) lines ((lowess(Semana, RESPMe5, f = 0.05)),lwd = 2) axis(1, at = seq(0,370, by = 52), pos = 0) axis(2, at = seq(0, 1100, by = 100), pos = 0) Ao invés de aparecer o número de semanas no eixo x, gostaria de escrever o nome dessas semanas: jan-2009, jan-2010, e assim por diante, como eu faço isso?Obrigada!Suzana

1 0

Criar uma coluna de lancamentos de pesca
by Humberto Hazin 19 Mai '16

19 Mai '16

Olá pessoal, Estou precisando criar uma coluna chamada lances. Onde cada lance corresponde a um dia, mês, ano e barco conforme o exemplo abaixo: dia mes ano boat 9 1 2005 AGIOS.NIKOLAUS 10 1 2005 AGIOS.NIKOLAUS 11 1 2005 AGIOS.NIKOLAUS 12 1 2005 AGIOS.NIKOLAUS 13 1 2005 AGIOS.NIKOLAUS Eu gostaria que a tabela ficasse dessa forma aqui dia mes ano boat Lance 9 1 2005 AGIOS.NIKOLAUS 1 10 1 2005 AGIOS.NIKOLAUS 2 11 1 2005 AGIOS.NIKOLAUS 3 12 1 2005 AGIOS.NIKOLAUS 4 13 1 2005 AGIOS.NIKOLAUS 5 Agradeço desde já Humberto

2 12

Script com for() e if {} else{}
by Yury Duarte 19 Mai '16

19 Mai '16

Boa tarde colegas programadores, como vão? Estou tendo alguma dificuldade em acertar o script para executar o calculo do balanço hídrico de Thornthwaite & Mather, mas para uma escala diária. Tenho essas condições que deve ser obedecida para execução dos cálculos: - Se P - ETP < 0: NAc = (P - ETP) + NAc[i-1] ARM = CAD * exp(NAc/CAD) - Se P - ETP >= 0: ARM = (P - ETP) + ARM[i-1] NAc = CAD * Ln(ARM/CAD) - ARM <= CAD Utilizei o comando for e o if else para estabelecer as condições e calcular cada variável para cada dia do período analisado. Entretanto, o script não está calculando as variáveis e retorna os vetores vazios. Segue o script no corpo do email e os dados de entrada em anexo. Agradeço de antemão a ajuda de todos! Abs #----------Remover Objetos----------# rm(list = ls()) #----------Selecionar arquivo de interesse----------# a = read.table("x.txt", sep = "\t", header = T); head(a) #----------Iniciar o BH----------# a$P.ETP = a$RAIN-a$ETP.Penman cad = 100 a$NAc = NA a$ARMi = NA a$ARMf = NA a$ALT = NA a$ETR = NA a$DEF = NA a$EXC = NA head(a) a$NAc[1] = 0 a$ARMi[1] = cad a$ARMf[1] = cad a$ALT[1] = 0 for(i in 2:length(a[,1])){ if(a$P.ETP[i]<0){ (a$NAc[i] = a$P.ETP[i]+a$NAc[i-1]) a$ARMi[i] = (cad*exp(-1*((abs(a$NAc[i]))/cad))) }else{ a$ARMf[i] = ifelse(a$ARMi[i] > cad, cad, a$ARMi[i]) a$ARMi[i] = a$ARMf[i-1]+a$P.ETP[i] a$NAc[i] = cad*log((a$ARMf[i])/cad) } a$ALT[i] = a$ARMf[i] - a$ARMf[i-1] a$ETR[i] = if(a$P.ETP[i]<0){a$RAIN[i]+abs(a$ALT[i])} else{a$ETP.Penman[i]} a$DEF[i] = a$ETP.Penman[i] - a$ETR[i] a$EXC[i] = a$ARMi[i] - a$ARMf[i] } View(a) Yury Duarte Engenheiro Agrônomo - ESALQ/USP

3 3

Separar o segundo nome
by Wagner Tassinari 18 Mai '16

18 Mai '16

Olá pessoal, tudo bem ? Gostaria de um help, estou querendo separar os nomes pelo: primeiro, segundo e último nome de uma pessoa. Para separar o primeiro e o último, eu utilizo os respectivos comandos: banco$firstname = sapply(strsplit(banco$nome, ' '), function(x) x[1]) banco$lastname = sapply(strsplit(banco$nome, ' '), function(x) x[length(x)]) Gostaria de ter uma idéia para separar o segundo nome baseado nesses comandos, Obrigado e abs 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 ----------------------------------------------------------------------

3 3

Modelo mixto - Erro
by Fernando Antonio de souza 16 Mai '16

16 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) ==============================

2 2

Tabela e Graficos no PDF
by Elias Carvalho 15 Mai '16

15 Mai '16

Pessoal Eu gerei uns gráficos e gostaria de imprimir uma tabela na primeira ou última página, estou usando o seguinte comando: par(mfrow=c(2,4)) imprimi gráficos ... par(mfrow=c(1,1)) addtable2plot(-5, -3, nmatrix, bty = "n", display.rownames = FALSE, hlines = TRUE, vlines = TRUE, title = "Normality Results") dev.off() e o resultado é esse abaixo, procurei várias maneiras de imprimir a tabela em pagina separada mas não estou conseguindo, alguém pode me ajudar por favor ? Obrigado -- Best regards... 8^) “The mind that is open to new ideas never come back to its original size” *Albert Einstein* -- Obrigado Elias

3 3

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

Uso de laços e while().
by Michelle Bau Graczyk 12 Mai '16

12 Mai '16

Caros, bom dia, Estou com um problema de lógica de programação. O que eu esperava era o seguinte: se ts[j]>4, eu abro uma operação e ela só para quando ts[j]<4. Enquanto uma operação está aberta, não se pode abrir mais nenhuma. Para isso eu tentei fazer de duas formas, usando e não usando o while, mas eu não consegui. Alguém teria alguma idéia sobre o que eu ando pensando errado? Muito obrigada, Michelle Primeira opção: > v<-c(1,2,5,4,8,9,10,4,2,7,6,8,2) > ts<-ts(v) > > for(j in 1:10){ + if(ts[j]>4){ + #vende + #compra + #compra + vc<-10*ts[j] + while(ts[j+1]>4) j=j+1 + for(k in j+1:10){ + if(ts[k]<4){ + vc2<-7*ts[k] + + } + } + } + else print("Não houve abertura de processo") + j=j+k + vc2 + vc + } [1] "Não houve abertura de processo" [1] "Não houve abertura de processo" [1] "Não houve abertura de processo" Error in if (ts[k] < 4) { : valor ausente onde TRUE/FALSE necessário > vc2 [1] 14 > > vc [1] 80 Segunda opção: > v<-c(1,2,5,4,8,9,10,4,2,7,6,8,2) > ts<-ts(v) > > for(j in 1:10){ + if(ts[j]>4){ + #vende + #compra + #compra + vc<-10*ts[j] + for(k in j+1:10){ + if(ts[k]<4){ + vc2<-7*ts[k] + + } + } + } + else print("Não houve abertura de processo") + j=j+k + vc2 + vc + } [1] "Não houve abertura de processo" [1] "Não houve abertura de processo" [1] "Não houve abertura de processo" Error in if (ts[k] < 4) { : valor ausente onde TRUE/FALSE necessário > vc2 [1] 14 > > vc [1] 80

3 5

Re: [R-br] Criar uma coluna de lancamentos de pesca
by sznelwar＠uol.com.br 12 Mai '16

12 Mai '16

2 1

Faores de risco de doenças - Referência
by Fernando Antonio de souza 11 Mai '16

11 Mai '16

Caros colegas, Estou com um aluno precisando fazer umas análises de dados. A análise que ele necessita é a seguinte: "Estou querendo tentar aprender e aplicar no programa R a estatística para indicar os fatores de risco para doenças, fazendo uma análise univariada, multivariada e por último a obtenção do odds ratio." Não possuo experiência em realizar esta análise. Também não possuo material para referenciar. Alguém pode me indicar alguma referência (livro, apostilas) para que eu possa repassar a este aluno. Att ======================================================================= 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/ ========================================================================

4 4

Manipulação de banco de dados
by Andre Oliveira 10 Mai '16

10 Mai '16

Bom dia,tenho um .txt com algumas linhas e colunas. Em uma coluna tem o fator doença anotados em A ou B; de onde veio a planilha, esqueceram na coluna de doenças células em branco ou ficaram com a minusculo. Como eu poderia programar em R para colocar nestas células vazia a palavra desconhecida por exemplo? E até mesmo trocar os a por A? setwd("/home/andre/Arquivos...") require(gdata) dados<-read.xls(".....xlsx",encoding = "latin1") str(dados) attach(dados)sexo=table(Sexo)sexo a A B 204 4 511 439 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

3 3

Ondulação Geoidal
by Andre Oliveira 10 Mai '16

10 Mai '16

Pessoal bom dia!Tem como pegar ondulação Geoidal usando algum pacote do R? 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

1 0

Números Extraídos do R
by Marcelo Menezes 09 Mai '16

09 Mai '16

Pessoal Boa Noite! Estou desenvolvendo minha monografia, meu co-orientador extraiu para mim dados de variações ambientais quantificadas do WorldClim a partir do R e os números vieram assim, por exemplo; Altitude da Cidade “A” = 82.000.000 m. Sei que na verdade a altitude da Cidade “A” é 820.000 m o ponto aí deveria vir uma casa decimal a frente, mas quando a variável ambiental está relacionada a temperatura, algumas cidades vieram com este valor ex; Cidade “B” 7.900.000 ºC, valor de Isothermality, este valor não pode ser 7,9 e nem muito menos 79 ºC. O que está acontecendo? Não estou entendendo, é algum padrão de número que desconheço? Como resolver isto? Alguém pode me ajudar? Desde já fico muito grato! Att Marcelo Menezes Alves Graduando em Ciências Biológicas pela Universidade Federal de Campina Grande Técnico em Suporte e Manutenção em Informática pelo Instituto Federal da Paraíba. *Lattes: http://lattes.cnpq.br/1992059133849245 <https://wwws.cnpq.br/cvlattesweb/PKG_MENU.menu?f_cod=413CE45ECE38ED23F1637F…>*

4 6

Pacote XTable - Erro de saida
by Elias Carvalho 09 Mai '16

09 Mai '16

Boa Tarde Pessoal Estou testando o pacote "xtable, e o meu ambiente é ubuntu 14.04, RStudio 0.99.896, porém estou tendo o seguinte problema: Quando executo o código abaixo: *# Load librarylibrary(xtable)## Load example datasetdata(tli)## Demonstrate data.frametli.table <- xtable(tli[1:20, ])xtable(mtcars)* O meu resultado não é uma tabela apresentável: % latex table generated in R 3.3.0 by xtable 1.8-2 package % Sun May 8 12:53:18 2016 \begin{table}[ht] \centering \begin{tabular}{rrrrrrrrrrrr} \hline & mpg & cyl & disp & hp & drat & wt & qsec & vs & am & gear & carb \\ \hline Mazda RX4 & 21.00 & 6.00 & 160.00 & 110.00 & 3.90 & 2.62 & 16.46 & 0.00 & 1.00 & 4.00 & 4.00 \\ Já busquei algumas soluções no google mas nada muito explicativo. Alguém pode me ajudar a solucionar o problema ? -- Obrigado Elias

4 5

Script para download automático
by Elias Carvalho 06 Mai '16

06 Mai '16

Boa Tarde Pessoal Alguém conhece algum pacote ou método para fazer download de artigos científicos, explicando: Tenho uma lista de títulos de artigos Gostaria de a partir do título fazer download do texto completo desde que disponíveis Obrigado -- Elias

6 9

Fwd: Script para download automático
by Alisson Lucrécio 06 Mai '16

06 Mai '16

Principio básico, acho que é parametrizar a aba de busca de uma revista ou do google com o título do artigo, e baixar os resultados. O pacote ghost tem funções que poderia te ajudar. Você teria que fazer um loop para os vários artigos. Entretanto, o ghost e python não funcionam muito bem no Windows Se você for usuário do Windows, use uma máquina virtual com o Ubuntu. Segue dois exemplos. Use o ipython notebook (*.ipynb*) (editor de scripts) para abrir, é semelhante ao Rstudio e R. Att., Alisson 2016-05-05 23:36 GMT-03:00 ecacarva [via R-br] < ml-node+s2285057n4666098h74(a)n4.nabble.com>: > Alisson, obrigado pela resposta, por favor me mande os exemplos, vou > tentar por ai. > > Leonardo, obrigado pela dica, conheço e uso também o zotero e o endnote, > mas acho que não fazem exatamente o que quero. > > Em 5 de maio de 2016 19:31, Elias Carvalho <[hidden email] > <http:///user/SendEmail.jtp?type=node&node=4666098&i=0>> escreveu: > >> Ola Leonardo >> >> Conheço o Mendeley, mas até onde sei ele não baixa automaticamente os >> pdfs a partir de uma lista de títulos, ou baixa ? >> >> Se sim você poderia me explicar como faz ? Isso resolveria meu problema. >> >> Obrigado >> >> Em 5 de maio de 2016 14:35, Elias Carvalho <[hidden email] >> <http:///user/SendEmail.jtp?type=node&node=4666098&i=1>> escreveu: >> >>> Boa Tarde Pessoal >>> >>> Alguém conhece algum pacote ou método para fazer download de artigos >>> científicos, explicando: >>> >>> Tenho uma lista de títulos de artigos >>> Gostaria de a partir do título fazer download do texto completo desde >>> que disponíveis >>> >>> Obrigado >>> >>> -- >>> Elias >>> >> >> >> >> -- >> Best regards... 8^) >> >> “The mind that is open to new ideas never come back >> to its original size” *Albert Einstein* >> >> >> _____________________________________________ >> Prof. Elias César Araújo de Carvalho >> CV: http://lattes.cnpq.br/4248328961021251 >> > > > > -- > Best regards... 8^) > > “The mind that is open to new ideas never come back > to its original size” *Albert Einstein* > > > _____________________________________________ > Prof. Elias César Araújo de Carvalho > CV: http://lattes.cnpq.br/4248328961021251 > > _______________________________________________ > R-br mailing list > [hidden email] <http:///user/SendEmail.jtp?type=node&node=4666098&i=2> > 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. > Ecacarva > > > ------------------------------ > If you reply to this email, your message will be added to the discussion > below: > > http://r-br.2285057.n4.nabble.com/R-br-Script-para-download-automatico-tp46… > To unsubscribe from R-br, click here > <http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…> > . > NAML > <http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> > -- Alisson Lucrecio da Costa -- Alisson Lucrecio da Costa

1 0

Usar \mathcal{C} na legenda de gráficos
by Michelle Bau Graczyk 05 Mai '16

05 Mai '16

Caros, bom dia, Estou escrevendo um artigo e preciso usar alguns símbolos nas legendas dos gráficos para ficar igual ao texto. Alguns eu consegui, mas outros não: Não consigo escrever $\widetilde{\mu }(t;1S13)$, tentei expression(widetilde(mu)(t;1S13)) e ão deu certo. Nem quando uso a função paste0 em conjunto. E não consigo escrever letras caligráficas do tipo $\mathcal{C} (s)$. Por favor, alguém teria uma idéia de como posso consegui-los? Muito obrigada, Michelle

2 3

Tabela descritiva
by Elias Carvalho 04 Mai '16

04 Mai '16

Boa Tarde Pessoal Existe algum pacote que auxilie / automatize a criação de uma tabela descritiva, similar a: Obrigado -- Best regards... 8^) “The mind that is open to new ideas never come back to its original size” *Albert Einstein* _____________________________________________ Prof. Elias César Araújo de Carvalho CV: http://lattes.cnpq.br/4248328961021251

3 3

importação de dados em arquivo csv
by Antonio Silva 04 Mai '16

04 Mai '16

Caros, Tenho um arquivo csv que possui dados de longitude e latitude no formato sexagesimal (46°00'00"W / 23°50'00"S) Estou tentando ler o arquivo com dados <- read.csv("dados.csv",sep=";") mas recebo a mensagem Error in make.names(col.names, unique = TRUE) : string multibyte inválida 3 Já verifiquei que se eu tirar os símbolos de grau, minuto e segundo (° ' " ) o arquivo é importado corretamente. Alguém sabem alguma forma de importar os dados sem a necessidade de retirar os símbolos? Obrigado. Antônio Olinto Ex: Lon;Lat;Cod 46°05'00"W;23°50'00"S;585300 46°05'00"W;23°50'00"S;585331 46°05'00"W;23°50'00"S;585514 46°05'00"W;23°50'00"S;585726 46°00'00"W;23°50'00"S;585761 46°05'00"W;23°50'00"S;585884

3 3

Exrpotar formula de modelo para planilha
by Pedro Emmanuel Alvarenga Americano do Brasil 04 Mai '16

04 Mai '16

Amigos de R, Ajustei um modelo de sobrevivencia e gostaria de exportar o modelo uma planilha xlsx ou csv, de tal forma que os coeficientes fiquem numa coluna e na coluna ao lado fique pelo o que o coeficiente está multiplicando. Já quebrei a cabeça aqui e não consegui. A pegadinha é que o modelo possui algumas transformações por splines e por isso preciso dos valores dos nós pra cada variável contínua. Achei uma função chamada Function que faz isso. Mas depois disso eu não avancei. Eu gostaria de pegar da saída função Function o que está entre as chaves, quebrar nos "+" e nos "-" pra separar os termos e depois quebrar nos "*" pra separar as colunas. Organizar em formato de planilha e exportar. library(rms) n <- 1000 age <- 50 + 12*rnorm(n) label(age) <- "Age" sex <- factor(sample(c('Male','Female'), n, rep=TRUE, prob=c(.6, .4))) cens <- 15*runif(n) h <- .02*exp(.04*(age-50)+.8*(sex=='Female')) dt <- -log(runif(n))/h label(dt) <- 'Follow-up Time' e <- ifelse(dt <= cens,1,0) dt <- pmin(dt, cens) units(dt) <- "Year" dd <- datadist(age, sex) options(datadist='dd') S <- Surv(dt,e) f <- cph(S ~ rcs(age,4) + sex, x=TRUE, y=TRUE) >Function(f) >function(age = 48.800654,sex = "Male") {0.28612987-0.017321951* age+0.00012586581*pmax(age-29.600456,0)^3-0.00046285671*pmax(age-44.273198,0)^3+0.00042111763*pmax(age-53.769262,0)^3-8.4126731e-05*pmax(age-69.855597,0)^3-0.64448916*(sex=="Male") } <environment: 0x0000000015fa7378> Pra ficar mais ou menos assim +0.28612987 -0.017321951 age +0.00012586581 pmax(age-29.600456,0)^3 -0.00046285671 pmax(age-44.273198,0)^3 +0.00042111763 pmax(age-53.769262,0)^3 -8.4126731e-05 pmax(age-69.855597,0)^3 -0.64448916 (sex=="Male") Alguma luz? Abraço forte, Pedro Brasil

3 5

Diagostico MLG
by Andre Oliveira 03 Mai '16

03 Mai '16

Boa tarde,alguém sugere algum pacote que pegue objetos do glm e faça direto envelope e outras avaliação do modelo? 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

4 4

Experimentos multi respostas - Otmização
by Andre Oliveira 30 Abr '16

30 Abr '16

Boa noite,alguém conhece algum pacote R para trabalhar com otimização em Experimentos multi respostas? Algo em R sobre técnicade sobreposição da contour individuais ou função de desejabilidade ou até mesmo PO. 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

3 2

Normalização
by Elias Carvalho 29 Abr '16

29 Abr '16

Pessoal To com dados de um banco aqui assim: Group of Age: 1,2,3,4,5 Smoke...........: Yes / No Diabetes........: 83, 79, 130 (Nível de Glicemia) BMI.................: 41.6, 21.6, 34.2 Systilic Pressure: 130, 1140, 160 Preciso rodar um algoritmo para gerar uma rede bayesiana a partir desses dados. Só gostaria de confirmar o seguinte: para obter um resultado mais acurado eu preciso normalizar esses dados certo ? Alguém poderia sugerir um bom pacote para esse tipo de normalização.

3 4

Evapotranspiração Penman-Monteith
by Yury Duarte 29 Abr '16

29 Abr '16

Boa tarde colegas programadores! Estou com algumas dificuldades em rodar o comando "ET.PenmanMonteith" do pacote "Evapotranspiration". Minha intenção é de calcular a evapotranspiração potencial para a escala diária. Segundo o pacote, as variáveis exigidas são: Tmax; Tmin; RHmax; RHmin; Rs e u2. Todas as variáveis exigidas pelo pacote estão compreendidas no meu arquivo de entrada (que segue em anexo), entretanto o comando segue dando erro. Caso alguém seja familiar com o pacote e/ou com o comando e puder me dar uma luz, seria de muita ajuda! Desde já agradeço pela colaboração de todos! Abraços Yury Duarte Engenheiro Agrônomo - ESALQ/USP

3 5

Como fazer gráficos com muitas coordenadas
by Paulo Henrique Pimenta 29 Abr '16

29 Abr '16

Olá Listeiros! Estou com uma problema que tenho certeza que poderão me ajudar! Gostaria de plotar gráficos de muitas curvas. Possuo 6 "x" e 6 "y" e gostaria de fazer como plotar em mesmo gráfico as curvas formadas por cada par de "x" e "y". Um grande abraço a todos! Paulo Henrique de A. S. Pimenta. Graduando em Meteorologia (Bacharelado) - IAG/USP. Fone: +5511981318435.

3 2

Erro em título do eixo y. (ylab/expression)
by Andre Oliveira 28 Abr '16

28 Abr '16

Boa tarde,não consigo plotar a razão (a_n+1)/(a_n) no eixo y, só plota a_n. Alguém poderia ajudar qual descobrir o erro que estou cometendo? plot(t, xlab=expression(a[n]), ylab=expression(frac(a[n+1], a[n])), main=expression(eta==10),col=4, pch=16,ylim=c(0.5,2.5),xlim=c(-0.5,8)) 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

2 1

Como excluir linhas repetidas sem excluir os missings
by Amikobh 28 Abr '16

28 Abr '16

Prezados, Tenho um banco com mais de 8000 linhas e 27 variáveis. Uma das variáveis é a "Processo". Nela, tenho várias linhas em branco e dados repetidos. Como faço para excluir as linhas duplicadas sem excluir as que estão em branco? De tal forma que depois eu tenha um novo banco com estas modificações. Agradecido pela ajuda! PROCESSO 5.155.855-8 5.153.696-8 5.154.220-6 5.154.220-6 5.155.960-6 5.157.157-7 5.154.220-6 Elerson

3 3

Excusão da lista
by Andre Oliveira 27 Abr '16

27 Abr '16

Prezad@s,de alguns dias atrás para cá venho recebendo e-mail dizendo que fui excluído da lista e com link para eu reinscrever na lista. Isto ocorreu 4 vezes, porém, na útima vez não consigo responder a quem responde minhas postagens. Meus tópicos estão na lista respondidos e não os recebo de volta no meu e-mail e sendo assim não há como interagir. Alguém saberia ajudar resolver este problema? Obs: Recebo no e-mail de discussões de outros, mas não recebo as que conversas que eu inicio! 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

3 3

Completar linhas de uma variável
by Amikobh 27 Abr '16

27 Abr '16

Prezados, Tenho um banco com mais de 8000 linhas e 27 variáveis. Uma delas é a "Processo". Nela, tenho várias linhas em branco as quais gostaria de completar com números sequenciais, como exemplo, de 1 até o numero total de linhas em branco (1,2,3,4,.....). De tal forma que depois do preenchimento eu tenha um novo banco com estas modificações. Agradecido pela ajuda! PROCESSO 5.155.855-8 5.153.696-8 5.154.220-6 5.155.401-1 5.155.960-6 5.157.157-7 5.157.658-4 Elerson

3 3

Como recriar um componente de PCA obtido com svyprcomp()
by Leonardo Ferreira Fontenelle 26 Abr '16

26 Abr '16

Boa tarde a todos! Estou tentando escrever os resultados de uma análise de componentes principais de forma que os leitores possam calcular o primeiro componente a partir dos dados sem a necessidade de utilizar um comando de análise de componentes principais propriamente dito. Minha ideia é informar a codificação das variáveis, as cargas, os desvios-padrão, e então orientar os leitores a multiplicar os valores das variáveis pelas cargas e dividir pelos desvios-padrão e então somar tudo. Minha expectativa é de que o resultado seja altamente correlacionado àquele obtido por uma análise de componentes principais propriamente dita. No entanto, como explico a seguir, não estou conseguindo, e de alguma forma isso tem a ver com eu estar usando svyprcomp() em vez de prcomp(). Para começar, mostro uma situação em que isso funciona: pca1 <- prcomp(USArrests, scale = TRUE) table1 <- data.frame(loadings = pca1$rotation[, "PC1"], scale = pca1$scale, coef = NA_real_, row.names = row.names(pca1$rotation)) table1$coef <- table1$loadings / table1$scale firstcomponent1 <- with(USArrests, Murder * table1["Murder", "coef"] + Assault * table1["Assault", "coef"] + UrbanPop * table1["UrbanPop", "coef"] + Rape * table1["Rape", "coef"]) cor(firstcomponent1, pca1$x[, "PC1"]) # resultado: 1,00 Agora, uma situação onde isso não funciona: library("survey") data(api) dclus2 <- svydesign(id = ~ dnum + snum, fpc = ~ fpc1 + fpc2, data = apiclus2) pca2 <- svyprcomp(~ api99 + api00 + ell, design = dclus2, scale = TRUE, scores = TRUE) table2 <- data.frame(loadings = pca2$rotation[, "PC1"], scale = pca2$scale, coef = NA_real_, row.names = row.names(pca2$rotation)) table2$coef <- table2$loadings / table2$scale firstcomponent2 <- with(apiclus2, api99 * table2["api99", "coef"] + api00 * table2["api00", "coef"] + ell * table2["ell", "coef"]) cor(firstcomponent2, pca2$x[, "PC1"]) # resultado: 0,506 Como eu consigo cargas ("loadings") e desvios-padrão ("scale") que eu possa usar de forma a conseguir calcular à mão algum vetor altamente correlacionado com o primeiro componente do svyprcomp()? Grato! Leonardo Ferreira Fontenelle[1] Links: 1. http://lattes.cnpq.br/9234772336296638

2 6

Levene test (original)
by Luiz Roberto Martins Pinto 26 Abr '16

26 Abr '16

Desejo utilizar o Levene test (o original, não o modificado). Alguém conhece alguma função? Não estou conseguindo utilizar a biblioteca 'car', que tem esta função. Luiz Roberto. Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC Laboratório de Estatística Computacional Universidade Estadual de Santa Cruz Ilhéus-Bahia-Brasil luizroberto.uesc(a)gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831

3 5

Estatísticas do R
by Andre Oliveira 26 Abr '16

26 Abr '16

Pessoal,preciso de estatísticas do R tais como nº de bibliotecas oficiais, pacotes que o R já vem na instalação. 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

8 8

Vale a pena tratar esses missings ?
by Elias Carvalho 26 Abr '16

26 Abr '16

Pessoal no meu caso tenho variáveis com alto indice de missing (53 A 79%) : [1] " variable: V1 number.na: 1456 percent.na: 0.73%" [1] " variable: V2 number.na: 1052 percent.na: 0.53%" [1] " variable: V3 number.na: 1577 percent.na: 0.79%" [1] " variable: V4 number.na: 1564 percent.na: 0.78%" Existe alguma regra ? Vale a pena aplicar um método pra tratar esses missing quando a quantidade é alta ?

2 1

Erro no integrate
by Daniela Recchia 26 Abr '16

26 Abr '16

Olá pessoal, quero comparar retas desenhadas por crianças (a mão, elas viam o gráfico com pontos e tinham que desenhar uma reta que melhor se ajustasse aos pontos) com 3 métodos diferentes que resultam na "reta correta". Para isso ploto as retas sempre de duas a duas (criança com método 1, criança com método 2 e criança com método 3) e calculo a área entre as duas retas (integral). Quero no fim poder olhar para os 3 valores de cada criança e achar o mínimo (melhor método). Para isso (tentei) criei um código que infelizmente não funciona... Os valores dos coeficientes das retas dos 3 métodos de referência são fixos (b e s), tenho uma matrix com os 75 valores dos coeficientes para as crianças, os coeficientes são m e t: Quero armazenar cada área clauclada em uma linha e na coluna específica do método de referênciaarea<-function(m,t,na.rm=TRUE) { ‪#‎na‬.rm não funciona u <- matrix(nrow=length(m),ncol=3) colnames(u)<-c("Method1","Method2","Method3") for(i in 1:length(m)){ for (j in 1:3){ f1<-function(x, m, t) {m*x+t} # m e t variáveis locais mt<-function(x) {b[j]*x+s[j]} # b e s variáveis globais dif<-function(x, m,t) {abs(f1(x,m,t)-mt(x))} int[i]<-integrate(dif, -1,1,m = m[i], t = t[i])$value u[i,j]<-int[i] }} return(u)} Recebo o erro de que o objeto "int" não foi encontrado Alguém teria uma sugestão? -- Daniela Recchia “The standard of success in life is absolutely the amount of joy you feel” Esther & Jerry Hick

1 0

Nuvem de Tags
by Elias Carvalho 25 Abr '16

25 Abr '16

Bom dia pessoal Tenho um arquivo formato .csv com 1800 publicações sobre redes bayesianas eu tenho uma coluna identificando o tipo de rede bayesiana, como por exemplo: Baysian Networks (BN), Dynamic Bayesian Networks (DBN), Gaussian Bayesian Networks (GBN), etc... Eu gostaria de ler essa coluna e montar uma nuvem de tags para representar graficamente a frequencia com que esses termos aparecem. Alguem pode dar uma dica ? Obrigado -- Best regards... 8^) “The mind that is open to new ideas never come back to its original size” *Albert Einstein* _____________________________________________ Prof. Elias César Araújo de Carvalho CV: http://lattes.cnpq.br/4248328961021251

3 5

Re: [R-br] Razão de prevalência
by sznelwar 25 Abr '16

25 Abr '16

3 2

Otimização - Fatorial
by Andre Oliveira 23 Abr '16

23 Abr '16

O meu e-mail foi respondido, porém, não sei o motivo de não ter recebido de volta a respostas em meu e-mail! Estou reenviando com a respostas e minha outra colocação. E-mail originalPrezad@s;estou iniciando um projeto em química onde será montado umexperimento no esquema fatorial onde tenho alguns fatores definidos a princípioe gostaria de otimizar este processo. A questão é que tenho 4respostas que serão medidas para cada tratamento. Gostaria desaber se a biblioteca rsm, lm e aov seriam a mais indicada para que eu possa analisar depois de montar o experimento e buscar o deslocamentosnecessários a fim de obter a otimização do processo nospróximos experimentos. Estou buscando material nesta linha e usando o R de preferência. Minha dúvida também se refere ao fato que tenho encontrado materiais com modelos com apenas uma resposta. Sealguém puder ajudar ficarei grato. O Ivan Bezerra Allaman respondeu da seguinte forma O é que eu entendi sua dúvida... vars=list()for(i in 1:30) vars[[i]] <- rnorm(48)names(vars) <- paste('var',1:30,sep='')vars varis <- do.call('cbind',vars) da <- expand.grid(FatA=rep(LETTERS[1:6],4), FatB=letters[1:2]) dad <- data.frame(da,varis) tabs <- apply(dad[,-c(1:2)], 2, function(x) aov(x ~ FatA*FatB,dad))lapply(tabs,summary) Oi Ivan,no caso o que procuro é uma função que faça o que você escreveu mesmo, mas no lugar de x e aov seria algo multivariado como abaixo.fit=manova(varis~da, data= dad) Neste caso a manova seria pouco informativo sobre o otimização do processo. Queria algum modelo com a estrutura do fit que escrevi, mas com mais informações sobre a otimização individual para cada y os fatores em questão. 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

1 0

[off] Ilustração de equações em uma lousa
by Listeiro 037 23 Abr '16

23 Abr '16

Olá a todos. Tenho uma dúvida que não sei onde posso esclarecer-me. Coloquei uma imagem no serviço do imgur para que possam ver. É só clicar no link. http://i.imgur.com/gBIwI1l.jpg A dúvida é se na imagem da lousa acima as equações fazem algum sentido, se no mínimo lembram alguma coisa com algum sentido, matematicamente falando. Se seriam, por exemplo, algo de Estatística ou Análise. Como temos bastante estatísticos aqui na lista creio encontrar a solução. Desde já agradeço.

1 0