- 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