- 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

Re: [R-br] Converter Shape de Pontos para Shape de Polígonos (Yury Duarte)
by Maurício Sangiogo 22 Ago '18

22 Ago '18

Não seria mais fácil tu fazer essa conversão pelo SIG? Eu faço essa conversão no QGis, mas outros programas da área também fazem essa conversão facilmente. Att, Maurício Sangiogo Obter o Outlook para Android<https://aka.ms/ghei36> ________________________________ From: R-br <r-br-bounces(a)listas.c3sl.ufpr.br> on behalf of r-br-request(a)listas.c3sl.ufpr.br <r-br-request(a)listas.c3sl.ufpr.br> Sent: Wednesday, August 22, 2018 2:23:48 PM To: r-br(a)listas.c3sl.ufpr.br Subject: Digest R-br, volume 92, assunto 14 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. Converter Shape de Pontos para Shape de Polígonos (Yury Duarte) ---------------------------------------------------------------------- Message: 1 Date: Wed, 22 Aug 2018 14:23:33 -0300 From: Yury Duarte <yurynepomuceno(a)gmail.com> To: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Converter Shape de Pontos para Shape de Polígonos Message-ID: <CADFS6E5bHfa1JhUawhruRXnRYCFB_e0UEdJwdey-WPLLv_DVyA(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" 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 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180822/648d2750/atta…> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: X.rar Tipo: application/octet-stream Tamanho: 142643 bytes Descrição: não disponível URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20180822/648d2750/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 92, assunto 14 *****************************************

1 0

1. Variáveis intervalares nos modelos mistos do R
by Thiago Wendling gonçalves de oliveira 20 Ago '18

20 Ago '18

Assunto: Variáveis intervalares nos modelos mistos do R Bom dia, meu nome é Thiago e estou precisando de ajuda na avaliação de alguns dados. Quero modelar a densidade básica da madeira de Eucalyptus de acordo com variáveis como: Espécies (3 espécies) Idade (contínua) Tipo de Solo (variável categórica - 7 classes) Textura do Solo (variável intervalar - 3 intervalos)Classes de Altitude (variável intervalar - 5 intervalos)Precipitação (variável inteira)Município (variável categórica - várias classes) e Estado (variável categórica - 3 classes). Eu estudei e penso em utilizar modelos mistos não lineares (do pacote nlme) visto que o comportamento da densidade básica perante as variáveis independentes não segue uma função linear, contudo me deparo com alguns problemas: 1) O primeiro é as variáveis intervalares pois não sei como trabalhar com elas no modelo.2) O segundo é que variáveis escolher para os efeitos aleatórios, eu primeiramente pensei no Estado e dentro dele o Município. Alguém poderia me auxiliar em algum caminho que poderia tomar para resolver esse problema? Quem sabe modelos mistos não sejam o caso e eu poderia utilizar alguma técnica de inteligência artificial, visto que, atualmente, a presença de variáveis intervalares e categóricas é comum em muitas bases de dados. Obrigado. Thiago Wendling Gonçalves de Oliveira Engenheiro FlorestalMestre em Recursos Florestais (ESALQ/USP)Doutorando em Engenharia Florestal (UFPR)Skype: Thiago W.G. de Oliveira+ 55 41 98858-2031

2 1

Ler arquivo JSON com diversos outros JSONs dentro dele
by Yury Duarte 17 Ago '18

17 Ago '18

Boa tarde colegas listeiros, como vão? Ultimamente venho encontrando certa dificuldade para trabalhar com arquivos JSON que contém outros JSON dentro deles. Sempre que precisei trabalhar com esse tipo de arquivo utilizei a expressão: library(jsonlite) raiz= ('C:\\Users\\Desktop\\') a = fromJSON(paste0(raiz,'document.json')) Dessa forma eu obtenho uma lista e consigo selecionar os valores dos campos que quero. Mas quando tento aplicar a mesma lógica para um JSON que contém diversos outros JSONs, o código quebra e retorna: Error in parse_con(txt, bigint_as_char) : parse error: trailing garbage 99809265137 } }, { "_index": "ch (right here) ------^ Meu interesse é acessar apenas as informações de alguns campos desses JSONs e construir uma tabela, onde nas colunas eu terei os campos que desejo e nas linhas eu terei os valores de cada um dos JSONs. A seguir, estou colando no corpo da mensagem um pedaço de um desses arquivos que estou tentando lidar no R. Desde já, agradeço pela ajuda e atenção. Att { "_index": "busca", "_type": "data_point", "_id": "interpolation:699646|mgper-0-20|2018-06-22|5x21sxq1km3b", "_score": 1, "_source": { "name": "Mg% 0-20", "data_type": "interpolation", "attribute": "Mg%", "farm_id": 22706, "depth": "0-20", "location": [ -48.91500982855756, -22.293821593345108 ], "inceres_id": 699646, "date_created": "2018-06-22T14:06:13.587631-03:00", "slug": "mgper-0-20", "quantity": 10.565964698791504 } }, { "_index": "busca", "_type": "data_point", "_id": "interpolation:701809|ph-20-40|2018-06-22|5x21uc94cc55", "_score": 1, "_source": { "name": "ph 20-40", "data_type": "interpolation", "attribute": "ph", "farm_id": 22706, "depth": "20-40", "location": [ -48.90241501215151, -22.28980716917012 ], "inceres_id": 701809, "date_created": "2018-06-22T14:06:35.183784-03:00", "slug": "ph-20-40", "quantity": 4.381259918212891 } }, { "_index": "busca", "_type": "data_point", "_id": "interpolation:699646|mgper-0-20|2018-06-22|5x21sz1rcjp5", "_score": 1, "_source": { "name": "Mg% 0-20", "data_type": "interpolation", "attribute": "Mg%", "farm_id": 22706, "depth": "0-20", "location": [ -48.915268161911925, -22.289490533842073 ], "inceres_id": 699646, "date_created": "2018-06-22T14:06:13.587631-03:00", "slug": "mgper-0-20", "quantity": 8.669662475585938 } }, { "_index": "busca", "_type": "data_point", "_id": "interpolation:701809|ph-20-40|2018-06-22|5x21uc3097d7", "_score": 1, "_source": { "name": "ph 20-40", "data_type": "interpolation", "attribute": "ph", "farm_id": 22706, "depth": "20-40", "location": [ -48.904161243141765, -22.289829735373146 ], "inceres_id": 701809, "date_created": "2018-06-22T14:06:35.183784-03:00", "slug": "ph-20-40", "quantity": 4.352087497711182 } } Yury Duarte Engenheiro Agrônomo - ESALQ/USP

3 8

Como aplicar um modelo em novos dados
by Edmar Caldas 13 Ago '18

13 Ago '18

Boa Tarde! Preciso aplicar um modelo em novos dados. Esse modelo é supervisionado. Novos dados não tem variável target, como proceder nesse caso, fazendo merge dos arquivos? Se fizer merge vai somar cerca de 10 milhões de registros. Edmar

2 1

Fatorial duplo e teste de tukey
by Geraldo Cabral Gouveia 13 Ago '18

13 Ago '18

Bom dia pessoal, estou precisando de uma ajuda Tenho que rodar a estatística de um fatorial duplo em tukey e para as interações quando significativas preciso do desdobramento Meus dados estão tabulados no excel Não tenho idéia de como proceder pois estou começando a trabalhar agora com o R. O fatorial é presença e ausência de um produto aplicado em 4 genótipos de arroz com 3 repetições para cada tratamento e 8 variáveis Obrigado Obter o Outlook para Android<https://aka.ms/ghei36>

2 1

Curvas não lineares - qual o melhor modelo?
by Thallita Grande 10 Ago '18

10 Ago '18

Thallita Grande compartilhou arquivos do OneDrive com você. Para exibi-los, clique nos links abaixo. <https://1drv.ms/u/s!AppBtJtj9_6eoi_LWHEfO01CSJ2u> [https://r1.res.office365.com/owa/prem/images/dc-png_20.png]<https://1drv.ms/u/s!AppBtJtj9_6eoi_LWHEfO01CSJ2u> IICxnativa_p1_2017 1.png<https://1drv.ms/u/s!AppBtJtj9_6eoi_LWHEfO01CSJ2u> <https://1drv.ms/u/s!AppBtJtj9_6eojB0ldbGtpxKiffC> [https://r1.res.office365.com/owa/prem/images/dc-png_20.png]<https://1drv.ms/u/s!AppBtJtj9_6eojB0ldbGtpxKiffC> IIC_fracoes_p1_2017_area 2.png<https://1drv.ms/u/s!AppBtJtj9_6eojB0ldbGtpxKiffC> Tenho essas 4 curvas não lineares entre os índices de conectividade e a porcentagem de área nativa remanescente. Quais são os modelos mais adequados para o ajuste das curvas nas imagens? Meu palpite é que para os 3 gráficos da primeira figura, seria um modelo logístico ou quadrático-platô e na outra figura, sera exponencial ou Michaelis-Manten. Thallita Grande

3 2

Dados de plano altimetria
by Ze Henrique 10 Ago '18

10 Ago '18

Bom dia Senhores, estou tendo que construir uma uma mapa com as altitudes dos municípios brasileiros e não estou conseguindo onde encontrar está informação para que consiga utilizar no R. Alguém conhece onde posso conseguir está informação. Desde já agradeço a todos. -- *José Henrique Soler Guilhen* Doutorando em Genética e Melhoramento Universidade Federal do Espírito Santo - CCAUFES

3 2

ANOVA com medidas repetidas no tempo
by Amanda Martini 09 Ago '18

09 Ago '18

Olá pessoal, Sou relativamento nova no e R e estou trabalhando em análises com medidas repetidas, porém estou insegura com o código que estou utilizando. Nós trabalhamos com avaliação de desempenho animal em diferentes períodos, sendo que nossas repetições são os animais que são aleatorizados nos tratamento e coletadas as informações em diferentes períodos. Para montar o código utilizei o pacote lme4 considerei como fixos: trat e período , aleatório: animais, ficando desse forma: a1 = lme(desempenho ~ trat * per, random = ~1|animal, method="REML",data=desempenho). Estou no caminho certo? Alguém tem alguma sugestão ? Outra duvida é: como eu indico que esses meus animais é dentro do tratamento? Montei o código para indicar essa questão da seguinte forma: a2= lme(desempenho ~ trat * per, random = ~1|animal/trat, method="REML",data=desempenho). Está certo ou estou seguindo um raciocínio errado? Além disso, como descubro qual a matriz de covariância melhor para os meus dados e onde indico no codigo? No SAS sempre realizamos essa analise para posteriormente fazer a analise de variância. Agradeço desde já, Att, Amanda Martini Mestranda em Produção Animal – Forragicultura Laboratório Pastos e Suplementos

1 0

Formatação de summary table no RStudio apresentando &nbsp
by Elias Carvalho 05 Ago '18

05 Ago '18

*Estou testando o pacote "qwraps2", então estou rodando o exemplo em: https://github.com/r-lib/xml2/issues/151 <https://github.com/r-lib/xml2/issues/151>, conforme abaixo.* set.seed(42) library(dplyr) ## ## Attaching package: 'dplyr' ## The following objects are masked from 'package:stats': ## ## filter, lag ## The following objects are masked from 'package:base': ## ## intersect, setdiff, setequal, union library(qwraps2) # define the markup language we are working in. # options(qwraps2_markup = "latex") is also supported. options(qwraps2_markup = "markdown") data(mtcars) mtcars2 <-dplyr::mutate(mtcars, cyl_factor = factor(cyl,levels = c(6, 4, 8), labels = paste(c(6, 4, 8), "cylinders")), cyl_character = paste(cyl, "cylinders")) our_summary1 <- list("Miles Per Gallon" = list("min" = ~ min(mpg), "max" = ~ max(mpg), "mean (sd)" = ~ qwraps2::mean_sd(mpg)), "Displacement" = list("min" = ~ min(disp), "max" = ~ max(disp), "mean (sd)" = ~ qwraps2::mean_sd(disp)), "Weight (1000 lbs)" = list("min" = ~ min(wt), "max" = ~ max(wt), "mean (sd)" = ~ qwraps2::mean_sd(wt)), "Forward Gears" = list("Three" = ~ qwraps2::n_perc0(gear == 3), "Four" = ~ qwraps2::n_perc0(gear == 4), "Five" = ~ qwraps2::n_perc0(gear == 5)) ) our_summary1 ### Overall summary_table(mtcars2, our_summary1) *O código foi executado com sucesso no RStudio, porém a tabela aparece como:* [image: Screenshot from 2018-08-04 18-29-38.png] *Como faço para aparecer a tabela formatada como no exemplo:* [image: Screenshot from 2018-08-04 18-31-38.png] -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

3 2

Dados ausaentes (NA)
by Edson Lira 03 Ago '18

03 Ago '18

Boa tarde, a todos, estou com um problema com essa tabela, dentre outras.Vejam que há 38 observações em branco (seriam NAs). $DispSecNao Qtde Perc 38 79.17DispAcionamento, DispParadaEmergencia, DispInterTrav 1 2.08DispParadaEmergencia 1 2.08DispParadaEmergencia, DispInterTrav 7 14.58DispPartida, DispParadaEmergencia, DispInterTrav 1 2.08Total 48 100.00 Usei a rotina abaixo para eliminar os NAs. b<-read.csv2("C:\\Users\\Edson Lira\\Google Drive\\CONSULTA\\CLAUDIO_SILVIO\\banco1.csv")b[b==""]<-NA Para outras tabelas que existem registro em branco, o NA foi eliminado das frequências, exceto na tabela acima. Alguém tem alguma sugestão? [ ]'s Prof. Edson Lira, Me Estatístico Manaus-Amazonas

3 3

Erro na leitura de dados
by Andre Oliveira 01 Ago '18

01 Ago '18

Boa noite! tenho feito a leitura de dados usando a função fread. Hoje rodei o script e deu este aviso abaixo. Parece ser memória, mas não mudei nada na máquina. System errno 22 unmapping file: Argumento inválido Error in fread("dados.csv", header = TRUE, select = c("ESCOLA", : Opened 4.917GB (5279747593 bytes) file ok but could not memory map it. This is a 64bit process. There is probably not enough contiguous virtual memory available. Se alguém puder ajudar ficarei grato! obg obgAndré

3 3

Re: [R-br] Analise Fatorial - Componentes Principais
by sznelwar＠uol.com.br 31 Jul '18

31 Jul '18

2 1

Analise Fatorial - Componentes Principais
by Edmar Caldas 31 Jul '18

31 Jul '18

Boa Tarde! Preciso de uma ajuda na Análise Fatorial, estou tentando salvar os scores dos fatores para colocar no aquirvo que estou trabalhando, mas não está rolando,estou comparando com outra ferramenta estatística SPSS e os valores são diferentes, somente os scores. Estou tentando salvar pelo método de regression. # AF pelo método componentes principais library(readxl)Olimpiadas <- read_excel("C:/Users/edmar/Desktop/FIAP Atualização BDT/Data Mining/Olimpiadas.xlsx")View(Olimpiadas)dim(Olimpiadas)CP <-Olimpiadas[,-11]CPdim(CP) View(CP) #Pedir as estatisticas descritivas e a matriz de correlação. cor(CP)sapply(CP,summary)sapply(CP,sd) # AF package psychinstall.packages("psych")library(psych) # kmoKMO(CP) # teste de esfericidade de Bartlettcortest.bartlett(CP) install.packages("GPArotation") library(GPArotation) # Tabela de variancia total explicad do spss pelos componentes principaisfit <-princomp(CP, cor=TRUE)fitsummary(fit) # Grafico scrren plotloadings(fit)plot(fit,type="lines") # extrair os fatores sem rotação - component matrix spssfit <-principal(CP,nfactors = 2, rotate = "none", scores = TRUE)fit # rotated componetn matrix spss fit <-principal(CP,nfactors=2,rotate="varimax",scores = TRUE, method = "regression")# EXTRAÇÃO PELO METODO VARIMAX.print(fit,digits = 3,cutoff=.0,sort = TRUE) #comunalidadesfit$communality #Grafico componentes espaço rotacionadoplot(load, type = "n")text(load, labels = names(CP),cex = .7)#nome nas variaveisabline(h=0)abline(v=0) Agora salvar os scores no arquivo não estou conseguindo. Os demais resultados batem todos. estou tentando: Mas sai diferente do SPSS?Tem algum segredo aqui????? > fit$scores RC1 RC2 [1,] -0.42227707 0.967455985 [2,] -1.08912141 0.462360637 [3,] -0.86924034 0.225857490 [4,] -0.53780514 0.919581430 [5,] -1.63914712 -0.616829273 [6,] -1.44129573 -0.088363255 [7,] -0.15298741 0.849870177 [8,] -0.73360090 0.216111850 [9,] -0.38707353 0.452628647[10,] 0.06303505 1.118762527[11,] 0.01554309 1.128912314[12,] -0.54253297 0.038456484[13,] -1.24977051 -0.689280878[14,] -0.83306807 0.006922647[15,] -0.48684252 0.033011053[16,] -0.25825680 0.210099942[17,] 2.08127102 2.190927462[18,] 0.86484820 1.186320909[19,] -0.82669302 -0.930504761[20,] 0.27811485 0.665010689[21,] 0.22335228 -0.092869066[22,] 0.20287071 -0.206901381[23,] 0.36309453 0.240739372[24,] 0.12631659 -0.329892015[25,] -0.17699769 -0.401440678[26,] -0.37557882 -1.074277378[27,] -0.60432737 -1.290110113[28,] 1.44661663 0.853691289[29,] -0.16292001 -1.175375590[30,] 0.78399429 -0.250845135[31,] 1.92828660 0.654161779[32,] 1.13276765 -0.405341447[33,] 0.51556035 -1.656032055[34,] 2.76386458 -3.212819656 Edmar

1 0

Decomposição de Cholesky para matriz não positiva definida. Exemplo geoestatístico
by Wagner Wolff 31 Jul '18

31 Jul '18

Olá pessoal da lista! Estou tentando estimar os parâmetros de um modelo geoestatistico no qual a matriz de covariancia é obtida por meio de uma matriz de distância entre linhas. Ou seja, linhas interceptadas tem distância zero, assim não somente a diagonal da matriz de distância é zero. Isso acaba acarretando em problema na hora de decompor a matriz de covariância usando Cholesky. Alguém tem alguma solucão ou dica? Já tentei algumas coisas como aproximacões para ser positiva definida, entre outras, mas na hora de estimar os parâmetros nada converge. Abraco

3 2

subset a dataframe
by Antonio Silva 28 Jul '18

28 Jul '18

Olá, Algo que deve ser simples mas no qual empaquei. Tenho um vetor com o nome de algumas linhas com dados suspeitos: a <- c(1,2,3,4,5) b <- c(2,8,6,4,10) dat <- data.frame(a,b) out <- c("2","4") # linhas com dados suspeitos Com dat[out,] eu consigo selecionar as linhas do dataframe com dados suspeitos. Mas eu quero selecionar as linhas sem dados suspeitos. Tentei dat[!out,] dat[-out,] e algumas outras coisas sem sucesso. Agradeço sugestões, obrigado. -- Antônio Olinto Ávila da Silva

4 5

Heatmaps para genes
by Paulo Sérgio Pereira de Amorim 25 Jul '18

25 Jul '18

Boa noite! Quais as etapas para elaboração de um Heatmaps para genes. Estou com dificuldades na utilização da função heatmap( ). Abraços.

3 3

Análise de amostra pequena
by Elias Carvalho 17 Jul '18

17 Jul '18

Prezados Tenho uma amostra pequena com apenas 17 registros, e 7 variáveis dicotômicas, cujo resumo é: A B C D E F G 0: 7 0:10 0: 6 0:11 0: 2 0: 1 0:15 1:10 1: 7 1:11 1: 6 1:15 1:16 1: 2 A variável G seria meu desfecho e o resto possíveis preditores Para veficiar a associação entre o desfecho e as outras variáveisexecutei o teste exato de fisher por ser uma amostra pequena, mas o resutados mostraram que o desfecho não depende de nenhuma variável: > fisher.test(data.to.work$G, data.to.work$A) Fisher's Exact Test for Count Data data: data.to.work$G and data.to.work$A p-value = 0.4853 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence interval: 0.1311443 Inf sample estimates: odds ratio Inf > fisher.test(data.to.work$G, data.to.work$B) Fisher's Exact Test for Count Data data: data.to.work$G and data.to.work$B p-value = 0.4853 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence interval: 0.000000 7.625189 sample estimates: odds ratio 0 > fisher.test(data.to.work$G, data.to.work$C) Fisher's Exact Test for Count Data data: data.to.work$G and data.to.work$C p-value = 0.5147 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence interval: 0.1003871 Inf sample estimates: odds ratio Inf > fisher.test(data.to.work$G, data.to.work$D) Fisher's Exact Test for Count Data data: data.to.work$G and data.to.work$D p-value = 0.5147 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence interval: 0.00000 9.96144 sample estimates: odds ratio 0 > fisher.test(data.to.work$G, data.to.work$E) Fisher's Exact Test for Count Data data: data.to.work$G and data.to.work$E p-value = 1 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence interval: 0.01971228 Inf sample estimates: odds ratio Inf > fisher.test(data.to.work$G, data.to.work$F) Fisher's Exact Test for Count Data data: data.to.work$G and data.to.work$F p-value = 1 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence interval: 0.003434375 Inf sample estimates: odds ratio Inf Em seguida calculei o V de Cramer, Coeficiente de Contingência e phi e todos mostram independência e associação fraca. summary(assocstats(tab_cont_G_A)) # Resultado: # X^2 df P(> X^2) # Likelihood Ratio 2.3071 1 0.12878 # P-Valor > 0.05 - independentes # Pearson 1.5867 1 0.20780 # P-Valor > 0.05 - independentes # # Phi-Coefficient : 0.306 # Coef de Phi - Associação fraca # Contingency Coeff.: 0.292 # Coef de Cont - Associação fraca # Cramer's V : 0.306 # V de Crammer - Associação fraca Pelo que pesquisei este seriam os métodos mais adequados para este tipo de análise Há um método melhor para avaliar esse tipo de dado? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

2 3

Re: [R-br] Pacote Chaid
by sznelwar＠uol.com.br 16 Jul '18

16 Jul '18

1 0

Pacote Chaid
by Edmar Caldas 16 Jul '18

16 Jul '18

Boa Tarde! Pessoal, não estou conseguindo executar esse código, podem me ajudar por favor? attach(credito1)str(credito1) install.packages("partykit")install.packages("CHAID", repos="http://R-Forge.R-project.org") library(partykit)library(CHAID) credit_rank<-as.factor(credit_rank)age <-as.factor(age)amex<-as.factor(amex)pay_week <-as.factor(pay_week)class <-as.factor(class) is.factor(credit_rank)is.factor(age)is.factor(class)is.factor(pay_week)is.factor(amex) > modelo <-chaid(credit_rank ~ ., data = credito1)Error: is.factor(x) is not TRUE

2 2

Eixo x
by Fabiana Goncalves Barbosa 16 Jul '18

16 Jul '18

Olá Pessoal Fiz um gráfico (plot) com os dados abaixo, porém no eixo x só apareceram alguns anos. Quero que apareçam todos os anos no eixo x. Como faço? qtidade=c(61,53,44,35,24,17,15,16,18,12,8,3,3,4,2,2,1,2,1,1,1,1,1) Anos=c(2018,2017,2016,2015,2014,2013,2012,2011,2010,2009,2008,2007,2006,2005,2003,2002,2001,1999,1998,1997,1992,1991,1981) plot(Anos,qtidade) Grata Fabiana

3 2

Modelos Mistos e Lattice (Grafico de predições)
by Fernando Souza 14 Jul '18

14 Jul '18

Caros colegas Desejo fazer um gráfico, semelhante ao retornado pela função augPred{nlme}, contendo os valores observados, os valores preditos pelo modelo fixo, e os valores individuais preditos. como esse exemplo (Só que para modelo não linear) https://rdrr.io/cran/nlme/man/plot.augPred.html por algum motivo a função augPred() não está funcionando com o meu modelo (Acredito por ser desbalanceado ) A função predict.nlme {nlme}, me retorna um dataframe com as informação que quero, porém como transformá-los em gráfico? Sei que é possível faze-lo através do xyplot{lattice} porém não estou conseguindo, devido a minha limitação no uso desse pacote. Alguém poderia me ajudar? segue CMR packages <- c("tidyr", "nlme""latticeExtra","lattice") #Instala pacotes necessários, se ja tiver instalado ativa-o for(p in packages){ if(!require(p,character.only = TRUE)) install.packages(p) library(p,character.only = TRUE) } ##==================================================================================================== ## **Preparação do banco de dados** ##==================================================================================================== DADOS.CURTO <- read.csv("https://www.dropbox.com/s/b4cckybgrcg86me/galinhas.csv?raw=1",head = TRUE) colnames(DADOS.CURTO) <- c("0","1","2","3","4","5","6","TRAT") DADOS.CURTO <-DADOS.CURTO[with(DADOS.CURTO,order(TRAT)),] DADOS.CURTO$ID <- c(1:93,94:280,281:370) DADOS<- DADOS.CURTO %>% gather('0','1','2','3','4','5','6',key ="SEMANA",value = "PESO", -TRAT,-ID) DADOS <- transform(DADOS, SEMANA = as.numeric(SEMANA),ID=factor(ID)) DADOS ##-------------------------------------------------------------------- GALINHA <- groupedData(PESO~SEMANA |ID,order.groups=FALSE,data = DADOS) ##------------------------------------------------------------------------------------------- modelo1.list <- nlsList(PESO ~ a * exp(-b * exp( -c * SEMANA)),start = c(a = 4.10,b = 4.53,c = 0.37 ),na.action=na.omit,data =GALINHA) ##--------------Definição do modelo de efeito Fixo-------------------------------- controle2<-nlmeControl(maxIter=300,msMaxIter = 300,niterEM= 100) modelo.nlme.0<-nlme(modelo1.list) modelo.nlme.2<-update(modelo.nlme.0,fixed=list(a~TRAT,b + c~1),start=c(4.10,0,0,4.77,0.34),control=controle2) ##------------------------------------------------------------------------------------------------------ controle3<-nlmeControl(maxIter=300,msMaxIter=1000) MM3.3<- update(modelo.nlme.2,weights = varExp(form = ~fitted(.)|SEMANA),control=controle3,correlation = corARMA(form = ~SEMANA |ID,q = 4)) summary(MM3.3) ##----------------------------------------------------------------- DADOS.ORIGINAL <- getData(MM3.3) predito<-predict(MM3.3) ##------------------------------------------------------------------ Sent from Mailspring (https://getmailspring.com/) the best free email app for work

2 1

Ajuda com Realocação de Inteiros
by Vanessa Felix 13 Jul '18

13 Jul '18

Pessoal, Será que podem me ajudar a pensar em um problema. Tenho um grupo de pessoas que preciso redistribuir de acordo com uma porcentagem, mas não estou conseguindo distribui-las de forma mais justa e que não desequilibre muito a relevância. Segue um exemplo da minha base de dados: Relevância Total Grupo Nova Redistribuição Arredondamento Arredondamento para baixo Arredondamento para cima 0,90% 80 0,721 1,0 - 1,0 1,80% 80 1,441 1,0 1,0 2,0 25,23% 80 20,180 20,0 20,0 21,0 27,93% 80 22,342 22,0 22,0 23,0 44,14% 80 35,315 35,0 35,0 36,0 SOMA 100% 80 79 78 83 relev=c(0.009009009,0.018018018,0.252252252,0.279279279,0.441441441) n=80 relev*n sum(relev*n) ceiling(relev*n) sum(ceiling(relev*n)) floor(relev*n) sum(floor(relev*n)) round(relev*n) sum(round(relev*n)) -- ------------------------------------------------------------------------------------------ Vanessa Felix do Nascimento Sergio Graduada em Estatística - ICE/UFJF Mestre em Saúde Coletiva - IMS/UERJ http://lattes.cnpq.br/2516291371035474 (21) 97557-3956 -------------------------------------------------------------------------------------------

2 1

Uniao por coluna - chave
by Edmar Caldas 12 Jul '18

12 Jul '18

Bom Dia! Preciso de outra ajuda. Estou tentando unir 3 arquivo usando a variavel chave id e não estou conseguindo, ele só uni e arquivos. chave<-merge(Apend_A,Apend_B,Apend_C, by.x=c("id"), by.y=c("id"), by.z=c("id"), all.x=TRUE, all.y=FALSE, all.z=FALSE) chave Edmar

1 0

Soma duas colunas pela função apply
by Edmar Caldas 11 Jul '18

11 Jul '18

Boa Tarde! Estou tentando somar duas colunas e coloca-lá na variável soma_debitos e não estou conseguindo. Não pretendo usar função! > Banco$soma_debitos <-apply(Banco$cartao_credito,Banco$emprestimos,2,sum)Error in match.fun(FUN) : '2' não é uma função, caracter ou símbolo Edmar

3 6

Dados DIC
by fúlvia cristina oliveira 10 Jul '18

10 Jul '18

Boa noite, estou rodando uns dados em DIC, onde tenho dois tratamentos (dietas) e algumas variáveis a saída saiu assim:O que tenho que fazer? > crd(TRAT, Hemoglobina, quali = FALSE, mcomp = "tukey", sigT = 0.05, sigF = 0.05)------------------------------------------------------------------------Analysis of Variance Table------------------------------------------------------------------------ DF SS MS Fc Pr>FcTreatament 1 59.806 59.806 14.495 0.00081149Residuals 25 103.150 4.126 Total 26 162.956 ------------------------------------------------------------------------CV = 25.89 % ------------------------------------------------------------------------Shapiro-Wilk normality testp-value: 0.07082243 According to Shapiro-Wilk normality test at 5% of significance, residuals can be considered normal.------------------------------------------------------------------------ Adjustment of polynomial models of regression------------------------------------------------------------------------Error in treat^2 : non-numeric argument to binary operator> Fúlvia Cristina Oliveira Bióloga - UFMS Mestre em Zootecnia - UEMSDoutoranda em Ciência Animal - UFMS(67) 99610-2812 / 99129-6612

2 1

Construção de gráfico em ggplot2
by Carlos Peixoto Neto 09 Jul '18

09 Jul '18

Olá, busco um auxílio para aprimorar um script em R. Tenho alguns erros de execução que apresento logo abaixo. No estudo, avalio o efeito da interação entre girinos, var. preditora "INT", no uso de recurso espacial através de duas variáveis resposta: deslocamento em Posição Vertical (PV) e Posição Horizontal (PH). Erros: - Na execução da função summarySE: "/Warning message:In qt(conf.interval/2 + 0.5, datac$N - 1) : NaNs produced" - Na execução de legend.title=element_blank(): Error: unexpected symbol in:"legend.title=element_blank() legend.title" #Importar dados. Pacotes para o gráfico pmed<-read.table("https://drive.google.com/open?id=1cET5cROSb-_D-cZYud3yKBHGR-qamind",header=T, sep=',') require(Rmisc) require(ggplot2) require(EnvStats) require(sciplot) require(dplyr) summarypmed <- summarySE(pmed, measurevar="PV", groupvars=c("ID","INT"),na.rm=T)# gráfico 1 graf1 <- ggplot(summarypmed, aes(x=INT, y= PV, fill=INT)+ facet_wrap(~INT) + geom_dotplot(binwidth=0.05,binaxis="y",stackdir = "center") + geom_errorbar(aes (ymin=PV.y-se,ymax=PV.y+se),width = 0.25,size=0.25)+ xlab("Tratamentos") + ylab(" Posição Vertical (cm)") + geom_text(aes(label = paste("N", "==",N,sep = "")), parse = TRUE, y=-0.15) + geom_point(aes(y=PV.y), size=1,show.legend = F) + theme_bw () + scale_fill_manual(values=c("grey85","grey30")) + theme(axis.ticks.x=element_blank(), axis.title.x=element_blank(), axis.text=element_text(size=7), axis.text.y=element_text(size=10), axis.text.x=element_blank(), axis.title.y=element_text(size=11), plot.title=element_text(size=12, hjust = 0.5), panel.grid=element_blank(), plot.title=element_text(vjust=2), legend.title=element_blank()) *Carlos Alberto Algarves Peixoto Neto* Biólogo | Universidade Federal do Maranhão - UFMA Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.do? id=K4303028A6 ___ "O amor de Deus: Cristo morrer em nosso favor quando ainda éramos pecadores." Romanos cap. 5: vers. 8 (NVI)

1 0

Convite para a comunidade de R: submissões de propostas para o Mozilla Festival 2018
by Marcos Vital 29 Jun '18

29 Jun '18

Oi, pessoal da lista! A nona edição do Mozilla Festival (o MozFest) será realizada em Londres, sexta-feira, 26 de outubro a domingo, 28 de outubro. O festival é um evento super diverso, que tenta acolher pessoas do mundo todo que estejam de alguma forma interessadas em uma internet livre e acessível para todas as pessoas - isso passa por vários aspectos, que incluem iniciativas ligadas à Open Science e Open Source, de forma que a comunidade ligada ao R pode se interessar. No ano passado eu tive a chance de participar como facilitador de uma sessão sobre ensino de R, e neste ano sou um dos "wranglers", a equipe que organiza o festival. No momento a chamada para sessões está aberta, e vocês podem submeter uma proposta aqui: https://mozillafestival.org/proposals As inscrições de sessões ficam abertas até o dia 1o de Agosto - depois disso só será possível se inscrever como participante. Há recursos para que algumas propostas de sessões tenham ajuda de custo, no caso de candidatos que não tenham condições de viajar para Londres e se hospedar lá sem isso. Caso alguém queira saber mais, fiquem à vontade para entrar em contato comigo, ok? Abraços! Marcos -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

1 0

Eliminação de observações que não variam ao longo do período
by Alexandre Loures 28 Jun '18

28 Jun '18

Boa noite! Tenho uma base com a seguinte característica: Ano origem destino HS04 valor 1989 ARG BRA 0101 1 1989 ARG BRA 0102 0 1989 ARG BRA 0103 0 1989 ARG BRA 0104 1 2010 ARG BRA 0101 1 2010 ARG BRA 0102 1 2010 ARG BRA 0103 1 2010 ARG BRA 0104 1 Preciso eliminar as observações HS04 que não variaram ao longo do período. Ou seja, HS04 0101 e a HS04 0104 pois tanto no ano inicial quanto no final tinham valor 1. Desde já, muito obrigado pela atenção! -- *Alexandre Loures* Bolsista de Pós-doutorado Júnior do CNPq Centro de Ciências Sociais Aplicadas Universidade Federal da Paraíba www.ccsa.ufpb.br/ppge/ www.rodriguesloures.com

1 0

Faixa etária
by Elias Carvalho 27 Jun '18

27 Jun '18

Boa noite pessoal Qual a melhor forma de categorizar a variável "idade"? Existe uma pré-definição de faixa etária no Brasil ? Coisa do tipo: criança, jovem, adulto, idoso? -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

3 4

Equalizar duas colunas em dois data frames a fim de realizar posterior merge
by Jobenil Luiz Magalhães Júnior 26 Jun '18

26 Jun '18

Boa tarde pessoal, Estou tentando realizar o merge em dois dataframes, porém as respectivas colunas são ligeiramente diferentes: df1 <- (concatenado = c("1 _ 307", "22 _ 105", "12 _ 303", "6 _ 201"), var1 = 1:4, var2 = rnorm(4)) df2 <- (concatenar = c("500002_01 _ 307", "5001245_05 _ 303", "500033_22 _ 105", "501256_12 _ 303", "500345_06 _ 201", var3 = runif(5), var4 = rnorm(5)) Quero usar os 8 últimos caracteres da variável df2$concatenar e comparar com as da variável df1$concatenado. Noto que tenho dois problemas: 1) completar com 0 "zeros" à esquerda os valores à esquerda em df1$concatenado, quando forem menores do que 10. 2) Extrair os 8 caracteres à direita da variável df2$concatenar. Como posso equalizar as duas colunas? Agradeço desde logo qualquer ajuda que puderem dispensar. Atenciosamente Jobenil Júnior [image: Mailtrack] <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=sig…> Remetente notificado por Mailtrack <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=sig…> 26/06/18 11:37:17

2 1

Transformação de dados em delineamento fatorial
by fúlvia cristina oliveira 22 Jun '18

22 Jun '18

Eu tenho alguns dados que rodei em delineamento fatorial, com dois fatores. Usei o pacote Expdes e o teste do Shapiro dos dados estão normais.Como transformo esses dados? Em DIC eu uso as transformações em Log, raiz quadrada e outros. Mas em fatorial qual valores que transformo? Do fator 1 ou fator 2.? Fúlvia Cristina Oliveira Bióloga - UFMS Mestre em Zootecnia - UEMSDoutoranda em Ciência Animal - UFMS(67) 99610-2812 / 99129-6612

3 2

dúvida no shiny
by Luiz Leal 21 Jun '18

21 Jun '18

PrezadosEstou rodando o código abaixo que peguei na internet, mas não funciona. Alguém saberia o motivo?Desde já agradeço.Luiz require(shiny) runApp( list(ui = pageWithSidebar( headerPanel(withMathJax("$$\\text{Display formula in heading }X_n=X_{n-1}-\\varepsilon$$")), sidebarPanel( ), mainPanel( ) ), server= function(input, output, session){ } ) )

2 1

Transformação dos dados
by fúlvia cristina oliveira 20 Jun '18

20 Jun '18

Boa noite Como eu faço para transformar dados que estão em esquema fatorial (com 2 fatores)? Fúlvia Cristina Oliveira Bióloga - UFMS Mestre em Zootecnia - UEMSDoutoranda em Ciência Animal - UFMS(67) 99610-2812 / 99129-6612

2 1

Erro na função da - asremlPlus
by Andre Oliveira 19 Jun '18

19 Jun '18

Boa tarde, estou com problema neste modelo. require(asremlPlus) data(Wheat.dat) current.asr <- asreml(yield ~ Rep + WithinColPairs + Variety, random = ~ Row + Column + units, rcov = ~ ar1(Row):ar1(Column), data=Wheat.dat) Error in asreml(yield ~ Rep + WithinColPairs + Variety, random = ~Row + : could not find function "asreml" Alguém saberia dizer alguma dica de solução? André Oliveira Souza. obg

1 0

Mauro Sznelwar compartilhou The R Book.pdf com você
by Mauro Sznelwar (via Dropbox) 17 Jun '18

17 Jun '18

Olá, Mauro Sznelwar (sznelwar(a)uol.com.br) convidou você para visualizar o arquivo " The R Book.pdf " no Dropbox. Visualizar arquivo[1] Aproveite! Equipe Dropbox Mauro e outras pessoas poderão ver quando você visualizar este arquivo. Outros arquivos compartilhados com você pelo Dropbox também podem mostrar esta informação. Saiba mais[2] em nossa Central de ajuda. [1]: https://www.dropbox.com/l/scl/AABKkLrU7BrGQMQ5jbYZa44eonjq2jqi1z0 [2]: https://www.dropbox.com/l/AABspd8wT0fNffNbgXZilx5XPqt2GJo3JIo

3 2

Mauro Sznelwar compartilhou The R Book.pdf com você
by Mauro Sznelwar (via Dropbox) 15 Jun '18

15 Jun '18

Olá, EDSON, Mauro Sznelwar (sznelwar(a)uol.com.br) convidou você para visualizar o arquivo " The R Book.pdf " no Dropbox. Visualizar arquivo[1] Aproveite! Equipe Dropbox [1]: https://www.dropbox.com/l/scl/AAB4hkRYDcrJ2TMz4AuhQrQaqaNzDeqa3QA

1 0

Re: [R-br] Planejamento e Análise de Experimentos
by sznelwar＠uol.com.br 15 Jun '18

15 Jun '18

2 1

Re: [R-br] Planejamento e Análise de Experimentos
by sznelwar＠uol.com.br 15 Jun '18

15 Jun '18

4 4

Planejamento e Análise de Experimentos
by Emerson Cotta Bodevan 14 Jun '18

14 Jun '18

Prezados, boa tarde. Estou estudando o assunto e gostaria de fazê-lo utilizando o R. Alguém poderia indicar um bom livro/material sobre planejamento e análise de experimentos (preferencialmente utilizando o R) voltado para área de Farmácia? Agradeço qualquer indicação. Att., *Emerson*

1 0

Uma solicitação
by Marcele de Jesus Correa 13 Jun '18

13 Jun '18

Olá, bom dia pessoal! Gostaria de saber se alguém tem um dos livros abaixo que possa compartilhar, porque tentei baixá-los mas, sem sucesso. Caso alguém tenha outros que possa compartilhar ou sugerir, aceito e desde já agradeço. A Modern Approach to Regression with R *ou* Modern Regression Techniques Using R: A Practical Guide. Att, 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.

2 2