Maio 2018 - R-br - DINF UFPR

Modelo linear de segundo grau
by Andre Oliveira 27 Mai '18

27 Mai '18

Sei que posso rodar dois modelos e fazer o LRT, entre outras análises, mas por algum motivo preciso da rsm! Alguém saberia o motivo deste modelo não "rodar"? x=c(1,1,1,2,2,2,3,3,3,4,4,4,5,5,5) y=c(4,8,5,9,12,9,10,13,12,12,14,13,11,11,12) plot(y~x) require(rsm) m3=rsm(y~ FO(x)+FO(x^2)) Error in dimnames(LM$B) <- list(fonm, fonm) : comprimento de 'dimnames' [1] não é igual ao tamanho do array Obg

2 2

Extraindo os raios e área de uma elipse ellipsoidhull() do pacote cluster
by ASANTOS 24 Mai '18

24 Mai '18

Caros Membros, Com fins ilustrativos e exploratório, tenho pontos correspondentes a dois tratamentos A e B espacialmente distribuídos no campo, estou delimitando o limite (zona de influência) de cada um como uma elipse, através da função ellipsoidhull() do pacote cluster conforme CRM abaixo. Gostaria de saber se alguém familiarizado com a função poderia me dizer como faço para extrair os raios e as áreas das elipses EpointsA e EpointsB? #Criação do data frame dados<-NULL #Coordenadas x dados$x<-c(0.350,0.487,0.637,0.775,0.825,0.087,0.237,0.400,0.575,0.737,0.062,0.212) #Coordenadas y dados$y<-c(0.025,0.087,0.050,0.025,0.125,0.187,0.150,0.162,0.212,0.237,0.362,0.337) #Classificação dados$Class<-c("A","A","A","A","A","A","B","B","B","B","B","B") #Visualização plot(dados$x,dados$y) text(dados$x,dados$y+0.01,labels=dados$Class,col="red") # #Usando ellipsoidhull require(cluster) PointsA<-unname(cbind(dados$x[1:6],dados$y[1:6])) EpointsA<-ellipsoidhull(PointsA) lines(predict(EpointsA), col="blue") PointsB<-unname(cbind(dados$x[7:12],dados$y[7:12])) EpointsB<-ellipsoidhull(PointsB) lines(predict(EpointsB), col="orange") # Abraços, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================

3 3

Transformação normal score para krigagem
by Wagner Wolff 22 Mai '18

22 Mai '18

Olá pessoal! Estou trabalhando com dados de chuva numa escala de 15min. Meu banco de dados tem muitos zeros o que deixa a distribuição dos dados assimétrica. Gostaria de interpolar espacialmente esses dados utilizando o geoR, ou seja, modelo geostatistico gaussiano. Assim, a normalização dos dados é necessaria, já tentei a boxcox e não funcionou. Pesquisando sobre o assunto encontrei uma solução na transformação normal score. O problema é como faço a retransformação pra escala original após a krigagem? Alguém tem alguma dica ou outra sugestão para esse problema? Outra característica é que os dados são zero inflacionados. Att.

3 3

TRABALHAR COM O CEP PARA CRIAR LATITUTDE E LONGITUDE EM UM BANCO
by Edson Lira 22 Mai '18

22 Mai '18

Boa tarde Listeiros, estou estudando a elaboração de mapas e me deparei com uma situação que ainda não consegui resolver. Se eu tenho um banco de dados com CEPs, como faço para a partir do CEP criar as variáveis lat lomg. No google mapas tem esse serviço, quando eu tenho um banco pequeno é fácil de resolver, mas quando o banco é grande, só para assinantes. Suponha um banco com 16000 linhas, ai é dificultoso. Eu pesquisei e encontrei um pacote chamado cepR, mas não consegui com o material que está disponível resolver o problema. Alguém poderia sugerir algo? [ ]'s. Prof. Edson Lira, Me Estatístico Manaus-Amazonas

3 3

Gerar pdf de um arquivo .rnw - Sweave
by Andre Oliveira 19 Mai '18

19 Mai '18

bom dia! ao processar meu arquivo .Rnw ocorre este erro. Pesquisei e me parece ser um erro comum, mas não achei a solução. Error in driver$finish(drobj) : the output file 'arquivo.tex' has disappeared Calls: <Anonymous> -> <Anonymous> Alguém tem alguma dica? Obrigado

1 0

Boxplot em grupo
by Marcele de Jesus Correa 09 Mai '18

09 Mai '18

Olá, boa tarde! Sou iniciante no R e preciso de muita ajuda. Tenho recorrido aos tutoriais, livros, algumas dúvidas que são postadas aqui e vídeos no youtube, que tem me ajudado bastante, porém, estou com dificuldade para fazer um gráfico boxplot envolvendo três variáveis juntas. Assim: Estou trabalhando com série temporal de três estações meteorológicas, num total de 27 variáveis, entre estas, estou analisando no momento 3 (precipitação, temperatura mínima e máxima) de cada estação. Para cada uma dessas variáveis e suas respectivas estações, elaborei um gráfico boxplot, conforme comando abaixo: boxplot(TempMin_Jagua ~ Mes, las=1, data = data_estacoes, col = "gray", notch = FALSE, ylab="Temperatura", xlab= "Mês", main= "Média Mensal da Temperatura Mínima da Estação de Jaguaruana"). Entretanto, gostaria de inserir num único gráfico a variável temperatura mínima das três estações, a fim de comparar numa mesma escala o "comportamento" dessa variável por estação, e não estou conseguindo. Li algo parecido aqui no grupo, em que foi utilizado o pacote ggplot2, tentei usar sem sucesso. Fiz-me entender? Espero que sim, porque esse negócio é bem complicado... mas estou gostando. Desde já agradeço. Cordialmente, Marcele Correa. Geographer Master Student in Climate Sciences Research Line - Bioclimatology and C limatic C hanges Advisor - Kellen Carla Lima Federal University of Rio Grande do Norte - UFRN / Brazil.

6 8

Instalação de pacote do github
by Andre Oliveira 06 Mai '18

06 Mai '18

Bom dia! faz algum tempo que este problema vem ocorrendo e dificultando as coisas para mim. Veja bem: require(devtools) Carregando pacotes exigidos: devtools install_github("tagteam/Publish") Installation failed: Problem with the SSL CA cert (path? access rights?) Pesquisando passei usar o comando! httr::set_config(httr::config(ssl_verifypeer = FALSE )) install_github("tagteam/Publish") Downloading GitHub repo tagteam/Publish@master from URL https://api.github.com/repos/tagteam/Publish/zipball/master Installation failed: Peer reports incompatible or unsupported protocol version. Não sei mais o que fazer, pois, preciso atualizar alguns pacotes e não consigo. Alguém teria uma sugestão? obg

2 1

Escalonar matrizes
by Andre Oliveira 06 Mai '18

06 Mai '18

boa tarde, alguém conhece uma função que escalona uma matriz? obg André Oliveira Souza

2 1

Emprego da PCA para seleção de variáveis mais correlacionadas no primeiro componente principal
by ASANTOS 05 Mai '18

05 Mai '18

Prezados Membros, Gostaria de selecionar as 5 variáveis resposta no primeiro componente principal mais correlacionadas usando a PCA, através da função prcomp(), conforme o post https://stats.stackexchange.com/questions/115032/how-to-find-which-variable…, sendo o meu CRM: ##Banco de dados set.seed(12345) mat <- matrix(rnorm(120,0,0.5),nrow=6,byrow=TRUE) rownames(mat) <- paste("s",1:6,sep="") colnames(mat) <- paste("g",1:20,sep="") head(mat) ## Espectros com maior correlação pca.object <- prcomp(mat,center=TRUE,scale.=FALSE) plot(pca.object) #Quero os cinco mais correlacionados topN <- 5 load.rot <- scale(pca.object$rotation) names(load.rot[,1][order(abs(load.rot[,1]),decreasing=TRUE)][1:topN]) Então cheguei as variáveis resposta, mas gostaria que alguém me explicasse o por que disso fazendo um favor, pois quando utilizo os meus dados originais, são sempre as variáveis resposta com maior valor numérico que são selecionadas, mesmo utilizando a função scale(). Obrigado, Alexandre -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================

3 7

Escala em gráficos
by Lucas Stempkowski 05 Mai '18

05 Mai '18

Bom dia, Plotei um gráfico boxplot com o seguinte comando: boxplot(prodt~sistema, xlab="Anos sem trigo", ylab="Produtividade kg ha-¹") Como a produtividade foi variável na cultivares testadas, gostaria de padronizar a escala no eixo y para que ficasse dessa forma 1500, 2000, 2500, 3000, 3500, 4000, como faço essa alteração? Muito obrigado. -- ----------------------------------------------------- Lucas Antonio Stempkowski Eng. Agronômo - UFFS Mestrando em Produção Vegetal CAV/UDESC +55 54 99153-0124/ 98448-4633 --------------------------------------------------

2 1