- R-br - DINF UFPR

subset de observações específicas de uma variável
by cayanportela＠hotmail.com 08 Jul '16

08 Jul '16

is.factor(dados1$HScode) [1] TRUE dados1$HScode <- as.character(dados1$HScode) dados2 <- subset(dados1,nchar(dados1$HScode)==4) Caso a sua variavel esteja definida como fator, transforme para character e use o nchar para contagem do numero de caracteres. <quote author='R-br mailing list'> Boa noite! Tenho uma base de dados em que uma das variáveis (HScode) possui três níveis. Por exemplo, 01 0101 010122 0102 010233 02 0201 020179 Contudo preciso de um subset apenas com as observações com quatro dígitos. Ou seja, 0101 0102 0201 structure(list(reporter = c("Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania"), HScode = c("01", "0101", "010110", "010190", "0102", "010210", "010290", "0103", "010310", "010391", "010392", "0104", "010410", "010420", "0105", "010511", "010512", "010519", "010594", "010599", "0106", "010611", "010612", "010619", "010620", "010631", "010632", "010639", "010690", "02", "0201", "020110", "020120", "020130", "0202"), boundtariff = c("6.571969696969697", "8.75", "7.5", "10", "7.291666666666667", "5", "9.583333333333334", "5", "5", "5", "5", "5", "5", "5", "9.5", "5", "5", "5", "12.5", "20", "5", "5", "5", "5", "5", "5", "5", "5", "5", "10.27914614121511", "10", "10", "10", "10", "10")), datalabel = "", time.stamp = " 7 Jul 2016 22:53", .Names = c("reporter", "HScode", "boundtariff"), formats = c("%54s", "%9s", "%20s"), types = c(54L, 7L, 20L), val.labels = c("", "", ""), var.labels = c("", "", ""), version = 12L, row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35"), class = "data.frame") Alguém saberia como fazer isso? Desde já muito obrigado! -- *Alexandre Rodrigues Loures* Doutorando em Economia Aplicada Universidade Federal da Paraíba - UFPB Centro de Ciências Sociais Aplicadas - CCSA Programa de Pós-Graduação em Economia - PPGE Site: www.ccsa.ufpb.br/ppge orcid www.orcid.org/0000-0002-1288-0135 _______________________________________________ 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. iD icon.gif (1K) <http://r-br.2285057.n4.nabble.com/attachment/4666373/0/iD%20icon.gif> </quote> Quoted from: http://r-br.2285057.n4.nabble.com/R-br-subset-de-observacoes-especificas-de… _____________________________________ Sent from http://r-br.2285057.n4.nabble.com

1 0

subset de observações específicas de uma variável
by Alexandre Loures 08 Jul '16

08 Jul '16

Boa noite! Tenho uma base de dados em que uma das variáveis (HScode) possui três níveis. Por exemplo, 01 0101 010122 0102 010233 02 0201 020179 Contudo preciso de um subset apenas com as observações com quatro dígitos. Ou seja, 0101 0102 0201 structure(list(reporter = c("Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania", "Albania"), HScode = c("01", "0101", "010110", "010190", "0102", "010210", "010290", "0103", "010310", "010391", "010392", "0104", "010410", "010420", "0105", "010511", "010512", "010519", "010594", "010599", "0106", "010611", "010612", "010619", "010620", "010631", "010632", "010639", "010690", "02", "0201", "020110", "020120", "020130", "0202"), boundtariff = c("6.571969696969697", "8.75", "7.5", "10", "7.291666666666667", "5", "9.583333333333334", "5", "5", "5", "5", "5", "5", "5", "9.5", "5", "5", "5", "12.5", "20", "5", "5", "5", "5", "5", "5", "5", "5", "5", "10.27914614121511", "10", "10", "10", "10", "10")), datalabel = "", time.stamp = " 7 Jul 2016 22:53", .Names = c("reporter", "HScode", "boundtariff"), formats = c("%54s", "%9s", "%20s"), types = c(54L, 7L, 20L), val.labels = c("", "", ""), var.labels = c("", "", ""), version = 12L, row.names = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35"), class = "data.frame") Alguém saberia como fazer isso? Desde já muito obrigado! -- *Alexandre Rodrigues Loures* Doutorando em Economia Aplicada Universidade Federal da Paraíba - UFPB Centro de Ciências Sociais Aplicadas - CCSA Programa de Pós-Graduação em Economia - PPGE Site: www.ccsa.ufpb.br/ppge orcid www.orcid.org/0000-0002-1288-0135

3 2

Extrair dados de ajuste
by Pedro Campelo 05 Jul '16

05 Jul '16

Pessoal, Estou rodando um script e preciso extrair os dados dos coeficientes ajustados para um modelo. No meu exemplo, quero extrair e separar os valores "k" e "n". *> data.frame(A)* * M R* *1 10 24.61* *2 20 29.41* *3 30 33.07* *4 40 37.37* *5 50 38.22* *6 60 41.14* *7 70 41.61* *8 80 43.41* *9 90 47.18* *10 100 48.77* *11 110 51.02* *12 120 52.95* *13 130 52.98* *14 140 52.92* *> model<-nls(R~k*M^n, start=list(k=0.5,n=1))* *> summary(model)* Obrigado *----* *Pedro Henrique Campelo Felix* *Bachelor of Food Engineering (UNI-BH)* *Msc. Chemical Engineering (UFMG)* *PhD Student of Food Science (UFLA)* *+55 (31) 97560-1545 (TIM)* *"**Nós somos aquilo que fazemos repetidamente. Excelência, então, não é um modo de agir, mas um hábito." Aristóteles*

3 3

Frequência de um registro específico
by Wagner Tassinari 01 Jul '16

01 Jul '16

Olá pessoal, Existe algum comando para verificarmos a frequencia de um registro em específico, por exemplo: Meu banco está dessa forma: var1 var2 .... vark 1 2 0 4 5 6 0 6 0 4 ... Gostaria de saber a frequencia de zeros da var1 obrigadoi ----- 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 ----------------------------------------------------------------------

5 4

Inputando registros
by Wagner Tassinari 30 Jun '16

30 Jun '16

Olá pessoal, tenho o seguinte banco de dados: var1 var2 20 25 40 40 45 40 50 55 50 60 60 Quero criar uma var3 que seja igual a var2 mas onde está em branco impute os valores de var1, ex: var3 25 40 45 40 55 50 60 Estou fazendo da seguinte forma: banco$var1=as.character(banco$var1) banco$var2=as.character(banco$var2) banco$var3= banco$var2 banco$var3[banco$var2 == " "] = banco$var1 Mas não esta dando certo !!!!!!!!!! Obrigado pessoal ----- 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 ----------------------------------------------------------------------

4 3

Dúvidas - Scott-Knott
by Pedro Campelo 28 Jun '16

28 Jun '16

Pessoal, Preciso de uma ajuda na execução desse script abaixo: *#Data.Frame* *> data.frame(A)* * M R* *1 1 0.150* *2 1 0.117* *3 1 0.125* *4 2 0.312* *5 2 0.288* *6 2 0.262* *7 3 0.318* *8 3 0.315* *9 3 0.299* *10 4 0.395* *11 4 0.312* *12 4 0.365* *#ANOVA One Way* *av<-aov(R ~ M, data.frame(A))* *summary(av)* *#Teste de Scott Knott* *sk <- SK(x=av, y=A$R, model='y~x', which='x', dispersion='se')* *summary(sk)* Quando vou executar a rotina, aparece esta mensagem: *> sk <- SK(x=av, y=A$R, model='y~x', which='x', dispersion='se')* *Error in SK.aov(x = av, y = A$R, model = "y~x", which = "x", dispersion = "se") : * * No factors in the fitted model!* *In addition: Warning message:* *In replications(paste("~", xx), data = mf) : non-factors ignored: M* Como poderia resolver isto? Obrigado *----* *Pedro Henrique Campelo Felix* *Bachelor of Food Engineering (UNI-BH)* *Msc. Chemical Engineering (UFMG)* *PhD Student of Food Science (UFLA)* *+55 (31) 97560-1545 (TIM)* *"**Nós somos aquilo que fazemos repetidamente. Excelência, então, não é um modo de agir, mas um hábito." Aristóteles*

3 5

obter data.frames com mesmas dimensões
by Heloíse Pavanato 28 Jun '16

28 Jun '16

Boa tarde, Eu tenho dois data.frames com número de linhas diferentes (nrow(rep1C) = 73, nrow(rep2C) = 45), e colunas iguais (19). Eu quero que ambos data.frames fiquem com a mesma dimensão. Para isso, rep2C deve ser aumentado com "NA". Porém, quero manter os mesmos índices "area", "season", "month" para ambos data.frames. Alguém tem alguma idéia de que função usar? Segue CMR: #---------------------------------------------------------------------- ## Load data # file: *https://www.datafilehost.com/d/64a6a8d9 <https://www.datafilehost.com/d/64a6a8d9>* load('dataC.inia.rda') ## Building matrix for replicates # Sampling 1 rep1C <- dataC.inia[dataC.inia$replicate == 1,] table(rep1C$area, rep1C$season, rep1C$transect) dim(rep1C) # Sampling 2 rep2 <- dataC.inia[dataC$replicate == 2,] table(rep2$area, rep2$season, rep2$transect) dim(rep2) #---------------------------------------------------------------------- Obrigada pela ajuda! Heloise

3 6

Inversa condicional
by Andre Oliveira 27 Jun '16

27 Jun '16

Bom dia, além da função ginv do pacote do MASS para calcula inversa condicional. Qual outro pacote faz esta mesma função só que outras inversas sem ser a de apenas a inversa condicional de Moore-Penrose? 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

Normalização de dados
by Elias Carvalho 24 Jun '16

24 Jun '16

Pessoal, Eu preciso fazer uma análise para verificar a relação das variáveis com a morte e meus dados são semelhantes a tabela abaixo: *CD_MUN* *MUNIC* *UF* *POPULAÇÃO* *UTI* *PIB* *AUTO* *MOTO* *ONIBUS* *CAMINHÃO* *MORTE* 1 A X 3577 10 18533,94 1549 661 40 360 2 *2* *B* *Y* *11958* *2* *12487,36* *5402* *2234* *185* *1401* *11* 3 C Z 6936 8 9947 2187 1478 60 666 6 Nessa tabela o município da linha 2 tem mais população, mais PIB, mais veículos e mais morte. Isso é óbvio porque a sua população e maior, mas isso não necessariamente quer dizer que esse município tenha mais mortes ou que elas são causadas mais provavelmente pela frota de carros ou motos. Então pensando em evitar viéses, e por isso eu pensei em normalizar os dados dividindo as colunas UTI, PIB, AUTO, MOTO, ONIBUS, CAMINHÃO e MORTE pela coluna da POPULAÇÃO e dividir por 10.000, por exemplo para MORTE da linha 2, ficaria: (11/11958)*1000 = 9,1989. Sessa forma a tabela ficaria assim: *CD_MUN* *MUNIC* *UF* *POPULAÇÃO* *UTI* *PIB* *AUTO* *MOTO* *ONIBUS* *CAMINHÃO* *MORTE* 1 A X 3577 27,9564 51814,2018 4330,4445 1847,9172 111,8256 1006,4300 5,5913 *2* *B* *Y* *11958* *1,6725* *10442,6827* *4517,4778* *1868,2054* *154,7081* *1171,6006* *9,1989* 3 C Z 6936 11,5340 14341,1188 3153,1142 2130,9112 86,5052 960,2076 8,6505 Ou ainda poderia usar alguma ferramento do R que fosse mais confiável ? -- Best regards... 8^) “The mind that is open to new ideas never come back to its original size” *Albert Einstein* -- Obrigado Elias

2 2

Operação entre datas
by Amikobh 24 Jun '16

24 Jun '16

Prezados, Desde já, obrigado pela ajuda. Tendo duas variáveis com datas referentes a determinada pessoa: NASCIMENTO DATA_FATO 26/11/97 23/02/15 20/03/99 28/03/15 14/05/99 08/04/15 06/05/98 08/04/15 03/01/00 08/04/15 19/05/97 08/05/15 01/02/01 17/05/15 28/11/97 03/06/15 10/02/00 03/06/15 Como procedo para criar mais uma coluna que calcula a idade no dia da ocorrência de determinado fato? De tal forma que obtenha: NASCIMENTO DATA_FATO IDADE 26/11/97 23/02/15 17 20/03/99 28/03/15 16 14/05/99 08/04/15 15 06/05/98 08/04/15 16 03/01/00 08/04/15 15 19/05/97 08/05/15 17 01/02/01 17/05/15 14 28/11/97 03/06/15 17 10/02/00 03/06/15 15 Valeu, gente!

4 7

Modelagem geoestatística para pixels de imagens de satelite
by Wagner Wolff 23 Jun '16

23 Jun '16

Caros membros. Estou precisando melhorar a resolução espacial de um raster de Temperatura da superfície oriundo de imagens Landsat8. A resolução das imagens Landsat8 são de 30 m, gostaria de saber se posso utilizar uma abordagem geoestatística, considerando os pixels como amostras pontuais? E assim, interpolar para afinar o grid de 30 para 5 m. Att Wagner Wolff

3 3

Georreferenciamento pelo R
by Wagner Tassinari 23 Jun '16

23 Jun '16

Saudações pessoal, Será que vocês saberiam me informar se existe algum script do R de georreferenciamento através do google maps ? Obrigado ----- 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 ----------------------------------------------------------------------

6 6

Modelagens de desfechos contínuos - quais alternativas?
by Pedro Emmanuel Alvarenga Americano do Brasil 22 Jun '16

22 Jun '16

Amigos de R, Estou pra iniciar o plano de análise de um projeto em que farei análise preditiva de um desfecho em escala continua. O que eu domino e farei de qualquer jeito é utilizar uma regressão linear simples (pelos minimos quadrados), possivelmente com alguma transformação do desfecho pra normalizar. Eu ja percebi que seria interessante comparar essa abordagem inicial com outros tipos de modelos. A minha pergunta é quais os modelos? POis na minha experiência, arvores, NN, bagging and boosting, random forest aparecem muito para desfechos binarios. Nao me recordo de exemplos pra desfechos contínuos. Eu achei uma função de NN que modela desfecho continuous no pacote AMORE, mas não achei nenhuma função com exemplo de de desfecho contínuo em outros tipos de modelo ainda. Se houvesse a indicação dos respectivos pacotes/funções pra eu estudar, ficaria bastante agradecido. Abraço forte a todos, Pedro Brasil

1 0

(sem assunto)
by Elias Carvalho 22 Jun '16

22 Jun '16

*Pessoal estou fazendo um grafo e usei os seguintes comandos:* graph.structure<- Rgraphviz::layoutGraph(graph.structure, edgeAttrs=list(label=labels)) graph::nodeRenderInfo(graph.structure) *e salvando ele:* png(filename=graph.name, units="in", width=10, height=8, pointsize=12, res=300) graph <- Rgraphviz::renderGraph(graph.structure) dev.off() *Quando abro o grafo vejo o problema abaixo onde o label extrapola o tamanho do nó* Alguem pode me dar uma dica ? -- Best regards... 8^) “The mind that is open to new ideas never come back to its original size” *Albert Einstein* -- Obrigado Elias

2 1

Aumentar nó do grafo
by Elias Carvalho 22 Jun '16

22 Jun '16

*Pessoal estou fazendo um grafo e usei os seguintes comandos:* graph.structure<- Rgraphviz::layoutGraph(graph.structure, edgeAttrs=list(label=labels)) graph::nodeRenderInfo(graph.structure) *e salvando ele:* png(filename=graph.name, units="in", width=10, height=8, pointsize=12, res=300) graph <- Rgraphviz::renderGraph(graph.structure) dev.off() *Quando abro o grafo vejo o problema abaixo onde o label extrapola o tamanho do nó* Também tentei fazer assim: Alguem pode me dar uma dica ? -- Best regards... 8^) “The mind that is open to new ideas never come back to its original size” *Albert Einstein* -- Obrigado Elias

2 1

Re: [R-br] regressão polinomial: gráfico e valor máximo (resolvido)
by Maurício Lordêlo 21 Jun '16

21 Jun '16

Muito obrigado a todos. Tudo esclarecido. Maurício Em 21 de junho de 2016 13:42, Éder Comunello <comunello.eder(a)gmail.com> escreveu: > Maurício e colegas, boa tarde! > > Mais uma opção de código para o problema. > [image: Imagem inline 1] > > # <code r> > # DIC 4 rep x 7 doses de gesso (trat): 0, 50, 100, 150, 200, 250, 300 > kg/ha > # Peso de 1.000 sementes (peso, gramas) > peso <- c(134.8, 139.7, 147.6, 132.3, 161.7, 157.7, 150.3, 144.7, > 160.7, 172.7, 163.4, 161.3, 169.8, 168.2, 160.7, 161.0, > 165.7, 160.0, 158.2, 151.0, 171.8, 157.3, 150.4, 160.4, > 154.5, 160.4, 148.8, 154.0) > trat <- rep(seq(0,300,50), each=4) > dados <- data.frame(peso, trat=as.factor(trat)) > > # Método dos Polinômios Ortogonais - Banzato & Kronka, p. 204 > contrasts(dados$trat) = contr.poly(levels(dados$trat)); > contrasts(dados$trat) > > fit1 = aov(peso ~ trat, dados) > summary(fit1) # anova(fit1) > summary(fit1, split = list(trat = list(LIN = 1, QUA = 2, CUB=3, > DESVIOS=4:6))) > # Df Sum Sq Mean Sq F value Pr(>F) > # trat 6 1941.8 323.6 7.668 0.000188 *** > # trat: LIN 1 423.2 423.2 10.026 0.004653 ** > # trat: QUA 1 1285.8 1285.8 30.465 1.78e-05 *** > # trat: CUB 1 155.0 155.0 3.673 0.068997 . > # trat: DESVIOS 3 77.8 25.9 0.614 0.613269 > # Residuals 21 886.3 42.2 > > model.tables(fit1, "means") > trat.m <- tapply(peso, trat, mean); trat.m > > # Uso do ajuste quadrático > fit2 <- lm(peso ~ I(trat)+I(trat^2)) > summary(fit2) > coef(fit2) > eq <- paste0("Y = ", round(coef(fit2)[1], 4), " + ", > round(coef(fit2)[2], 4), "X - ", abs(round(coef(fit2)[3], > 6)), "X²"); eq > > # Coeficiente de determinação > summary(fit2)$r.sq > cor(fitted(fit2), peso)**2 > sumsq <- summary(fit1, split = list(trat = list(LIN = 1, QUA = 2, CUB=3, > RES=4:6)))[[1]][2] > rsq <- sum(sumsq[2:3,])/sumsq[1,]; rsq # Banzato & Kronka, p. 209 > eq.r <- paste0("R² = ", round(rsq, 4)); eq.r > > # Predições > pred <- data.frame(trat = seq(0,300,1)) > pred$val <- predict(fit2, pred); head(pred) > > # Gráfico > par(xpd=T) > plot(c(0,300), c(135,170), type="n", xlab="Doses de gesso (kg/ha)", > ylab="Peso de 1.000 sementes (g)", bty="n", xaxs="i", yaxs="i", las=1) > points(as.numeric(names(trat.m)), trat.m, pch=18, cex=1) > lines(pred$trat, pred$val, lwd=2, lty=1, col=1) > text(300, 129, "X") > text(-27, 170, "Y") > text(150, 153, eq, adj=c(0.5,0.5)) > text(150, 151, eq.r, adj=c(0.5,0.5)) > > # Max - Método 1 > ymax <- max(pred$val); ymax # 164.7042 gramas > xmax <- pred[which.max(pred$val), "trat"]; xmax # com 175 Kg/ha de gesso > segments(xmax, 135, xmax, ymax, lty=3, col=2) > segments(0, ymax, xmax, ymax, lty=3, col=2) > points(xmax, ymax, pch=20, col=2) > > # Max - Método 2 > # ax^2 + bx + c > coefs <- data.frame(t(coef(fit2))); names(coefs) <- letters[3:1]; coefs > with(coefs, -b/(2*a)) # 174.8403 > with(coefs, -(b^2 - 4*a*c)/(4*a)) # 164.7043 > > # Max - Método 3 > fun <- function(x) 140.7839 + 0.2736*x - 0.000782*x**2 > fun(175) > # optimize(fun, interval=c(0, 300), maximum=F) > optimize(fun, interval=c(0, 300), maximum=T) > # $maximum > # [1] 174.9361 > # > # $objective > # [1] 164.7152 > > # Max - Derivadas > > D1 <- D(expression(140.7839 + 0.2736*x - 0.000782*x**2), "x"); D1 > # 0.2736 - 0.000782 * (2 * x) > # 0.2736 - 0.001564x > > D2 <- D(D1, "x"); D2 > # -(0.000782 * 2) > # -0.001564 > > # 0.2736 - 0.001564x = 0 > solve(0.001564, 0.2736) # Máximo > # [1] 174.9361 > fun(174.9361) # 164.7152 > > deriva <- deriv(~140.7839 + 0.2736*x - 0.000782*x**2,"x"); deriva > x <- 175; eval(deriva) > x <- 174.8403; eval(deriva) > x <- 174.9361; eval(deriva) > # </code> > > > ================================================ > Éder Comunello > Researcher at Brazilian Agricultural Research Corporation (Embrapa) > DSc in Agricultural Systems Engineering (USP/Esalq) > MSc in Environ. Sciences (UEM), Agronomist (UEM) > --- > Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| > ================================================ > GEO, -22.2752, -54.8182, 408m > UTC-04:00 / DST: UTC-03:00 > > > > > Em 20 de junho de 2016 19:34, Maurício Lordêlo <r-br(a)listas.c3sl.ufpr.br> > escreveu: > >> Saudações a todos desta lista, >> Estou tentando reproduzir um exemplo de regressão polinomial do livro do >> Banzatto e Kronka. >> Minhas dúvidas são: >> 1. Como reproduzir exatamente o gráfico ilustrado no exemplo? O gráfico >> que consegui fazer está parecido porém não está igual. >> O gráfico apresentado no livro encontra-se neste link: >> https://www.dropbox.com/s/qyv7ofwryegcu7m/figura.jpg?dl=0 >> 2. Como encontrar o valor de X que anula a derivada primeira? E depois >> como encontrar o máximo da função? >> Grato, >> Maurício >> >> Segue o CMR: >> >> ################################################################################ >> #Experimento inteiramente casualizado com 4 repetições para estudar os >> efeitos de >> #7 doses de gesso (tratamentos): >> #0, 50, 100, 150, 200, 250 e 300 kg/ha sobre diversas características do >> feijoeiro >> #Para a característica "peso de 1000 sementes", tem-se os resultados >> obtidos em gramas >> peso = c(134.8, 139.7, 147.6, 132.3, >> 161.7, 157.7, 150.3, 144.7, >> 160.7, 172.7, 163.4, 161.3, >> 169.8, 168.2, 160.7, 161.0, >> 165.7, 160.0, 158.2, 151.0, >> 171.8, 157.3, 150.4, 160.4, >> 154.5, 160.4, 148.8, 154.0) >> tratamento = factor(rep(seq(0,300,50), each=4)) >> dados = data.frame(peso, tratamento) >> attach(dados) >> >> niveis.tratamento = seq(0,300,50) >> contrasts(dados$tratamento) = contr.poly(niveis.tratamento) >> contrasts(dados$tratamento) >> modelo.aov = aov(peso ~ tratamento, dados) >> summary(modelo.aov) >> summary(modelo.aov, split = list(tratamento = list(L = 1, Q = 2, C=3, >> Dev=c(4,5,6)))) >> >> plot(c(0,300), c(130,175), type="n", xlab="Doses de gesso (kg/ha)", >> ylab="Peso (g) de 1000 sementes") >> d = rep(seq(0,300,50), each=4) >> med_t = rep(tapply(peso,tratamento,mean),each=4) >> points(d,med_t) >> x = seq(0,300,0.1) >> lines(x, (predict(lm(peso ~ I(d)+I(d^2)+I(d^3)),data.frame(d=x))),lty=2, >> col="blue") >> >> >> >> _______________________________________________ >> 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. >> > >

1 0

regressão polinomial: gráfico e valor máximo
by Maurício Lordêlo 21 Jun '16

21 Jun '16

Saudações a todos desta lista, Estou tentando reproduzir um exemplo de regressão polinomial do livro do Banzatto e Kronka. Minhas dúvidas são: 1. Como reproduzir exatamente o gráfico ilustrado no exemplo? O gráfico que consegui fazer está parecido porém não está igual. O gráfico apresentado no livro encontra-se neste link: https://www.dropbox.com/s/qyv7ofwryegcu7m/figura.jpg?dl=0 2. Como encontrar o valor de X que anula a derivada primeira? E depois como encontrar o máximo da função? Grato, Maurício Segue o CMR: ################################################################################ #Experimento inteiramente casualizado com 4 repetições para estudar os efeitos de #7 doses de gesso (tratamentos): #0, 50, 100, 150, 200, 250 e 300 kg/ha sobre diversas características do feijoeiro #Para a característica "peso de 1000 sementes", tem-se os resultados obtidos em gramas peso = c(134.8, 139.7, 147.6, 132.3, 161.7, 157.7, 150.3, 144.7, 160.7, 172.7, 163.4, 161.3, 169.8, 168.2, 160.7, 161.0, 165.7, 160.0, 158.2, 151.0, 171.8, 157.3, 150.4, 160.4, 154.5, 160.4, 148.8, 154.0) tratamento = factor(rep(seq(0,300,50), each=4)) dados = data.frame(peso, tratamento) attach(dados) niveis.tratamento = seq(0,300,50) contrasts(dados$tratamento) = contr.poly(niveis.tratamento) contrasts(dados$tratamento) modelo.aov = aov(peso ~ tratamento, dados) summary(modelo.aov) summary(modelo.aov, split = list(tratamento = list(L = 1, Q = 2, C=3, Dev=c(4,5,6)))) plot(c(0,300), c(130,175), type="n", xlab="Doses de gesso (kg/ha)", ylab="Peso (g) de 1000 sementes") d = rep(seq(0,300,50), each=4) med_t = rep(tapply(peso,tratamento,mean),each=4) points(d,med_t) x = seq(0,300,0.1) lines(x, (predict(lm(peso ~ I(d)+I(d^2)+I(d^3)),data.frame(d=x))),lty=2, col="blue")

4 3

Análise de dados medidas repetidas no tempo
by Luiz Roberto Martins Pinto 17 Jun '16

17 Jun '16

Eu enviei para várias pessoas um script R para análise de dados de medidas repetidas no tempo, elaborado pelo prof. Ivan Bezerra (UESC). Solicito o favor de informar: 1) Se o script foi relevante para as análises. 2) Se podem informar o script que efetivamente foi utilizado. 3) Novas sugestões de análises. Aos que não solicitaram o script, mas que podem colaborar, agradeço se o fizerem. Abraços, 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

2 1

Manipular dados em grupos e mesmo data frame
by Edimeire Alexandra Pinto 16 Jun '16

16 Jun '16

OI pessoal. Estou com uma base de dados imensa que quase estoura a memória do R e por isso, parafacilitar o entendimento de minha pergunta, suponha que a base de dados seja amtcars que existe no próprio R. Segue um pedaço dela abaixo. > mtcars mpg cyl disp hp drat wt qsec vs am gear carb Mazda RX4 21.0 6 160.0 110 3.90 2.62016.46 0 1 4 4 Mazda RX4 Wag 21.0 6 160.0 110 3.90 2.87517.02 0 1 4 4 Datsun 710 22.8 4 108.0 93 3.85 2.320 18.61 1 1 4 1 Hornet 4 Drive 21.4 6 258.0 110 3.08 3.21519.44 1 0 3 1 Hornet Sportabout 18.7 8 360.0 175 3.15 3.44017.02 0 0 3 2 Valiant 18.1 6 225.0 105 2.76 3.46020.22 1 0 3 1 Duster 360 14.3 8 360.0 245 3.21 3.57015.84 0 0 3 4 Como faço para agrupar a variável gear em grupos 3, depoisde 4, etc, conforme a ordenação das variáveis cyl ( em 4,6,8,..) e VS ( em 0,1),ou seja, criar grupo da variável gear=3 onde a variáveis cyl e VS estejamordenadas, depois criar grupo da variável gear=4 com variáveis cyl e VS ordenadas,e assim por adiante....quero todos em uma mesma base de dados, ou seja, mesmodata frame. Peço isso a vocês, pois meu objetivo final é pegar variáveissempre que na linha i se tenha VS=1 e na linha i+2 cyl=2 com gear=4, e assimsucessivamente, ou seja, para cada grupo obtido para gear, intercruzar, VS=1 nalinha i e cyl=2 na linha i+2.Tentei usar o pacote dylpr e usar depois os comandos:selecao<-arrange(mtcars,cyl,VS,mpg, disp, hp,drat,wt,qsec,am,gear,carb) g<-group_by(selecao,cyl,vs,gear) conj<-filter(select(g,mpg,vs,cyl,gear),(vs==1 & lead(cyl,2)==4) | gear==4) Mas ele não separou por grupos, pois ao usar lead (lags parafrente) ele não fica restrito a cada grupo da variável formado com a gear, ouseja, o R desconsidera os grupos formados e não limita os “lags” para frente emapenas os grupos formados. Gente, se alguém puder ajudar, agradeço muito.

2 1

Teste em SPARC
by Daniela Petruzalek 13 Jun '16

13 Jun '16

Bom dia pessoal! Alguém aqui tem acesso a um servidor SPARC e poderia fazer um pequeno teste para mim? Eu preciso validar o meu código antes de submeter uma nova release para o CRAN. O objetivo era corrigir um bug que apareceu no Solaris. Eu já testei no Solaris x86 mas gostaria de validar no Solaris SPARC antes de submeter de novo. O código que eu quero validar eu já submeti ao repo no GitHub: https://github.com/danicat/read.dbc/ Obrigada! Daniela

1 0

Windows Server 2008 (64-bit)
by Daniel Marcelino 09 Jun '16

09 Jun '16

Olá, será que alguém da lista tem acesso a um Windows Server 2008 (64-bit)? Eu mantenho um pacote que está apresentando um erro de instalação nesse tipo de máquina, mas eu não consigo explorar o problema porque não tem acesso a uma máquina semelhante. Eu acredito que o problema esteja relacionado a versão GCC, mas é apenas uma intuição. Se alguém tiver a máquina e puder tentar instalar e me mandar o log de erro eu serei muito grato. https://cran.r-project.org/web/packages/SciencesPo/index.html Daniel

3 4

Re: [R-br] Estatística circular - gráfico
by Cassiano 09 Jun '16

09 Jun '16

Segue: Dados AMBIENTE mess angle DISTANCIA ABUND RIQUEZA NASC 1 0 0 2 2 NASC 2 30 0 3 1 NASC 3 60 0 4 3 NASC 4 90 0 8 5 NASC 5 120 0 0 0 NASC 6 150 0 0 0 NASC 7 180 0 31 1 NASC 8 210 0 10 3 NASC 9 240 0 18 3 NASC 10 270 0 7 5 NASC 11 300 0 2 2 NASC 12 330 0 1 1 ############################################################### ################### início da análise ############################## ############################################################### library(circular) nascente=read.table("nasc.csv",h=T) nascente nascente0=nascente[which(nascente$DISTANCIA==0),] nascente0 meses<-c("Mai","Jun","Jul","Ago","Set","Out","Nov","Dez","Jan","Fev","Mar","Abr") #convertendo os dados para circular, usando radianos nasc.c0<-circular(rep(rad(nascente0$angle),nascente0$ABUND),units=c("radians"),rotation="clock") #%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Angulo médio + Desvio padrão angular%%%%%%%%%%%%%%%%%%%% # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% media=mean(nasc.c0) sqrt(-2*log(rho.circular(nasc.c0))) #%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% #### Gráfico ############ par(mfrow=c(1,3)) par(mar=c(0,0,2,0)) plot(nasc.c0,axes=F,main="0") axis.circular(at=circular(rad(nascente0$angle)),labels=meses,units=c("radians"), rotation="clock") rose.diag(nasc.c0,add=T,axes=F) arrows.circular(median(nasc.c0),y=rho.circular(nasc.c0),lwd=4,lty=1,angle=90,len=0.0,col="green") #aqui eu inseri a mediana arrows.circular(mean(nasc.c0),y=rho.circular(nasc.c0),lwd=4,lty=1,angle=90,len=0.0,col="red") #aqui eu inseri a média #agora preciso inserir os desvios, ou erros padrão. ########################### FIM ###################################### Obrigado Em 8 de junho de 2016 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > 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. Re: Estatística circular - gráfico (Éder Comunello) > 2. Erro modelo - glm (Andre Oliveira) > 3. Re: Erro modelo - glm (Marcus Nunes) > 4. Re: Erro modelo - glm (Leonardo Ferreira Fontenelle) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 7 Jun 2016 14:39:28 -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] Estatística circular - gráfico > Message-ID: > <CABmC8gn-f8sdoHs0-Hc=zAGTmd=iYjMHizFe= > X_HK6JYa4Vobw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Cassiano, boa tarde! > > Talvez seja o caso de postar seu código com um dataset de exemplo. > > ================================================ > Éder Comunello > Researcher at Brazilian Agricultural Research Corporation (Embrapa) > DSc in Agricultural Systems Engineering (USP/Esalq) > MSc in Environ. Sciences (UEM), Agronomist (UEM) > --- > Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>| > ================================================ > GEO, -22.2752, -54.8182, 408m > UTC-04:00 / DST: UTC-03:00 > > > > > Em 7 de junho de 2016 10:11, Cassiano <cassianosr(a)gmail.com> escreveu: > > > Obrigado pela ajuda. > > Realmente o ggplot2 produz gráficos mais interessantes. Mas ainda não > > achei a opção para plotar a média e os desvios angulares. > > Se puderem me dar mais instruções, ficarei agradecido. > > > > Abraços > > Cassiano > > > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160607/2bc39a6a/atta… > > > > ------------------------------ > > Message: 2 > Date: Tue, 7 Jun 2016 23:26:58 +0000 (UTC) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br> > To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Erro modelo - glm > Message-ID: > <1694617424.92644.1465342018577.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > boa noite, alguém poderia ajudar descobrir onde está o erro? > > remove(list=ls())diurno=c(68,36,40,30) > noturno=c(46,26,55,22) > quadra= c("A", "B", "G","D") > dados=data.frame(diurno,noturno,quadra) > fit = glm(cbind(diurno, noturno) ~ factor(quadra), family = poisson, > data=dados) > Error in x[good, , drop = FALSE] : > (subscript) subscrito lógico muito longo > 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 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160607/7d2e5cef/atta… > > > > ------------------------------ > > Message: 3 > Date: Tue, 7 Jun 2016 22:37:17 -0300 > From: Marcus Nunes <marcus.nunes(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Andre Oliveira > <andreolsouza(a)yahoo.com.br> > Subject: Re: [R-br] Erro modelo - glm > Message-ID: > < > CA+QGQvvGr8Rg5w+mO4tNmRiG1du+nF0Z83XCZUrDfYyhgZdTww(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Ajuste os modelos individualmente: > > diurno=c(68,36,40,30) > noturno=c(46,26,55,22) > quadra= c("A", "B", "G","D") > dados=data.frame(diurno,noturno,quadra) > fit1 = glm(diurno ~ factor(quadra), family = poisson, data=dados) > fit2 = glm(noturno ~ factor(quadra), family = poisson, data=dados) > > Isto deve resolver o teu problema. > > > > 2016-06-07 20:26 GMT-03:00 Andre Oliveira <andreolsouza(a)yahoo.com.br>: > > > boa noite, alguém poderia ajudar descobrir onde está o erro? > > > > remove(list=ls()) > > diurno=c(68,36,40,30) > > noturno=c(46,26,55,22) > > quadra= c("A", "B", "G","D") > > dados=data.frame(diurno,noturno,quadra) > > fit = glm(cbind(diurno, noturno) ~ factor(quadra), family = poisson, > > data=dados) > > > > *Error in x[good, , drop = FALSE] : * > > * (subscript) subscrito lógico muito longo* > > > > 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. > > > > > > -- > Marcus Nunes > http://marcusnunes.me/ > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160607/a1cd75b1/atta… > > > > ------------------------------ > > Message: 4 > Date: Wed, 08 Jun 2016 11:16:04 -0300 > From: Leonardo Ferreira Fontenelle <leonardof(a)leonardof.med.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Erro modelo - glm > Message-ID: > <1465395364.4116631.631632273.0B1F3BA2(a)webmail.messagingengine.com > > > Content-Type: text/plain; charset="utf-8" > > Em Ter 7 jun. 2016, às 20:26, Andre Oliveira escreveu: > > boa noite, alguém poderia ajudar descobrir onde está o erro? > > > > remove(list=ls()) > > diurno=c(68,36,40,30) > > noturno=c(46,26,55,22) > > quadra= c("A", "B", "G","D") > > dados=data.frame(diurno,noturno,quadra) > > fit = glm(cbind(diurno, noturno) ~ factor(quadra), family = poisson, > > data=dados) > > *Error in x[good, , drop = FALSE] : > * > > * (subscript) subscrito lógico muito longo* > > Andre, conforme está escrito na documentação do glm, a resposta pode ser > fornecida como duas colunas (sucessos e fracassos) *apenas para > regressões binomiais e quasi-binomiais*; não para Poisson. > > Supondo que seus dados sejam o resumo de 68 + 46 = 114 observações para > a quadra A, 36 + 26 = 62 observações para a quadra B, e por aí em > diante; e que você queira calcular uma razão de prevalência / risco / > outra proporção entre as quadras, o jeito mais óbvio de fornecer seus > dados para a função é em formato longo: > > dados <- data.frame(horariodiurno = c(1, 0, 1, 0, 1, 0, 1, 0), quadra = > c("A", "A", "B", "B", "C", "C", "D", "D")) > dados <- dados[rep(1:nrow(dados), c(68, 46, 36, 26, 40, 55, 30, 22)), ] > fit = glm(horariodiurno ~ quadra, family = quasipoisson, data = dados) > > Repare também que, se você quer trabalhar com razão de proporções, a > distribuição Poisson em princípio não se aplica, daí a necessidade de > você especificar uma quasipoisson ou então um estimador de variância > consistente com heterocedasticidade (pacote sandwich). > > Por favor, use o código com cuidado porque não estou com meu computador > aqui para testá-lo. > > Atenciosamente, > > Leonardo Ferreira Fontenelle[1] > > Links: > > 1. http://lattes.cnpq.br/9234772336296638 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160608/7af199f3/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 66, assunto 7 > **************************************** > -- ======================================= Cassiano S. Rosa CRBio: 104117/04-D Universidade Federal do Triângulo Mineiro (UFTM), *campus* Iturama-MG Avenida Rio Paranaíba, 1241, Centro, 38280-000 Tel: (34)3415-2512 Lattes <http://lattes.cnpq.br/3030825329622014> Scholar Google <http://scholar.google.com/citations?user=JmELzOMAAAAJ> Academia.edu <https://uftm.academia.edu/CassianoRosa> Researchgate <https://www.researchgate.net/profile/Cassiano_Rosa> =======================================

2 2

I found it!
by david_ufvjm＠yahoo.com.br 09 Jun '16

09 Jun '16

Hello, I know you were looking for that stuff for a long time, so I guess I found it, take a look here <http://synkyxejo.elearning.xyz/xykdyitd> Take care, david_ufvjm(a)yahoo.com.br

1 0

Erro modelo - glm
by Andre Oliveira 08 Jun '16

08 Jun '16

boa noite, alguém poderia ajudar descobrir onde está o erro? remove(list=ls())diurno=c(68,36,40,30) noturno=c(46,26,55,22) quadra= c("A", "B", "G","D") dados=data.frame(diurno,noturno,quadra) fit = glm(cbind(diurno, noturno) ~ factor(quadra), family = poisson, data=dados) Error in x[good, , drop = FALSE] : (subscript) subscrito lógico muito longo 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 2

Estatística circular - gráfico
by Cassiano 07 Jun '16

07 Jun '16

Obrigado pela ajuda. Realmente o ggplot2 produz gráficos mais interessantes. Mas ainda não achei a opção para plotar a média e os desvios angulares. Se puderem me dar mais instruções, ficarei agradecido. Abraços Cassiano Em 6 de junho de 2016 16:08, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > 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. Re: Leitura de Arquivo (Anselmo Alves de Sousa) > 2. R e DATASUS (Daniela Petruzalek) > 3. Re: R e DATASUS (Leonardo Ferreira Fontenelle) > 4. Estatística circular - gráfico (Cassiano) > 5. Re: Estatística circular - gráfico (Augusto Ribas) > 6. Re: R e DATASUS (Daniela Petruzalek) > 7. Re: Estatística circular - gráfico (Edimeire Alexandra Pinto) > 8. Re: R e DATASUS (Elias Teixeira Krainski) > 9. Re: R e DATASUS (Julimar Pinto) > 10. Re: R e DATASUS (Leonardo Ferreira Fontenelle) > 11. Re: R e DATASUS (Daniela Petruzalek) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 6 Jun 2016 12:28:41 -0300 > From: Anselmo Alves de Sousa <startermath(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Leitura de Arquivo > Message-ID: > <CABb2rDq9cxbdfTH8WDPjorPRxnDckDxSte4ziy8zJuRPjgF7= > w(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Muito obrigado. > > Não atentei que um simples read.table() resolveria o problema. Fui traído > pelo .TL > > Restou uma dúvida: preciso transformar esse dado numa variável série > temporal da taxa de envio do pacote por segundo. Como posso agregar esses > pacotes por segundo? > []'s > Anselmo > > Em 4 de junho de 2016 02:33, salah <salah3.1416(a)gmail.com> escreveu: > > > sugestão: > > > > dataset = read.csv("C:\\Users\\Anselmo\\Desktop\\trab\\BC-pOct89.TL-1", > > sep="", header = FALSE) > > > > saudações > > > > Em 02/06/2016 00:15, Anselmo Alves de Sousa escreveu: > > > > Boa noite listeiros > > > > Estou com a missão de auxiliar um colega no trabalho escolar. A primeira > > missão é ler o dado. Como begginer no R não estou conseguindo ler o > arquivo > > ftp://ita.ee.lbl.gov/traces/BC-pOct89.TL.Z. Sabemos que tem duas > colunas, > > uma com o timestamp e outra com o tamanho em bytes. > > > > Usando > > > > > dataset<-scan("C:\\Users\\Anselmo\\Desktop\\trab\\BC-pOct89.TL\\BC-pOct89.TL") > > > > no diretório apropriado > > > > setwd("C://Users//Anselmo//Desktop//trab_vinicius//BC-pOct89.TL") > > e str(dataset) > > num [1:2000000] 0.0177 87 0.0368 142 0.0368 ... > > > > Ou seja, obtenho apenas as duas colunas numa só. Como posso ler > > adequadamente esse arquivo? > > > > []'s > > Anselmo > > > > > > _______________________________________________ > > R-br mailing listR-br@listas.c3sl.ufpr.brhttps:// > 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/20160606/e3cb746d/atta… > > > > ------------------------------ > > Message: 2 > Date: Mon, 06 Jun 2016 16:35:19 +0000 > From: Daniela Petruzalek <daniela.petruzalek(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] R e DATASUS > Message-ID: > <CA+=iygsQDOaso54n-RpodxGW= > EkAyYNoAMS_x28SbVoxDyhhxw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Boa tarde pessoal! > > Eu sou nova na lista mas tive a oportunidade de ler algumas discussões pelo > histórico na WEB e resolvi me juntar ao grupo. Por isso, antes de mais > nada, uma rápida introdução: eu sou nova na linguagem R (e em Data Science > de modo geral) mas tenho estudado bastante sobre o tema, visto que este se > tornou um dos principais temas em tecnologia nos últimos dois anos. > > Eu sou especialista em banco de dados relacionais (principalmente Oracle) e > também tenho experiência com linguagens de programação, então embora nova > no R acredito que eu tenha uma bagagem para contribuir com a comunidade. > > Uma das minhas áreas de interesse além da TI é saúde, e eu vi que várias > pessoas já perguntaram nesta lista sobre uma forma de abrir os arquivos do > DATASUS no R nativamente. Eu já tive oportunidade de trabalhar com o TabWin > no passado e não era muito feliz, sempre fui a favor de plataformas mais > flexíveis e abertas. > > Por este motivo eu criei uma package para o R chamada "read.dbc" que > permite abrir os arquivos DBC nativamente, sem depender do ferramental que > o DATASUS disponibiliza. Inclusive, por ser compilada diretamente no R, ela > pode ser utilizada em Linux sem precisar de Wine ou qualquer outra > ferramenta. > > Esta package é completamente open source e pode ser utilizada por qualquer > pessoa que se interessar. Caso a utilizem nos seus trabalhos, porém, só > solicito que o meu nome seja creditado de acordo com a licença AGPL-3. > > A quem interessar, segue o link para o código fonte/download: > https://github.com/danicat/read.dbc > > Att., > > Daniela > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160606/4b9e4f9f/atta… > > > > ------------------------------ > > Message: 3 > Date: Mon, 06 Jun 2016 13:48:58 -0300 > From: Leonardo Ferreira Fontenelle <leonardof(a)leonardof.med.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] R e DATASUS > Message-ID: > <1465231738.3516313.629468385.2970D8EF(a)webmail.messagingengine.com > > > Content-Type: text/plain; charset="utf-8" > > Daniela, tenho um interesse preciso em utilizar bancos de dados do > DATASUS, e agradeço sua iniciativa! > > Vou instalar e testar assim que possível. > > Leonardo Ferreira Fontenelle[1] > > > Em Seg 6 jun. 2016, às 13:35, Daniela Petruzalek escreveu: > > Boa tarde pessoal! > > Eu sou nova na lista mas tive a oportunidade de ler algumas discussões > > pelo histórico na WEB e resolvi me juntar ao grupo. Por isso, antes de > > mais nada, uma rápida introdução: eu sou nova na linguagem R (e em > > Data Science de modo geral) mas tenho estudado bastante sobre o tema, > > visto que este se tornou um dos principais temas em tecnologia nos > > últimos dois anos. > > > > Eu sou especialista em banco de dados relacionais (principalmente > > Oracle) e também tenho experiência com linguagens de programação, > > então embora nova no R acredito que eu tenha uma bagagem para > > contribuir com a comunidade. > > > > Uma das minhas áreas de interesse além da TI é saúde, e eu vi que > > várias pessoas já perguntaram nesta lista sobre uma forma de abrir os > > arquivos do DATASUS no R nativamente. Eu já tive oportunidade de > > trabalhar com o TabWin no passado e não era muito feliz, sempre fui a > > favor de plataformas mais flexíveis e abertas. > > > > Por este motivo eu criei uma package para o R chamada "read.dbc" que > > permite abrir os arquivos DBC nativamente, sem depender do ferramental > > que o DATASUS disponibiliza. Inclusive, por ser compilada diretamente > > no R, ela pode ser utilizada em Linux sem precisar de Wine ou qualquer > > outra ferramenta. > > > > Esta package é completamente open source e pode ser utilizada por > > qualquer pessoa que se interessar. Caso a utilizem nos seus trabalhos, > > porém, só solicito que o meu nome seja creditado de acordo com a > > licença AGPL-3. > > > > A quem interessar, segue o link para o código fonte/download: > > https://github.com/danicat/read.dbc > > > > Att., > > Daniela > > _________________________________________________ > > 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. > > > Links: > > 1. http://lattes.cnpq.br/9234772336296638 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160606/43f0ec03/atta… > > > > ------------------------------ > > Message: 4 > Date: Mon, 6 Jun 2016 13:56:58 -0300 > From: Cassiano <cassianosr(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, r-br-request(a)listas.c3sl.ufpr.br > Subject: [R-br] Estatística circular - gráfico > Message-ID: > < > CAF78VjJG0v81_knO1YUCgDFUQ+B7tbVHdOj13+voD70ouWEB_w(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Prezados, > > Tenho uns dados coletados durante o ano e fiz análise circular. > Quanto à análise, está tranquilo, mas não consigo plotar no gráfico a > variação em torno da média. > Com o "arrows.circular" eu consegui inserir a média no gráfico, mas não > encontrei nada para incluir a variação em torno da média. > > > Desde já agradeço. > cassiano > > > > > -- > ======================================= > Cassiano S. Rosa > CRBio: 104117/04-D > > Universidade Federal do Triângulo Mineiro (UFTM), > *campus* Iturama-MG > Avenida Rio Paranaíba, 1241, Centro, 38280-000 > Tel: (34)3415-2512 > > Lattes <http://lattes.cnpq.br/3030825329622014> > Scholar Google <http://scholar.google.com/citations?user=JmELzOMAAAAJ> > Academia.edu <https://uftm.academia.edu/CassianoRosa> > Researchgate <https://www.researchgate.net/profile/Cassiano_Rosa> > > ======================================= > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160606/aece1d81/atta… > > > > ------------------------------ > > Message: 5 > Date: Mon, 6 Jun 2016 13:41:06 -0400 > From: Augusto Ribas <ribas.aca(a)gmail.com> > To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Estatística circular - gráfico > Message-ID: > < > CACMkfRwLNnMcrKSXrVes_sK+LZi0D6krC3xT68J9x3s+v8q2Ew(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Usando o ggplot2, eu achei relativamente mais fácil fazer gráficos > circulares que com outros pacotes, como o circ e o circular. Basta você > fazer um gráfico normal e depois adicionar + coord_polar() > Segue uns exemplos que eu estava testando, talvez você ache interessante. > > library("ggplot2") > > > dados<-data.frame(preditor=1:20,resposta=runif(20,5,10),grupo=rep(1:4,each=5),linha=runif(20,8,12))head(dados) > ## Gráfico de barras > ggplot(data=dados, aes(x=preditor, y=resposta)) + geom_bar(stat="identity") > ## Separando por grupos > ggplot(data=dados, aes(x=preditor, y=resposta,fill=grupo)) + > geom_bar(stat="identity")##ggplot(data=dados, aes(x=preditor, > y=resposta)) + geom_bar(aes(fill=grupo),stat="identity") > ## Adicionando um contorno preto, e retirando a legenda > ggplot(data=dados, aes(x=preditor, y=resposta,fill=grupo)) + > geom_bar(colour="black",stat="identity") + guides(fill=FALSE) > ##Gráfico de linhas > ggplot(data=dados, aes(x=preditor, y=linha)) + geom_line() + > guides(fill=FALSE) > ggplot(data=dados, aes(x=preditor, y=linha)) + geom_line() + > geom_point() + guides(fill=FALSE) > ##Combinando os dois > ggplot(data=dados, aes(x=preditor, y=resposta,fill=grupo)) + > geom_bar(colour="black",stat="identity") + guides(fill=FALSE) + > geom_line(data=dados, aes(x=preditor, y=linha)) + > geom_point(data=dados, aes(x=preditor, y=linha)) > > ##Transformando em uma figura circular > ggplot(data=dados, aes(x=preditor, y=resposta,fill=grupo)) + > geom_bar(colour="black",stat="identity") + guides(fill=FALSE) + > geom_line(data=dados, aes(x=preditor, y=linha)) + > geom_point(data=dados, aes(x=preditor, y=linha))+coord_polar() > > ggplot(data=dados, aes(x=preditor, y=linha)) + geom_line() + > geom_point() + guides(fill=FALSE) + coord_polar() > ggplot(data=dados, aes(x=preditor, y=linha)) + geom_line() + > geom_point() + guides(fill=FALSE) + coord_polar() + ylim(0,12) > > > > Em 6 de junho de 2016 12:56, Cassiano <cassianosr(a)gmail.com> escreveu: > > > Prezados, > > > > Tenho uns dados coletados durante o ano e fiz análise circular. > > Quanto à análise, está tranquilo, mas não consigo plotar no gráfico a > > variação em torno da média. > > Com o "arrows.circular" eu consegui inserir a média no gráfico, mas não > > encontrei nada para incluir a variação em torno da média. > > > > > > Desde já agradeço. > > cassiano > > > >> > >> -- > > ======================================= > > Cassiano S. Rosa > > CRBio: 104117/04-D > > > > Universidade Federal do Triângulo Mineiro (UFTM), > > *campus* Iturama-MG > > Avenida Rio Paranaíba, 1241, Centro, 38280-000 > > Tel: (34)3415-2512 > > > > Lattes <http://lattes.cnpq.br/3030825329622014> > > Scholar Google <http://scholar.google.com/citations?user=JmELzOMAAAAJ> > > Academia.edu <https://uftm.academia.edu/CassianoRosa> > > Researchgate <https://www.researchgate.net/profile/Cassiano_Rosa> > > > > ======================================= > > > > _______________________________________________ > > 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. > > > > > > -- > Grato > Augusto C. A. Ribas > > Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> > Github: https://github.com/Squiercg > Lattes: http://lattes.cnpq.br/7355685961127056 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160606/027ed7ab/atta… > > > > ------------------------------ > > Message: 6 > Date: Mon, 06 Jun 2016 18:20:30 +0000 > From: Daniela Petruzalek <daniela.petruzalek(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] R e DATASUS > Message-ID: > <CA+=iygssvOQRg78iwsas2sWnBzn= > 03TfE6_98akSX100ZBdqFw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Leonardo, > > Fico feliz em poder ajudar. Se precisar de alguma coisa fique a vontade. > > Aproveitando, talvez este outro projeto meu te interesse: > > https://github.com/danicat/Rdatasus > > Trata-se de uma interface para as bases de dados do DATASUS, para que você > possa carregar os arquivos direto do ftp para o R, sem ter que passar pelo > site. Por enquanto é apenas uma prova de conceito e o uso não é muito > intuitivo, mas já adicionei o suporte ao SIM com todos os seus tipos de > arquivo (DO, DOFET, DOEXT, DOINF e DOMAT). > > Você pode carregar os dados com uma única chamada de função, por exemplo: > > df <- datasus.mortality(2012, states = c("PR","SC")) > > ou > > df <- datasus.mortality(2013, type = "DOINF") > > Eu optei por escrever o projeto em inglês para poder ser utilizado também > internacionalmente, mas o meu objetivo no futuro é disponibilizar o suporte > às duas linguas. > > Este projeto está bem no começo, mas em breve devo trabalhar em melhorar a > usabilidade dele e também disponibilizar as outras bases. > > Att., > > Daniela > > Em seg, 6 de jun de 2016 às 13:49, Leonardo Ferreira Fontenelle < > leonardof(a)leonardof.med.br> escreveu: > > > Daniela, tenho um interesse preciso em utilizar bancos de dados do > > DATASUS, e agradeço sua iniciativa! > > > > Vou instalar e testar assim que possível. > > > > Leonardo Ferreira Fontenelle <http://lattes.cnpq.br/9234772336296638> > > > > > > Em Seg 6 jun. 2016, às 13:35, Daniela Petruzalek escreveu: > > > > Boa tarde pessoal! > > Eu sou nova na lista mas tive a oportunidade de ler algumas discussões > > pelo histórico na WEB e resolvi me juntar ao grupo. Por isso, antes de > mais > > nada, uma rápida introdução: eu sou nova na linguagem R (e em Data > Science > > de modo geral) mas tenho estudado bastante sobre o tema, visto que este > se > > tornou um dos principais temas em tecnologia nos últimos dois anos. > > > > Eu sou especialista em banco de dados relacionais (principalmente Oracle) > > e também tenho experiência com linguagens de programação, então embora > nova > > no R acredito que eu tenha uma bagagem para contribuir com a comunidade. > > > > Uma das minhas áreas de interesse além da TI é saúde, e eu vi que várias > > pessoas já perguntaram nesta lista sobre uma forma de abrir os arquivos > do > > DATASUS no R nativamente. Eu já tive oportunidade de trabalhar com o > TabWin > > no passado e não era muito feliz, sempre fui a favor de plataformas mais > > flexíveis e abertas. > > > > Por este motivo eu criei uma package para o R chamada "read.dbc" que > > permite abrir os arquivos DBC nativamente, sem depender do ferramental > que > > o DATASUS disponibiliza. Inclusive, por ser compilada diretamente no R, > ela > > pode ser utilizada em Linux sem precisar de Wine ou qualquer outra > > ferramenta. > > > > Esta package é completamente open source e pode ser utilizada por > qualquer > > pessoa que se interessar. Caso a utilizem nos seus trabalhos, porém, só > > solicito que o meu nome seja creditado de acordo com a licença AGPL-3. > > > > A quem interessar, segue o link para o código fonte/download: > > https://github.com/danicat/read.dbc > > > > Att., > > Daniela > > > > *_______________________________________________* > > 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/20160606/a458de1f/atta… > > > > ------------------------------ > > Message: 7 > Date: Mon, 6 Jun 2016 18:29:35 +0000 (UTC) > From: Edimeire Alexandra Pinto <economatistica(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Estatística circular - gráfico > Message-ID: > <219253963.1926673.1465237775267.JavaMail.yahoo(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Isso mesmo, usa o ggplot2 > > > http://www.ceb-institute.org/bbs/wp-content/uploads/2011/09/handout_ggplot2… > > Em Segunda-feira, 6 de Junho de 2016 14:41, Augusto Ribas < > ribas.aca(a)gmail.com> escreveu: > > > Usando o ggplot2, eu achei relativamente mais fácil fazer gráficos > circulares que com outros pacotes, como o circ e o circular. Basta você > fazer um gráfico normal e depois adicionar + coord_polar() > Segue uns exemplos que eu estava testando, talvez você ache interessante. > library("ggplot2") > > > dados<-data.frame(preditor=1:20,resposta=runif(20,5,10),grupo=rep(1:4,each=5),linha=runif(20,8,12)) > head(dados) > > ## Gráfico de barras > ggplot(data=dados, aes(x=preditor, y=resposta)) + geom_bar(stat="identity") > > ## Separando por grupos > ggplot(data=dados, aes(x=preditor, y=resposta,fill=grupo)) + > geom_bar(stat="identity") > ##ggplot(data=dados, aes(x=preditor, y=resposta)) + > geom_bar(aes(fill=grupo),stat="identity") > > ## Adicionando um contorno preto, e retirando a legenda > ggplot(data=dados, aes(x=preditor, y=resposta,fill=grupo)) + > geom_bar(colour="black",stat="identity") + guides(fill=FALSE) > > ##Gráfico de linhas > ggplot(data=dados, aes(x=preditor, y=linha)) + geom_line() + > guides(fill=FALSE) > ggplot(data=dados, aes(x=preditor, y=linha)) + geom_line() + geom_point() > + guides(fill=FALSE) > > ##Combinando os dois > ggplot(data=dados, aes(x=preditor, y=resposta,fill=grupo)) + > geom_bar(colour="black",stat="identity") + guides(fill=FALSE) + > geom_line(data=dados, aes(x=preditor, y=linha)) + > geom_point(data=dados, aes(x=preditor, y=linha)) > > > ##Transformando em uma figura circular > ggplot(data=dados, aes(x=preditor, y=resposta,fill=grupo)) + > geom_bar(colour="black",stat="identity") + guides(fill=FALSE) + > geom_line(data=dados, aes(x=preditor, y=linha)) + > geom_point(data=dados, aes(x=preditor, y=linha))+coord_polar() > > ggplot(data=dados, aes(x=preditor, y=linha)) + geom_line() + geom_point() > + guides(fill=FALSE) + coord_polar() > ggplot(data=dados, aes(x=preditor, y=linha)) + geom_line() + geom_point() > + guides(fill=FALSE) + coord_polar() + ylim(0,12) > > Em 6 de junho de 2016 12:56, Cassiano <cassianosr(a)gmail.com> escreveu: > > Prezados, > Tenho uns dados coletados durante o ano e fiz análise circular.Quanto à > análise, está tranquilo, mas não consigo plotar no gráfico a variação em > torno da média.Com o "arrows.circular" eu consegui inserir a média no > gráfico, mas não encontrei nada para incluir a variação em torno da média. > > Desde já agradeço.cassiano > > > -- > =======================================Cassiano S. RosaCRBio: 104117/04-D > Universidade Federal do Triângulo Mineiro (UFTM), campus Iturama-MGAvenida > Rio Paranaíba, 1241, Centro, 38280-000 > Tel: (34)3415-2512 > LattesScholar GoogleAcademia.edu > Researchgate > > ======================================= > > _______________________________________________ > 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. > > > > > -- > Grato > Augusto C. A. Ribas Site Pessoal: http://recologia.com.br/Github: > https://github.com/SquiercgLattes: http://lattes.cnpq.br/7355685961127056 > > _______________________________________________ > 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/20160606/5524dce0/atta… > > > > ------------------------------ > > Message: 8 > Date: Mon, 6 Jun 2016 20:35:27 +0200 > From: Elias Teixeira Krainski <eliaskrainski(a)yahoo.com.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] R e DATASUS > Message-ID: <5755C26F.8090802(a)yahoo.com.br> > Content-Type: text/plain; charset="utf-8"; Format="flowed" > > Parabéns Daniela! Já testei! > > Elias > > On 06/06/16 18:35, Daniela Petruzalek wrote: > > Boa tarde pessoal! > > > > Eu sou nova na lista mas tive a oportunidade de ler algumas discussões > > pelo histórico na WEB e resolvi me juntar ao grupo. Por isso, antes de > > mais nada, uma rápida introdução: eu sou nova na linguagem R (e em > > Data Science de modo geral) mas tenho estudado bastante sobre o tema, > > visto que este se tornou um dos principais temas em tecnologia nos > > últimos dois anos. > > > > Eu sou especialista em banco de dados relacionais (principalmente > > Oracle) e também tenho experiência com linguagens de programação, > > então embora nova no R acredito que eu tenha uma bagagem para > > contribuir com a comunidade. > > > > Uma das minhas áreas de interesse além da TI é saúde, e eu vi que > > várias pessoas já perguntaram nesta lista sobre uma forma de abrir os > > arquivos do DATASUS no R nativamente. Eu já tive oportunidade de > > trabalhar com o TabWin no passado e não era muito feliz, sempre fui a > > favor de plataformas mais flexíveis e abertas. > > > > Por este motivo eu criei uma package para o R chamada "read.dbc" que > > permite abrir os arquivos DBC nativamente, sem depender do ferramental > > que o DATASUS disponibiliza. Inclusive, por ser compilada diretamente > > no R, ela pode ser utilizada em Linux sem precisar de Wine ou qualquer > > outra ferramenta. > > > > Esta package é completamente open source e pode ser utilizada por > > qualquer pessoa que se interessar. Caso a utilizem nos seus trabalhos, > > porém, só solicito que o meu nome seja creditado de acordo com a > > licença AGPL-3. > > > > A quem interessar, segue o link para o código fonte/download: > > https://github.com/danicat/read.dbc > > > > Att., > > > > Daniela > > > > > > _______________________________________________ > > 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/20160606/8ee16b0d/atta… > > > > ------------------------------ > > Message: 9 > Date: Mon, 6 Jun 2016 15:40:32 -0300 > From: Julimar Pinto <julimarsp.r(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] R e DATASUS > Message-ID: > <CAM=r93wpP3Zm8= > ebfr48pDZ+_XBwzGujWeausJ45ZfRrdBXcRw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Ótimo desenvolvimento! Sobretudo especificamente quanto ao "read.dbc", é > uma imensa mão na roda uma outra solução que não o wine. Parabéns!!! > > Julimar > Em 06/06/2016 15:20, "Daniela Petruzalek" <daniela.petruzalek(a)gmail.com> > escreveu: > > > Leonardo, > > > > Fico feliz em poder ajudar. Se precisar de alguma coisa fique a vontade. > > > > Aproveitando, talvez este outro projeto meu te interesse: > > > > https://github.com/danicat/Rdatasus > > > > Trata-se de uma interface para as bases de dados do DATASUS, para que > você > > possa carregar os arquivos direto do ftp para o R, sem ter que passar > pelo > > site. Por enquanto é apenas uma prova de conceito e o uso não é muito > > intuitivo, mas já adicionei o suporte ao SIM com todos os seus tipos de > > arquivo (DO, DOFET, DOEXT, DOINF e DOMAT). > > > > Você pode carregar os dados com uma única chamada de função, por exemplo: > > > > df <- datasus.mortality(2012, states = c("PR","SC")) > > > > ou > > > > df <- datasus.mortality(2013, type = "DOINF") > > > > Eu optei por escrever o projeto em inglês para poder ser utilizado também > > internacionalmente, mas o meu objetivo no futuro é disponibilizar o > suporte > > às duas linguas. > > > > Este projeto está bem no começo, mas em breve devo trabalhar em melhorar > a > > usabilidade dele e também disponibilizar as outras bases. > > > > Att., > > > > Daniela > > > > Em seg, 6 de jun de 2016 às 13:49, Leonardo Ferreira Fontenelle < > > leonardof(a)leonardof.med.br> escreveu: > > > >> Daniela, tenho um interesse preciso em utilizar bancos de dados do > >> DATASUS, e agradeço sua iniciativa! > >> > >> Vou instalar e testar assim que possível. > >> > >> Leonardo Ferreira Fontenelle <http://lattes.cnpq.br/9234772336296638> > >> > >> > >> Em Seg 6 jun. 2016, às 13:35, Daniela Petruzalek escreveu: > >> > >> Boa tarde pessoal! > >> Eu sou nova na lista mas tive a oportunidade de ler algumas discussões > >> pelo histórico na WEB e resolvi me juntar ao grupo. Por isso, antes de > mais > >> nada, uma rápida introdução: eu sou nova na linguagem R (e em Data > Science > >> de modo geral) mas tenho estudado bastante sobre o tema, visto que este > se > >> tornou um dos principais temas em tecnologia nos últimos dois anos. > >> > >> Eu sou especialista em banco de dados relacionais (principalmente > Oracle) > >> e também tenho experiência com linguagens de programação, então embora > nova > >> no R acredito que eu tenha uma bagagem para contribuir com a comunidade. > >> > >> Uma das minhas áreas de interesse além da TI é saúde, e eu vi que várias > >> pessoas já perguntaram nesta lista sobre uma forma de abrir os arquivos > do > >> DATASUS no R nativamente. Eu já tive oportunidade de trabalhar com o > TabWin > >> no passado e não era muito feliz, sempre fui a favor de plataformas mais > >> flexíveis e abertas. > >> > >> Por este motivo eu criei uma package para o R chamada "read.dbc" que > >> permite abrir os arquivos DBC nativamente, sem depender do ferramental > que > >> o DATASUS disponibiliza. Inclusive, por ser compilada diretamente no R, > ela > >> pode ser utilizada em Linux sem precisar de Wine ou qualquer outra > >> ferramenta. > >> > >> Esta package é completamente open source e pode ser utilizada por > >> qualquer pessoa que se interessar. Caso a utilizem nos seus trabalhos, > >> porém, só solicito que o meu nome seja creditado de acordo com a licença > >> AGPL-3. > >> > >> A quem interessar, segue o link para o código fonte/download: > >> https://github.com/danicat/read.dbc > >> > >> Att., > >> Daniela > >> > >> *_______________________________________________* > >> 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. > > > > > > _______________________________________________ > > 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/20160606/927af098/atta… > > > > ------------------------------ > > Message: 10 > Date: Mon, 06 Jun 2016 15:50:17 -0300 > From: Leonardo Ferreira Fontenelle <leonardof(a)leonardof.med.br> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] R e DATASUS > Message-ID: > <1465239017.789326.629595105.21E3703E(a)webmail.messagingengine.com> > Content-Type: text/plain; charset="utf-8" > > Obrigado novamente, vou testar também. > > Leonardo Ferreira Fontenelle[1] > > > Em Seg 6 jun. 2016, às 15:20, Daniela Petruzalek escreveu: > > Leonardo, > > Fico feliz em poder ajudar. Se precisar de alguma coisa fique a > > vontade. > > > > Aproveitando, talvez este outro projeto meu te interesse: > > > > https://github.com/danicat/Rdatasus > > > > Trata-se de uma interface para as bases de dados do DATASUS, para que > > você possa carregar os arquivos direto do ftp para o R, sem ter que > > passar pelo site. Por enquanto é apenas uma prova de conceito e o uso > > não é muito intuitivo, mas já adicionei o suporte ao SIM com todos os > > seus tipos de arquivo (DO, DOFET, DOEXT, DOINF e DOMAT). > > > > Você pode carregar os dados com uma única chamada de função, por > > exemplo: > > df <- datasus.mortality(2012, states = c("PR","SC")) > > ou > > df <- datasus.mortality(2013, type = "DOINF") > > > > Eu optei por escrever o projeto em inglês para poder ser utilizado > > também internacionalmente, mas o meu objetivo no futuro é > > disponibilizar o suporte às duas linguas. > > > > Este projeto está bem no começo, mas em breve devo trabalhar em > > melhorar a usabilidade dele e também disponibilizar as outras bases. > > > > Att., > > > > Daniela > > > > Em seg, 6 de jun de 2016 às 13:49, Leonardo Ferreira Fontenelle > > <leonardof(a)leonardof.med.br> escreveu: > >> __ > >> Daniela, tenho um interesse preciso em utilizar bancos de dados do > >> DATASUS, e agradeço sua iniciativa! > >> > >> Vou instalar e testar assim que possível. > >> > >> Leonardo Ferreira Fontenelle[2] > >> > >> > >> Em Seg 6 jun. 2016, às 13:35, Daniela Petruzalek escreveu: > >>> Boa tarde pessoal! > >>> Eu sou nova na lista mas tive a oportunidade de ler algumas > >>> discussões pelo histórico na WEB e resolvi me juntar ao grupo. Por > >>> isso, antes de mais nada, uma rápida introdução: eu sou nova na > >>> linguagem R (e em Data Science de modo geral) mas tenho estudado > >>> bastante sobre o tema, visto que este se tornou um dos principais > >>> temas em tecnologia nos últimos dois anos. > >>> > >>> Eu sou especialista em banco de dados relacionais (principalmente > >>> Oracle) e também tenho experiência com linguagens de programação, > >>> então embora nova no R acredito que eu tenha uma bagagem para > >>> contribuir com a comunidade. > >>> > >>> Uma das minhas áreas de interesse além da TI é saúde, e eu vi que > >>> várias pessoas já perguntaram nesta lista sobre uma forma de abrir > >>> os arquivos do DATASUS no R nativamente. Eu já tive oportunidade de > >>> trabalhar com o TabWin no passado e não era muito feliz, sempre fui > >>> a favor de plataformas mais flexíveis e abertas. > >>> > >>> Por este motivo eu criei uma package para o R chamada "read.dbc" que > >>> permite abrir os arquivos DBC nativamente, sem depender do > >>> ferramental que o DATASUS disponibiliza. Inclusive, por ser > >>> compilada diretamente no R, ela pode ser utilizada em Linux sem > >>> precisar de Wine ou qualquer outra ferramenta. > >>> > >>> Esta package é completamente open source e pode ser utilizada por > >>> qualquer pessoa que se interessar. Caso a utilizem nos seus > >>> trabalhos, porém, só solicito que o meu nome seja creditado de > >>> acordo com a licença AGPL-3. > >>> > >>> A quem interessar, segue o link para o código fonte/download: > >>> https://github.com/danicat/read.dbc > >>> > >>> Att., > >>> Daniela > >>> _________________________________________________ > >>> 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. > > _________________________________________________ > > 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. > > > Links: > > 1. http://lattes.cnpq.br/9234772336296638 > 2. http://lattes.cnpq.br/9234772336296638 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160606/94f8e0bf/atta… > > > > ------------------------------ > > Message: 11 > Date: Mon, 06 Jun 2016 19:07:49 +0000 > From: Daniela Petruzalek <daniela.petruzalek(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] R e DATASUS > Message-ID: > <CA+=iygvqs-Lgu9f7M8OMCY= > xEiGtpoB_it4TGqxtxT7C8HtQvg(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Obrigada, Julimar! :) > > Eu estou tentando subir este projeto (read.dbc) na CRAN para ficar > acessível para todos os usuários do R. Se tudo der certo logo vamos ter a > package disponível com um simples install.packages, simplificando a vida de > todos os usuários de linux :) > > Atualmente, está quase assim, mas ainda tem a dependência da devtools para > usar o install_github. > > Att., > > Dani > > Em seg, 6 de jun de 2016 às 15:40, Julimar Pinto <julimarsp.r(a)gmail.com> > escreveu: > > > Ótimo desenvolvimento! Sobretudo especificamente quanto ao "read.dbc", é > > uma imensa mão na roda uma outra solução que não o wine. Parabéns!!! > > > > Julimar > > Em 06/06/2016 15:20, "Daniela Petruzalek" <daniela.petruzalek(a)gmail.com> > > escreveu: > > > >> Leonardo, > >> > >> Fico feliz em poder ajudar. Se precisar de alguma coisa fique a vontade. > >> > >> Aproveitando, talvez este outro projeto meu te interesse: > >> > >> https://github.com/danicat/Rdatasus > >> > >> Trata-se de uma interface para as bases de dados do DATASUS, para que > >> você possa carregar os arquivos direto do ftp para o R, sem ter que > passar > >> pelo site. Por enquanto é apenas uma prova de conceito e o uso não é > muito > >> intuitivo, mas já adicionei o suporte ao SIM com todos os seus tipos de > >> arquivo (DO, DOFET, DOEXT, DOINF e DOMAT). > >> > >> Você pode carregar os dados com uma única chamada de função, por > exemplo: > >> > >> df <- datasus.mortality(2012, states = c("PR","SC")) > >> > >> ou > >> > >> df <- datasus.mortality(2013, type = "DOINF") > >> > >> Eu optei por escrever o projeto em inglês para poder ser utilizado > também > >> internacionalmente, mas o meu objetivo no futuro é disponibilizar o > suporte > >> às duas linguas. > >> > >> Este projeto está bem no começo, mas em breve devo trabalhar em melhorar > >> a usabilidade dele e também disponibilizar as outras bases. > >> > >> Att., > >> > >> Daniela > >> > >> Em seg, 6 de jun de 2016 às 13:49, Leonardo Ferreira Fontenelle < > >> leonardof(a)leonardof.med.br> escreveu: > >> > >>> Daniela, tenho um interesse preciso em utilizar bancos de dados do > >>> DATASUS, e agradeço sua iniciativa! > >>> > >>> Vou instalar e testar assim que possível. > >>> > >>> Leonardo Ferreira Fontenelle <http://lattes.cnpq.br/9234772336296638> > >>> > >>> > >>> Em Seg 6 jun. 2016, às 13:35, Daniela Petruzalek escreveu: > >>> > >>> Boa tarde pessoal! > >>> Eu sou nova na lista mas tive a oportunidade de ler algumas discussões > >>> pelo histórico na WEB e resolvi me juntar ao grupo. Por isso, antes de > mais > >>> nada, uma rápida introdução: eu sou nova na linguagem R (e em Data > Science > >>> de modo geral) mas tenho estudado bastante sobre o tema, visto que > este se > >>> tornou um dos principais temas em tecnologia nos últimos dois anos. > >>> > >>> Eu sou especialista em banco de dados relacionais (principalmente > >>> Oracle) e também tenho experiência com linguagens de programação, então > >>> embora nova no R acredito que eu tenha uma bagagem para contribuir com > a > >>> comunidade. > >>> > >>> Uma das minhas áreas de interesse além da TI é saúde, e eu vi que > várias > >>> pessoas já perguntaram nesta lista sobre uma forma de abrir os > arquivos do > >>> DATASUS no R nativamente. Eu já tive oportunidade de trabalhar com o > TabWin > >>> no passado e não era muito feliz, sempre fui a favor de plataformas > mais > >>> flexíveis e abertas. > >>> > >>> Por este motivo eu criei uma package para o R chamada "read.dbc" que > >>> permite abrir os arquivos DBC nativamente, sem depender do ferramental > que > >>> o DATASUS disponibiliza. Inclusive, por ser compilada diretamente no > R, ela > >>> pode ser utilizada em Linux sem precisar de Wine ou qualquer outra > >>> ferramenta. > >>> > >>> Esta package é completamente open source e pode ser utilizada por > >>> qualquer pessoa que se interessar. Caso a utilizem nos seus trabalhos, > >>> porém, só solicito que o meu nome seja creditado de acordo com a > licença > >>> AGPL-3. > >>> > >>> A quem interessar, segue o link para o código fonte/download: > >>> https://github.com/danicat/read.dbc > >>> > >>> Att., > >>> Daniela > >>> > >>> *_______________________________________________* > >>> 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. > >> > >> > >> _______________________________________________ > >> 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/20160606/9e94a683/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 66, assunto 5 > **************************************** > -- ======================================= Cassiano S. Rosa CRBio: 104117/04-D Universidade Federal do Triângulo Mineiro (UFTM), *campus* Iturama-MG Avenida Rio Paranaíba, 1241, Centro, 38280-000 Tel: (34)3415-2512 Lattes <http://lattes.cnpq.br/3030825329622014> Scholar Google <http://scholar.google.com/citations?user=JmELzOMAAAAJ> Academia.edu <https://uftm.academia.edu/CassianoRosa> Researchgate <https://www.researchgate.net/profile/Cassiano_Rosa> =======================================

2 1

Leitura de Arquivo
by Anselmo Alves de Sousa 07 Jun '16

07 Jun '16

Boa noite listeiros Estou com a missão de auxiliar um colega no trabalho escolar. A primeira missão é ler o dado. Como begginer no R não estou conseguindo ler o arquivo ftp://ita.ee.lbl.gov/traces/BC-pOct89.TL.Z. Sabemos que tem duas colunas, uma com o timestamp e outra com o tamanho em bytes. Usando dataset<-scan("C:\\Users\\Anselmo\\Desktop\\trab\\BC-pOct89.TL\\BC-pOct89.TL") no diretório apropriado setwd("C://Users//Anselmo//Desktop//trab_vinicius//BC-pOct89.TL") e str(dataset) num [1:2000000] 0.0177 87 0.0368 142 0.0368 ... Ou seja, obtenho apenas as duas colunas numa só. Como posso ler adequadamente esse arquivo? []'s Anselmo

4 4

R e DATASUS
by Daniela Petruzalek 06 Jun '16

06 Jun '16

Boa tarde pessoal! Eu sou nova na lista mas tive a oportunidade de ler algumas discussões pelo histórico na WEB e resolvi me juntar ao grupo. Por isso, antes de mais nada, uma rápida introdução: eu sou nova na linguagem R (e em Data Science de modo geral) mas tenho estudado bastante sobre o tema, visto que este se tornou um dos principais temas em tecnologia nos últimos dois anos. Eu sou especialista em banco de dados relacionais (principalmente Oracle) e também tenho experiência com linguagens de programação, então embora nova no R acredito que eu tenha uma bagagem para contribuir com a comunidade. Uma das minhas áreas de interesse além da TI é saúde, e eu vi que várias pessoas já perguntaram nesta lista sobre uma forma de abrir os arquivos do DATASUS no R nativamente. Eu já tive oportunidade de trabalhar com o TabWin no passado e não era muito feliz, sempre fui a favor de plataformas mais flexíveis e abertas. Por este motivo eu criei uma package para o R chamada "read.dbc" que permite abrir os arquivos DBC nativamente, sem depender do ferramental que o DATASUS disponibiliza. Inclusive, por ser compilada diretamente no R, ela pode ser utilizada em Linux sem precisar de Wine ou qualquer outra ferramenta. Esta package é completamente open source e pode ser utilizada por qualquer pessoa que se interessar. Caso a utilizem nos seus trabalhos, porém, só solicito que o meu nome seja creditado de acordo com a licença AGPL-3. A quem interessar, segue o link para o código fonte/download: https://github.com/danicat/read.dbc Att., Daniela

5 8

Estatística circular - gráfico
by Cassiano 06 Jun '16

06 Jun '16

Prezados, Tenho uns dados coletados durante o ano e fiz análise circular. Quanto à análise, está tranquilo, mas não consigo plotar no gráfico a variação em torno da média. Com o "arrows.circular" eu consegui inserir a média no gráfico, mas não encontrei nada para incluir a variação em torno da média. Desde já agradeço. cassiano > > -- ======================================= Cassiano S. Rosa CRBio: 104117/04-D Universidade Federal do Triângulo Mineiro (UFTM), *campus* Iturama-MG Avenida Rio Paranaíba, 1241, Centro, 38280-000 Tel: (34)3415-2512 Lattes <http://lattes.cnpq.br/3030825329622014> Scholar Google <http://scholar.google.com/citations?user=JmELzOMAAAAJ> Academia.edu <https://uftm.academia.edu/CassianoRosa> Researchgate <https://www.researchgate.net/profile/Cassiano_Rosa> =======================================

3 2

Trabalhar com planilha de dados muito grande
by Andre Oliveira 06 Jun '16

06 Jun '16

bom dia,estou precisando filtrar informação de arquivos de microdados do INEP, tem arquivos csv e txt com até 6gb a título basicamente de estudos. Sou leigo neste tema e gostaria de saber qual melhor forma de proceder. Me recomendaram usar um banco de dados neste caso (SQLite, MySQL, PostgreSQL ou Oracle), mas não sei por onde começar. Alguém tem material que julgue importante e que possa compartilhar comigo? 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

6 7

Espacialização randomica de pontos
by Rudiney S Pereira 04 Jun '16

04 Jun '16

Prezados colegas, Solicito auxilio no seguinte assunto: Possuo uma planilha contendo contendo as informações de N.do campo, Área e Número de amostras. Preciso de uma rotina que espacialize os pontos amostrais previamente determinados para cada área de campo de forma aleatória. Estrutura basica dos dados segue nessa forma: N.Campo; Área(m2); Quant.de Pontos Amostrais 1 43,5 8 2 44 7 3 36 5 ... ... ... 10 87 14 NOTA: tenho as coordenadas limites de cada polígono que define a área de cada campo. Desde já agradeço! -- Prof. Dr. Rudiney Soares Pereira UFSM - CCR - DER - Núcleo de Geotecnologias Prédio 44 J, Sala 214 RS 509 KM 9 - Campus Universitário 97.105-900 Santa Maria RS Brasil site: http://www.ufsm.br/geosere skype: rudiney.s.pereira fone: 55 3220-9468 55 9118-4142

3 2

Teste
by Rodrigo Coster 04 Jun '16

04 Jun '16

Alguém recebeu?

4 3

[OFF-TOPIC] R-consortium
by Felipe 30 Mai '16

30 Mai '16

Prezados colegas, o que acham deste post: https://www.r-consortium.org/news/blogs/2015/08/best-practices-using-r-secu… -- Atenciosamente Felipe E. Barletta Mendes Estatístico(UFPR) - Conre3 9766-A Mestrando em Bioestatística(UEM) +55 (41)-92077191 +55 (41)-33287216

2 1

Output de um loop em uma matriz vazia
by Yury Duarte 30 Mai '16

30 Mai '16

Boa tarde colegas programadores! Estou tendo dificuldades em acertar o formato que devo organizar os outputs de um looping que fiz. Tenho um banco de dados com datas e algumas variáveis climáticas. Dessas datas, preciso simular alguns cálculos para dias específicos. São 36 simulações por ano, durante um período de 34 anos, totalizando 1224 simulações. Para alcançar esse objetivo estou trabalhando com 3 comandos for, um dentro do outro. O primeiro seleciona as datas/posições que devem ser feitos os cálculos. O segundo executa os devidos cálculos para as datas selecionadas e o ultimo integra os valores calculados, gerando a informação de interesse. Acredito que a logica de cálculos dentro de cada for() está ok, mas a saída não está puxando os valores que o looping calcula. Segue o script no corpo do email e, em anexo, um exemplo de um ano do banco de dados que está sendo utilizado. Desde já, agradeço a colaboração dos colegas. Abs #--------------------Remover objetos salvos na memória do R--------------------# rm(list = ls()) #--------------------Seleção de diretorios e arquivos de interesse--------------------# a = read.table("x.txt", sep = "\t", header = T) #--------------------Calculo das Variaves Astronomicas--------------------# head(a) a$DATE = as.Date(levels(a$DATE))[a$DATE] ano = format(a$DATE, trim = T, '%Y'); a$ano=ano mes = format(a$DATE, trim = T, '%m') ; a$mes=mes dia = format(a$DATE, trim = T, '%d') ; a$dia=dia dece = ifelse(dia<11, 1,ifelse(dia<21,2,3)); a$decendio=dece lat = -7.53 corr = pi/180 decl = 23.45*sin(corr*((a$DiaJuliano-80)*360/365)); decl hn = 1/corr*acos(-tan(corr*lat)*tan(corr*decl)) a$N = 2*hn/15 s = (1+0.033*cos(a$DiaJuliano*360/365)) t = corr*hn*sin(corr*lat)*sin(corr*decl) u = cos(corr*lat)*cos(corr*decl)*sin(corr*hn) a$Qo = 37.6*s*(t+u) head(a) x = 0.25 b = 0.50 n.est = (a$SRAD/a$Qo-x)*a$N/b; n.est a$n<- ifelse((n.est)<0,1, n.est); head(a) a$n.N = a$n/a$N; head(a) #--------------------Transformar MJ em caloria--------------------# a$Srad.cal = round((a$SRAD*1000000)/41800,1); head(a) #round é para delimitar o numero de casas apos a virgula a$sem = paste0(dia, '-', mes) a$TMED = (((a$TMAX+a$TMIN))/2) head(a) #-------------------- Determinacao das datas de simulacao--------------------# ndc = 140 datas = 1:366 # referentes aos dia juliano dim(a) anos = table(a$ano); anos potencial = matrix(NA, nrow = 366, ncol = 34 ) rownames(potencial) = 1:366 colnames(potencial) = c(1980:2013) head(potencial); dim(potencial) head(a) for (i in (a$DiaJuliano)){ epocas = a[which( a$sem == '01-01'|a$sem == '11-01' |a$sem =='21-01'| a$sem == '01-02'|a$sem == '11-02' |a$sem =='21-02'| a$sem == '01-03'|a$sem == '11-03' |a$sem =='21-03'| a$sem == '01-04'|a$sem == '11-04' |a$sem =='21-04'| a$sem == '01-05'|a$sem == '11-05' |a$sem =='21-05'| a$sem == '01-06'|a$sem == '11-06' |a$sem =='21-06'| a$sem == '01-07'|a$sem == '11-07' |a$sem =='21-07'| a$sem == '01-08'|a$sem == '11-08' |a$sem =='21-08'| a$sem == '01-09'|a$sem == '11-09' |a$sem =='21-09'| a$sem == '01-10'|a$sem == '11-10' |a$sem =='21-10'| a$sem == '01-11'|a$sem == '11-11' |a$sem =='21-11'| a$sem == '01-12'|a$sem == '11-12' |a$sem =='21-12'& a$ano>=1980),] semear = length(epocas[,1]) for (j in length(semear)){ ini = as.numeric(rownames(epocas[j])) fim = (ini+ndc)-1 pac = a[(ini:fim),] #--------------------Inserindo as Produtividades Parciais--------------------# #----------Inserindo influencia da nebulosidade----------# pac$ctc[j] = -4.16+(0.4325*pac$TMED[j])-(0.00725*((a$TMED[j])^2)) pac$ctn[j] = -1.064+(0.173*a$TMED[j])-(0.0029*((a$TMED[j])^2)) pac$PPBc[j] = (107.2+(0.36*pac$Qo[j]*1000000)/41800)*pac$ctc[j]*pac$n.N[j] pac$PPBn[j] = (31.7+(0.219*pac$Qo[j]*1000000)/41800)*pac$ctn[j]*(1-(pac$n.N[j])) pac$PPB[j] = round((pac$PPBc[j])+(pac$PPBn[j]),1); head(pac) #--------------------Calculo da Produtividade - base seca --------------------# pac$c.resp[j] = ifelse(a$TMED[j] < 20, 0.6, 0.5) pac$c.colheita[j] = 0.5 umidade = 13 c.um = (1/(1-0.01*umidade)) pac$iaf[j] = 5 pac$c.iaf[j] = (0.0093+0.185*pac$iaf[j])-(0.0175*pac$iaf[j]^2) pac$PP.Liquida_seca = NA pac$PP.Liquida_seca[1] = round(pac$PPB[1]*pac$c.resp[1]*pac$c.colheita[1]*pac$c.iaf[1],1); head(pac) # em kg/ha for (w in 2:ndc){ pac$PP.Liquida_seca[w] = round(pac$PPB[w]*pac$c.resp[w]*pac$c.colheita[w]*pac$c.iaf[w]+ pac$PP.Liquida_seca[w-1],1) } #--------------------Soma da PP--------------------# pac$PP_final = round(pac$PP.Liquida_seca[w]*c.um,1) dim(pac) PP = pac[ndc,32] potencial[i,j] = PP } print(i) } View(pac) View(potencial) Yury Duarte Engenheiro Agrônomo - ESALQ/USP

1 0

aceitação
by Edimeire Alexandra Pinto 30 Mai '16

30 Mai '16

Bom dia. Obrigada pela aceitação

1 0

Contar casos em cada variável e em todas as variáveis
by Amiko Bh 28 Mai '16

28 Mai '16

Pessoal, Dado o banco: *ATO_1* *ATO_2* *ATO_3* ART 121 ART 155 ART 180 ART 33 ART 28 ART 10 ART 121 ART 15 ART 155 ART 155 ART 33 ART 121 ART 28 ART 33 ART 33 ART 121 ART 121 ART 121 ART 15 ART 121 ART 180 a) Preciso contar cada elemento em cada variável, para que vislumbre o seguinte resultado: *ATO_1* ART 121 = 3 ART 33 = 1 ART 155= 1 ART 28= 1 ART 15 = 1 e assim com *ATO_2 e **ATO_3.* b) preciso saber o total de cada elemento em todos as colunas, assim: ART 121 = 7 ART 33 = 4 ART 155= 3 ART 28= 2 ART 15 = 2 ART 180=2 ART 10=1 Muito obrigado.

4 5

Colocar legenda sem borda
by Andre Oliveira 27 Mai '16

27 Mai '16

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

3 6

Como colocar barras justapostas com o barplot
by Andre Oliveira 27 Mai '16

27 Mai '16

Pessoal já tentei, mas não consegui! Queria algo do tipo! barplot(dadosXX,beside=T,col = heat.colors(21),space=1) dadosXX= data.frame(Faixa, TFD,TFDMG)dadosXX Faixa TFD TFDMG 1 0-4 anos 3.35 6.52 2 5-9 anos 4.61 7.35 3 10-14 anos 4.09 8.42 4 15-19 anos 5.03 8.37 5 20-24 anos 2.73 7.52 6 25-29 anos 4.19 7.36 7 30-34 anos 6.08 7.53 8 35-39 anos 4.51 7.06 9 40-44 anos 9.22 6.80 10 45-49 anos 7.13 6.72 11 50-54 anos 7.65 5.93 12 55-59 anos 11.11 5.13 13 60-64 anos 7.76 4.21 14 65-69 anos 8.07 3.41 15 70-74 anos 6.60 2.90 16 75-79 anos 3.67 2.22 17 80-84 anos 2.41 1.43 18 85-89 anos 1.15 0.65 19 90-94 anos 0.42 0.31 20 95-99 anos 0.10 0.14 21 100 anos ou mais 0.10 0.03 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 2

Moving Average, Modelos mistos
by Fernando Antonio de souza 27 Mai '16

27 Mai '16

Boa tarde, colegas Estou trabalhando com modelos mistos, com dados de medidas repetidas no tempo. Ao definir a estrutura da matriz de covariância, testei a estrutura autoregressiva de ordem 1 AR1(), Moving average ARMA(p=0,q=?), onde testei para ordem de 1 a 6 (acima de 6 o modelo não convergiu). Acontece que com a ordem 6 obtive o melhor modelo (segundo anova()), embora nem todas as correlações foram resolvidas, conforme o gráfico (anexo) Ainda não me sinto confortal em trabalhar com estas estruturas. Faz sentido uma ordem de 6 em ARMA? Vocês teriam algum material que me explique o uso destas estruturas? 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

Erro ao fazer download de PDFs: failed, no valid url links detected
by Elias Carvalho 27 Mai '16

27 Mai '16

Bom dia Pessoal Estou baixando tentando baixar alguns pdf de artigos científicos pelo DOI usando o pacote metagear, porém a maioria apresenta a mensagem "failed, no valid url links detected", ou "cannot open: HTTP status was '404 Not Found' PDF download... skipped", mesmo ligando o proxy da universidade. De 1800 artigos que preciso consegui por volta de uns 200 download apenas. Alguém poderia me dar uma ajuda em como resolver este problema ? -- Best regards... 8^) “The mind that is open to new ideas never come back to its original size” *Albert Einstein* -- Obrigado Elias

3 6

ajuda com implementação de Programação Linear com o R - constrOptim
by Cleber N.Borges 26 Mai '16

26 Mai '16

Prezados, Bom dia! Como na lista tem pessoas com conhecimento e experiência, peço algumas orientações... :-) Imagino que o método de "Programação Linear (LP)" deve servir para resolver meu problema de achar os parâmetros de uma dada função objetivo. Sei que no R tem a opção nativa "constrOptim" (poderia ser outra, posso instalar pacotes) Entretanto eu não sei como elaborar minhas condições de contorno em formato matricial que é o modo como o comando constrOptim trabalha (até onde consegui entender) Agradeceria muito se alguém puder me ajudar, respondendo se eu realmente estou certo em pensar que a LP é a ferramenta certa para meu caso e /(ou que) me indique um material básico (introdutório). Abaixo segue um exemplo simplificado do tipos de função que pretendo obter parâmetros. Agradeço antecipadamente pela ajuda e por qualquer dica Cleber Borges ##################### # fobjetivo <- function( x = 0.1 ){ k1 <- 0.007585776 k2 <- 6.16595e-08 k3 <- 2.137962e-13 k4 <- 1e-14 k <- k1 * k2 * k3 * k4 # Determinar as 6 variáveis: a, b, c, d, e, f ################### # 6 Condições de contorno: ( Igualdades ) ################### k1 == a * d / c # 1 k2 == a * e / d # 2 k3 == a * f / e # 3 k4 == a * b # 4 4 * a == b + d + e + f # 5 x == c + d + e + f # 6 return( c( a, b, c, d, e, f ) ) } --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

3 4

Substituindo caracteres em palavras
by Wagner Tassinari 24 Mai '16

24 Mai '16

Olá pessoal, Tenho a seguinte amostra com os nomes: [1] "PIETRO LUCAS DA SILVA NAO MORA" "ROSECLEIA DE JESUS PORTUGAL NAO MORA MA" [3] "TAINA OLIVEIRA NAO MORA" "DANIELLE DE ALMEIDA MELO DA SILVA NAO E" [5] "DAVI FERREIRA DE CARVALHO NAO RESIDE NE" "EDNALDO DA SILVA NAO RESIDE MAIS NA ARE" [7] "MARCIA FERREIRA DE MORAIS NAO RESIDE NE" "MARIA JOANA LUDGERIO DE PAULO NAO MORA" [9] "ROSANGELA DE BRITO NAO R" "NATANAEL DOS SANTOS MACENA NAO MORA MAIS" Esses nomes apresentam imperfeições no preenchimento, daí eu gostaria de excluir tudo que está após a palavra NAO e criar outra variável para os nomes ficarem limpos dessa forma: [1] "PIETRO LUCAS DA SILVA" "ROSECLEIA DE JESUS PORTUGAL" [3] "TAINA OLIVEIRA" "DANIELLE DE ALMEIDA MELO DA SILVA" [5] "DAVI FERREIRA DE CARVALHO" "EDNALDO DA SILVA" [7] "MARCIA FERREIRA DE MORAIS" "MARIA JOANA LUDGERIO DE PAULO" [9] "ROSANGELA DE BRITO " "NATANAEL DOS SANTOS MACENA" Obrigado pessoal ----- 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

como fazer plot de componentes principais em cores?
by Jackson Rodrigues 24 Mai '16

24 Mai '16

Oi pessoal, Eu preciso fazer um PCA que varie as cores em função dos anos. O mais antigo em vermelho e a medida que chega as datas mais recentes vai ficando roxo (uma escala de cor). eu tenho esse exemplo hipotético. mat5 <- matrix(rnorm(2000), ,4) # minha matriz mat5 ages<-seq(1, 50000, by = 100) #idades ages library(plotrix) agecol<-color.scale(ages,extremes=c("purple","red")) # encaixando a escala de cor nas idades rownames(mat5) <- c(seq ) # fazendo as idades os nomes da linhas da matriz library("vegan") test.pca<-rda(mat5) # calculando componentes principais plot(test.pca,display="sites", cex=1, type="p",scaling=-3) # plotando componentes princiapais. O problema é: Como fazer os pontos do gráfico ficarem coloridos como estabelecido em " agecol<-color.scale(ages,extremes=c("purple","red")) " Agradeço desde já a ajuda Abraços Jackson

3 2

Subistituindo Caracteres em palavras
by Wagner Tassinari 24 Mai '16

24 Mai '16

Olá pessoal, Tenho a seguinte amostra com os nomes: [1] "PIETRO LUCAS DA SILVA NAO MORA" "ROSECLEIA DE JESUS PORTUGAL NAO MORA MA" [3] "TAINA OLIVEIRA NAO MORA" "DANIELLE DE ALMEIDA MELO DA SILVA NAO E" [5] "DAVI FERREIRA DE CARVALHO NAO RESIDE NE" "EDNALDO DA SILVA NAO RESIDE MAIS NA ARE" [7] "MARCIA FERREIRA DE MORAIS NAO RESIDE NE" "MARIA JOANA LUDGERIO DE PAULO NAO MORA" [9] "ROSANGELA DE BRITO NAO R" "NATANAEL DOS SANTOS MACENA NAO MORA MAIS" Esses nomes apresentam imperfeições no preenchimento, daí eu gostaria de excluir tudo que está após a palavra NAO e criar outra variável para os nomes ficarem limpos dessa forma: [1] "PIETRO LUCAS DA SILVA" "ROSECLEIA DE JESUS PORTUGAL" [3] "TAINA OLIVEIRA" "DANIELLE DE ALMEIDA MELO DA SILVA" [5] "DAVI FERREIRA DE CARVALHO" "EDNALDO DA SILVA" [7] "MARCIA FERREIRA DE MORAIS" "MARIA JOANA LUDGERIO DE PAULO" [9] "ROSANGELA DE BRITO " "NATANAEL DOS SANTOS MACENA" Obrigado pessoal ----- 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 2

invalid multibyte string
by Wagner Tassinari 24 Mai '16

24 Mai '16

Estou tentando verificar o número de caracteres de uma variável em doi bancos de dados. Um banco eu consigo e no outro me dá o erro: invalid multibyte string, element 593 Alguém saberia como posso lidar com isso ? Obrigado e abs --- Segue o comando reproduzido abaixo: > table(nchar(bancolink3$nom_pes)) # Observando o numero de caracteres do nome da pessoa 3 4 5 6 7 8 9 10 11 12 13 14 15 16 2 1 1 4 6 21 209 843 2156 4112 6802 9330 11849 15355 17 18 19 20 21 22 23 24 25 26 27 28 29 30 20044 32021 49598 76109 105911 129330 143014 144552 131724 117983 103914 93237 85936 78879 31 32 33 34 35 36 37 38 39 40 41 42 43 44 70489 60862 49874 37813 28071 20123 14360 10916 7915 5641 3430 2594 1868 1225 45 46 47 48 49 50 51 52 53 54 55 56 57 58 860 583 392 242 185 118 63 51 20 18 14 6 6 6 59 60 61 62 63 69 70 1 5 1 1 1 1 2 > table(nchar(fichaAlink3$NOME)) # Observando o numero de caracteres do nome da pessoa Error in nchar(fichaAlink3$NOME) : invalid multibyte string, element 593 ----- 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

Re: [R-br] out of topic
by sznelwar＠uol.com.br 23 Mai '16

23 Mai '16

2 1

ANALISe DE EXPERIMENTO DE ADUBAÇÃO COM COLHEITAS ANUAIS
by Wenceslau Teixeira 23 Mai '16

23 Mai '16

Caros um colega me pediu auxilio para analisar um experimento cujas colheitas foram repetidas por nove anos. O experimento tinha cinco tratamentos (não eram doses) de fertilizantes, que foram aplicados anualmente nas mesmas parcelas, em um delinamento em blocos. Era feito plantio e colheita sucessivas (em todos os anos) de milho e de soja (por nove anos). Os parâmetros medidos foram colheitas anuais (kg ha-1), teores foliares, analise de solos. Alguem tem alguma literatura ou exemplo (script) sobre experimentos com repetição no tempo ? Grato, Wenceslau Teixeira Embrapa Solos -- Wenceslau Teixeira wgt007(a)gmail.com skype: wgt007

2 1

Enc: Re: Eixo x
by suzabio 22 Mai '16

22 Mai '16

Oi Éder, Mais cedo, eu mandei esse email com um pedaço do banco de dados em txt, mas não foi. Eu estou reenviando por que o comando labels parece suficiente para mim, mas não funcionou. Eu coloquei os comandos que vc mandou antes do comando do axis 1 e dentro dele coloquei labels = T, mas não deu certo. O outro comando deu certo, mas o gráfico ficou estranho. No comando tá dizendo que os dados têm distribuição normal? Pq eles não têm. Achei que poderia ser por isso que o gráfico ficou estranho. Obrigada e um abraço,Suzana Enviado do meu dispositivo Samsung -------- Mensagem original -------- De : suzana de souza rodrigues <suzabio(a)yahoo.com.br> Data: 21/05/2016 09:24 (GMT-03:00) Para: r-br(a)listas.c3sl.ufpr.br, comunello.eder(a)gmail.com Assunto: Re: [R-br] Eixo x Oi Éder,Vou mandar em anexo um pedaço do banco em txt. Em Sábado, 21 de Maio de 2016 8:29, Éder Comunello <comunello.eder(a)gmail.com> escreveu: Suzana, bom dia! É o tipo de situação em que você deveria fornecer uma amostra dos seus dados ou pelo menos informar a natureza deles. A julgar pelo gráfico, seu objeto "Semana" é um vetor de inteiros. Nesse caso, uma opção seria criar um novo objeto do tipo "Date", onde cada observação poderia ser a data central da semana em questão. Nesse caso, acredito que baste utilizar esse novo objeto no lugar de "Semana" com a opção axes=T, sem precisar recorrer ao comando axis(). Mas é só um palpite. ================================================ É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 20 de maio de 2016 21:22, suzabio <suzabio(a)yahoo.com.br> escreveu: 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 Enviado do meu dispositivo Samsung _______________________________________________ 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.

2 1

Criar duas novas variáveis a partir de outra
by Alexandre Loures 21 Mai '16

21 Mai '16

Boa tarde! Tenho uma variável como abaixo e quero criar duas novas variáveis com a seguinte condição: Primeira variável a ser criada (exp): para todo F_3idoXidd_2_ * a nova variável assumirá BRA e para todo F_3idoXidd_3_* a nova variável assumirá ARG e assim por diante. Segundo variável a ser criada (imp): para todo F_3idoXidd_*_2 a nova variável assumirá BRA e para todo F_3idoXidd_*_3 a nova variável assumirá ARG e assim por diante. P.S.: O asterisco (*) indica qualquer outro número, ou seja, para todo F_3idoXidd_2_* sempre a nova variável será BRA. F_3idoXidd_2_2 F_3idoXidd_2_3 F_3idoXidd_2_4 F_3idoXidd_2_5 . . . F_3idoXidd_2_27 F_3idoXidd_2_28 F_3idoXidd_2_29 F_3idoXidd_2_30 F_3idoXidd_3_2 F_3idoXidd_3_3 F_3idoXidd_3_4 F_3idoXidd_3_5 . . . F_3idoXidd_3_27 F_3idoXidd_3_28 F_3idoXidd_3_29 F_3idoXidd_3_30 . . . Alguém poderia me ajudar? Desde já muito obrigado! -- *Alexandre Rodrigues Loures*

1 0

Re: [R-br] Eixo x
by suzabio 21 Mai '16

21 Mai '16

Bom dia Éder, Realmente Semana é um vetor de números inteiros. O gráfico é número de semanas por número de internações por problemas respiratórios. Vc acha que eu devo criar esse novo objeto com o comando as.Date? É a melhor solução? Um abraço e obrigada! Enviado do meu dispositivo Samsung -------- Mensagem original -------- De : Éder Comunello <comunello.eder(a)gmail.com> Data: 21/05/2016 08:28 (GMT-03:00) Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Eixo x Suzana, bom dia! É o tipo de situação em que você deveria fornecer uma amostra dos seus dados ou pelo menos informar a natureza deles. A julgar pelo gráfico, seu objeto "Semana" é um vetor de inteiros. Nesse caso, uma opção seria criar um novo objeto do tipo "Date", onde cada observação poderia ser a data central da semana em questão. Nesse caso, acredito que baste utilizar esse novo objeto no lugar de "Semana" com a opção axes=T, sem precisar recorrer ao comando axis(). Mas é só um palpite. ================================================ É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 20 de maio de 2016 21:22, suzabio <suzabio(a)yahoo.com.br> escreveu: 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 Enviado do meu dispositivo Samsung _______________________________________________ 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.

2 2

Eixo x
by suzabio 21 Mai '16

21 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 Enviado do meu dispositivo Samsung

3 2