- R-br - DINF UFPR

SOMA DE LINHAS
by Edson Lira 27 Set '18

27 Set '18

Bom dia a todos! Tenho um banco aonde da coluna 19 a coluna 29 tenho respostas do tipo: Concordância, Indiferente e Discordância, com um banco de 29 linhas. O meu objetivo é somar por linha (por cada indivíduo, a quantidade e o percentual de Concordância), ou seja acrescentar mais duas colunas por indivíduo. Abaixo as respostas de um indivíduo (uma linha). Concordância Concordância Concordância Concordância Concordância Concordância Concordância Discordância Concordância Concordância Concordância Discordância Indiferente Indiferente Indiferente Concordância Concordância Discordância Indiferente Indiferente Indiferente [ ]'sProf. Edson Lira, Me Estatístico Manaus-Amazonas

3 4

"variance function coefficient" pacote nlme
by Luiz Leal 18 Set '18

18 Set '18

Prezados, boa tardeAlguém sabe como se estima o "variance function coefficient" no pacote nlme?Desde já agradeçoLuiz

1 0

Perdoem-me. Preciso do e-mail de Benilton Carvalho
by Luiz Roberto Martins Pinto 18 Set '18

18 Set '18

Caros amigos, Perdoem-me solicitar este favor, mas preciso do e-mail do Benilton Carvalho. Já enviei vários e-mails mas parece que ele não utiliza mais os que conheço. Abraços, 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 "*The s**cience exists because there are patterns. * * The patterns exist because God created them*. * The statistic exists to research the patterns that God created.*"

1 0

Aplicar msm intervalos para 2 rasters
by Jackson Rodrigues 17 Set '18

17 Set '18

Bom dia pessoal! Estou gerando mapas de temperaturas mínimas e máximas e gostaria de usar apenas 1 intervalo de valores/escala de cores que se aplique tanto a às mínimas quanto às máximas. O problema é que raster é S4 e o intervalo (que criei) é vetor. A menor temperatura é 3 e a maior é 31. Alguém pode me ajudar associar os valores do raster ao intervalo que criei? cols<-colorRampPalette(c("royalblue","springgreen","yellow","orange","red"))(29) intervalo = seq(from=3, to=31) # criando o vetor de cores cores1 <- cols[findInterval(TMin_masked, vec = intervalo)] Error in as.double(x) : cannot coerce type 'S4' to vector of type 'double' Obrigado Jackson

3 2

Generelized least squares - Random effects model
by Luiz Leal 11 Set '18

11 Set '18

PrezadosEstou usando o pacote nlme e estimando o modelo utilizando a função gls().O modelo é de efeitos aleatórios e preciso estimar as componentes de variância (variance between treatments e variance within treatments), mas não sei como.Alguém poderia me ajudar?Desde já agradeçoLuiz

1 0

Dúvidas de SIG no R
by Yury Duarte 06 Set '18

06 Set '18

Boa tarde colegas listeiros, como vão? Há algum tempo atras recorri a lista com algumas duvidas de manipulação do SIG com o R e me levaram bem longe, mas ainda sigo com algumas dificuldades. Estou fazendo uma analise de agrupamento para um local com base em informações georreferenciadas e espacializadas numa grade de 20x20 metros. Inicialmente tenho essas informações num shape de pontos e também tenho o shape do contorno da minha área. O resultado da minha analise de agrupamento é uma lista indicando o cluster de cada conjunto lat;long. Após esse passo, gostaria de criar um filtro que fosse capaz de identificar e substituir o valor de um pixel caso ele estiver cercado por valores de cluster diferente. existe alguma ferramenta que faça algo do tipo? Paralelo a isso, gerei um shape de polígonos com a mesma dimensão da grade de pontos (20x20 metros). Nesse momento ocorre a quebra do script, pois o número de polígonos gerado é menor do que o número de informações que eu tenho no meu shape de pontos original e na minha lista de clusters. Minha ideia foi de fazer um merge (ou coisa do tipo) entre minha lista de cluster e uma lista com os centroides de cada polígono, assim, eu eliminaria os pontos "excedentes" em relação aos polígonos, mas não consegui fazer com que funcionasse. Segue script ao final da mensagem. Desde já, agradeço pela ajuda de todos! library(spatstat) library(maptools) library(shapefiles) library(foreign) library(magrittr) library(formattable) library(e1071) library(rgeos) library(rgdal) library(rgeos) library(raster) library(sp) library(RcppCNPy) library(deldir) library(dismo) library(dplyr) library(ggplot2) library(gstat) #------Ler shape de contorno e shape de dados + transformar coordenadas em UTM------# #-----Lendo shape de dados------# shapefile_dados_das_camadas_utilizadas = readOGR(paste0('diretorio', 'dados.shp')) #-----Lendo shape de contorno------# shapefile_contorno_do_talhao = readOGR(paste0('diretorio', 'contorno.shp')) #------Clusterizacao------# shapefile_dados_padronizados = scale(shapefile_dados_das_camadas_utilizadas@data) shapefile_dados_geolocalizacao = shapefile_dados_das_camadas_utilizadas@coords wss = (nrow(shapefile_dados_padronizados)-1)*sum(apply(shapefile_dados_padronizados,2,var)) cluster_maximo = 30 for (i in 1:cluster_maximo) wss[i] = sum(kmeans(shapefile_dados_padronizados, centers=i)$withinss) plot(1:cluster_maximo, wss, type="b", xlab="Número of Clusters", ylab="Soma dos quadrados dentro dos clusteres") #dendograma = shapefile_dados_padronizados %>% dist %>% hclust #plot(dendograma) #rect.hclust(dendograma, k = 3, border = "blue") #rect.hclust(dendograma, k = 4, border = "red") escolher_num_cluster = 6 lista_clusteres = data.frame((kmeans(shapefile_dados_padronizados, centers = escolher_num_cluster)$cluster),shapefile_dados_geolocalizacao) #------Transformando em poligono------# x = shapefile_dados_das_camadas_utilizadas@coords[,1] y = shapefile_dados_das_camadas_utilizadas@coords[,2] df_xy = data.frame(shapefile_dados_das_camadas_utilizadas@coords [,1],shapefile_dados_das_camadas_utilizadas@coords[,2]) #coordinates(df_xy) <- ~x+y z = deldir(x,y,plot=TRUE,wl='triang',wp="none") w = tile.list(z) polys = vector(mode="list", length=length(w)) for (i in seq(along=polys)) { pcrds = cbind(w[[i]]$x, w[[i]]$y) pcrds = rbind(pcrds, pcrds[1,]) polys[[i]] = Polygons(list(Polygon(pcrds)), ID=as.character(i)) } SP = SpatialPolygons(polys) plot(SP) celulas_dentro_do_talhao = raster::intersect(SP, shapefile_contorno_do_talhao) plot(celulas_dentro_do_talhao) #------Adicionando Informacao de cluster ao poligono------# length(celulas_dentro_do_talhao) confirmando_pontos_dentro_da_area = data.frame(coordinates(celulas_dentro_do_talhao)) length(lista_clusteres[,1]) colnames(lista_clusteres) = c('Cluster', 'Long', 'Lat') colnames(confirmando_pontos_dentro_da_area) = c('Long', 'Latitude') #confirmando_pontos_dentro_da_area$Long = round(confirmando_pontos_dentro_da_area$Long, 5) #lista_clusteres$Long = round(lista_clusteres$Long, 5) #teste = merge(confirmando_pontos_dentro_da_area, lista_final, all.x = TRUE) #celulas_dentro_do_talhao@data$Cluster = '???' plot(celulas_dentro_do_talhao) writeOGR(celulas_dentro_do_talhao, 'C:\\Users\\Yury\\Desktop\\InCeres\\UGD', layer = 'sss', driver = 'ESRI Shapefile') Yury Duarte Engenheiro Agrônomo - ESALQ/USP

1 1

erro: funcao autoplot (ggplot2)
by Alex Brito 03 Set '18

03 Set '18

Prezados, sou iniciante no R, e ao elaborar uma apresentação com o pacote ggplot2, através da função autoplot, apareceu esta mensagem "Error: Objects of type ts not supported by autoplot." Achei estranho pq esse erro. Segue abaixo descrição do script e também o aquivo anexo que pode ser executado em qualquer RStudio. Agradeço, de antemão, a ajuda de vocês!! # Extrair dados do SIDRA-IBGE, a partir do API - Pequisa contas trimestrais - CNT ### tabelas 1621(c/ajuste sazonal) # library(sidrar) tabela1621=get_sidra(api = "/t/1621/n1/all/v/all/p/all/c11255/90687,90691,90696,90707,93404,93405,93406,93407,93408/d/v584%202") # criar vetores para as séries coletadas # series = c(90687,90691,90696,90707,93404,93405,93406,93407,93408)names = c("Agro", "Ind", "Serv", "PIB", "Consumo", "Governo", "FBKF", "Exportacao", "Importacao") # construir uma matriz onde as variáveis dispostas em linhas nas tabelas correspondam a colunas na matriz e os valores das linhas correspondam às observações# pib_sa = matrix(NA, ncol = length(series), nrow = nrow(tabela1621)/length(series)) ## fazer loop para levar as observações de cada série da tabela para a matriz ## for (i in 1:length(series)) { pib_sa[,i]<- tabela1621$Valor[tabela1621$`Setores e subsetores (Código)` ==series[i]] } # dar nomes às colunas # colnames(pib_sa) = names ## transformar em séries temporais ## pib_sa <- ts(pib_sa, start=c(1996,01), freq=4) ## apresentação gráfico ## library("ggplot2")autoplot(pib_sa[,4]) "Error: Objects of type ts not supported by autoplot." ## Alex BritoDepartamento de Economia - UFMAGrupo de Estudos e Pesquisa em Macroeconomia - GRAMMACidade Universitária Dom Delgado, Av. dos Portugueses, 1966, Bacanga,São Luís - MA, 65065-545

2 1

sobre conjuntos
by Antonio Silva 31 Ago '18

31 Ago '18

Car@s Gostaria de saber se há alguma função no R, em algum pacote, que retorne as interseções e valores únicos de diversos conjuntos ou grupos de subconjuntos. Por exemplo para: SP1 <- c("a","b","c","d","g","h","i") SP2 <- c("a","b","c","d","e","f","j") SP3 <- c("a","b","g","h","e","f","k") teríamos a saída: SP1: i SP2: j SP3: k SP1-SP2: c,d SP1-SP3: g,h SP2-SP3: e,f SP1-SP2-SP3: a,b Encontrei funções como intersect, setdiff e Reduce, Nem mesmo os valores exclusivos de cada conjunto eu consegui obter. Agradeço qualquer indicação. Obrigado, -- Antônio Olinto Ávila da Silva Instituto de Pesca São Paulo, Brasil

3 3

Resposta binária
by ana paula coelho madeira 28 Ago '18

28 Ago '18

Prezados colegas, Boa tarde! Estou avaliando um experimento com 5 tratamentos e a varável resposta é binária ( conta a presença ou ausência de uma dada característica). Sei (espero não estar enganada) que se trata de uma regressão logística, visto que a variável reposta é binária e a explicativa categórica. Não estou conseguindo tratar esses dados no R. Encontrei um script no livro do Crawley porém não consegui adaptar. Peço indicações. Desde já agradeço Att., [https://ipmcdn.avast.com/images/icons/icon-envelope-tick-round-orange-anima…]<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…>.

3 2

Converter Shape de Pontos para Shape de Polígonos
by Yury Duarte 24 Ago '18

24 Ago '18

Boa tarde colegas listeiros, como vão? Estou começando a trabalhar com dados georreferenciados e geoestatistica no R e me deparei com uma dificuldade que não estou conseguindo contornar. Tenho um shapefile de pontos que eu gostaria de transformar em um shape de polígonos e não estou sabendo como fazer essa conversão. Depois de fazer essa transformação de pontos para polígonos, queria saber se também é possível eu unir os polígonos adjacentes quando seus valores forem iguais, afim de reduzir o numero total de polígonos do meu shapefile. Desde já, agradeço pela ajuda e atenção de todos! Att Yury Duarte Engenheiro Agrônomo - ESALQ/USP

3 6