- R-br - DINF UFPR

Problema numérico no gnm
by Davi 08 Mai '13

08 Mai '13

Boa noite a todos! Estou tentando simular de um modelo não-linear generalizado e depois ajustar esse modelo aos dados gerados na simulação. A variável resposta pode seguir Poisson ou Binomial, não tem problema. A média é modelada por: n*(1-exp(-(b0+b1*x))*(exp(-gx) ), em que x é uma variável quantitativa, n pode ser uma variável offset, no caso da Poisson, ou o número de ensaios Bernoulli, no caso da binomial. Aí vai o exemplo que vem dando erro, que eu suponho ter a ver com o optim, devido aos valores muito pequenos dos parâmetros b0, b1 e g e do valor da offset n ser muito grande. Agradeço a todos pela atenção! Abraços. # no R require(gnm) x <- sort( rep(0:5,3) ) b0 <- 1.5*1e-6 b1 <- 3*1e-6 g <- 3*1e-1 n <- rep(1e+8,length(x)) p <- (1 - exp(-(b0+b1*x)) ) * exp(-g*x) ; p mu <- n*p ## Funções não lineares # Parametrização original margo.betas <- function(x){ list(predictors=list(beta0=1,beta1=1,gamma=1), variables=list(substitute(x)), term=function(predictors,variables){ pred<-paste("(1-exp(",variables[1],"+",variables[2],"*",predictors[1],"))*(e xp(", variables[3],"*",predictors[1],"))",sep="") }) } class(margo.betas) <- "nonlin" # Reparametrização, com b0 = 1/(ln k0), b1=1/(ln k1) e g=1/(ln k2) margo.ks <- function(x){ list(predictors=list(k0=1,k1=1,k2=1), variables=list(substitute(x)), term=function(predictors,variables){ pred<-paste("(1-",variables[1],"*(",variables[2],"^",predictors[1],"))*(", variables[3],"^",predictors[1],")",sep="") }) } class(margo.ks) <- "nonlin" ######## ABORDAGEM POISSON set.seed(0119) y <-rpois(length(x),lambda=mu) eta.margo1 <- formula(y ~ -1 + Mult(offset(n), margo.ks(x) ) ) eta.margo2 <- formula(y ~ -1 + Mult(offset(n), margo.betas(x) ) ) gnm(eta.margo1,family=poisson,tolerance=1e-50) gnm(eta.margo2,family=poisson,tolerance=1e-50) ######## ABORDAGEM BINOMIAL set.seed(0119) y <-rbinom(length(x),size=n,prob=p) eta.margo1 <- formula(cbind(y,n-y) ~ -1+ margo.ks(x) ) eta.margo2 <- formula(cbind(y,n-y) ~ -1+ margo.betas(x) ) gnm(eta.margo1,family=binomial) gnm(eta.margo2,family=binomial) ___________________________________________________ Davi Butturi-Gomes Doutorando em Estatística e Experimentação Agronômica - ESALQ /USP - Piracicaba Mestre em Biometria - IB/UNESP - Botucatu Ecólogo - IB/UNESP - Rio Claro Skype: davibg87 (19)8206.6565

1 0

Krigagem
by Arthur Telles Calegario 08 Mai '13

08 Mai '13

Bom dia Senhores,Estou com um pequeno problemas em relação a Krigagem e talvez os senhores tenham uma resposta.Na verdade o problema esta em preparar os dados. Estou trabalhando com evapotranspiração no Rio Doce e possuo as coordenadas das estações meteorológicas em graus. Contudo, para a construção do semivariograma, necessito que a distância entre as mesmas seja linear.Assim devo projeta-las. É um procedimento simples quando todos as estações estão em um mesmo 'datum'. Contudo, as estações que estão do 'datum' 23S são insuficientes para a construção do semivariograma. Assim necessito saber como faço para linearizar as distâncias entre estações meteorológicas cujos limites extrapolam um 'datum' apenas. Entenderam? rs. Agradeço desde já,Arthur Calegario.

11 33

Krigagem
by Tito Conte 03 Mai '13

03 Mai '13

Bom dia, Alguém poderia me ajudar com krigagem? Estou usando o fields, e tenho uma matriz com 57703, com x, y, z pontos, quando uso o comando Krig ele da vários erros. O comando é K = Krig(xy,z) xy são as posições e z é a altura Tito Conte

4 3

pp-plot
by Julia Maria 02 Mai '13

02 Mai '13

Pessoal, o R tem alguma forma pré definida para se fazer um Probability-Probability plot? Por exemplo, o pacote circular tem a função: pp.plot(x), que verifica o ajuste a uma distribuição von-mises. E se eu quiser fazer o mesmo com outras ditribuições? Grata

1 0

set.seed
by Daniel Marcelino 02 Mai '13

02 Mai '13

Caros, estou realizando um experimento e acabei constatando que a semente que eu distribuir é fundamental para eu obter as amostras binárias. Eu estou interessado em distribuir 25 unidades para tratamento e o restante para controle. Há alguma maneira de saber, sem ser por tentativa e erro, que minha distribuição terá 25 unidades em cada categoria? Por exemplo, > # 50 Bernoulli trials > set.seed(51) > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50)) > table(T) T 0 1 22 28 > set.seed(12345) > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50)) > table(T) T 0 1 23 27 > set.seed(123) > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50)) > table(T) T 0 1 25 25 > set.seed(2) > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50)) > table(T) T 0 1 25 25 -- \begin{signature} Daniel Marcelino ☁ dm.silva(a)umontreal.ca ☎ (514) 343 6111 #3799 Skype: d.marcelino ✎ 3200 Jean Brillant, Office C5071 Montreal, QC; H3T 1N8 Canada \end{signature} "Small steps toward a much better world"

4 5

excluir linhas que tenham códigos repetidos
by Fátima Lima Paula 01 Mai '13

01 Mai '13

Pessoal, existe algum comando que exclua as linhas que apresentam códigos repetidos? Não seria deixar apenas um não. É excluir todas. Exemplo nome=c("Ricardo","Maria","Tereza","Pedro","José", "Germano","Paulo","Ana") sobrenome=c("Cunha","Andrade","Silva","Soares","Silva","Lima","Pitanga","Amora") aih=c(6867,6798,8903,7865,8903,6745,6867,8903) ano=c(1944,1945,1923,1954,1947,1945,1946,1944) dataint=c(20121004,20120907,20120703,20120301,20120805,20121209,20120907,20120703) sexo=c("M","F","F","M","M","M","M","F") banco1=cbind(nome,sobrenome,aih,sexo,ano,dataint) banco1=data.frame(banco1) banco1 nome sobrenome aih sexo ano dataint 1 Ricardo Cunha 6867 M 1944 20121004 2 Maria Andrade 6798 F 1945 20120907 3 Tereza Silva 8903 F 1923 20120703 4 Pedro Soares 7865 M 1954 20120301 5 José Silva 8903 M 1947 20120805 6 Germano Lima 6745 M 1945 20121209 7 Paulo Pitanga 6867 M 1946 20120907 8 Ana Amora 8903 F 1944 20120703 Quero excluir as linhas que apresentam as aihs repetidas (6867 e 8903), de forma que o banco fique assim: nome sobrenome aih sexo ano dataint 1 Maria Andrade 6798 F 1945 20120907 2 Pedro Soares 7865 M 1954 20120301 3 Germano Lima 6745 M 1945 20121209 Obrigada -- "Minha felicidade depende da qualidade dos meus pensamentos"

5 5

Função maxDrawdown
by Marcos Vagner Macedo 01 Mai '13

01 Mai '13

Boa noite pessoal, Não estou conseguindo utilizar a função maxDrawdown do pacote PerformanceAnalytics. Alguém pode me ajudar? Estou utilizando dias úteis no formato zoo. Quando utilizo a função tenho o seguinte resultado: > t(round(maxDrawdown(fi_97_12[,2:4]),4)) Worst Drawdown Ibovespa NaN ALFA.DINAMICO.FIA NaN ALFA.FIA NaN Warning messages: 1: In merge.zoo(fx, .xts(, .index(x))) : Index vectors are of different classes: integer POSIXct 2: In merge.zoo(fx, .xts(, .index(x))) : Index vectors are of different classes: integer POSIXct 3: In merge.zoo(fx, .xts(, .index(x))) : Index vectors are of different classes: integer POSIXct > Muito Obrigado, Marcos Macedo

1 0

Juntar colunas
by Edson Lira 01 Mai '13

01 Mai '13

Caros amigos, tenho um banco de dados aonde as quatro primeiras colunas, por problema de importação é o nome do indivíduo do banco, qual a rotina que posso juntar as colunas com os respectivos espaços. Exemplo: C1 C2 C3 C4 ADRIANA DE SOUZA MALHEIROS [ ]'S. Edson Lira Estatístico Manaus-Amazonas

4 5

package Babel com knitr [Resolvido]
by Fernando Colugnati 30 Abr '13

30 Abr '13

Oi Benilton, resolveu com o latin1 mesmo. Pouco antes de sua resposta tinha conseguido um exemplo! Abs e obrigado -- Fernando A.B. Colugnati

1 0

package Babel com knitr
by Fernando Colugnati 30 Abr '13

30 Abr '13

Prezados, estou com um problema na questão de internacionalização com o Latex, via R, usando o knitr (não o sweave). Apesar de incluir \usepackage[portuguese]{babel} não esto tendo sucesso na compilação para PDF com o pdfLaTex, os caracteres saem errados. Em tempo, estou rodando pelo Rstudio, em um MacPro (pensei em ser algum problema de codificação...). A instalação que tenho do Tex é a padrão para Mac. Alguém já passou por isso? Segue um pedaço de código, mas com todo o preambulo que estou usando: \documentclass{article} \usepackage[portuguese]{babel} \usepackage[sc]{mathpazo} \usepackage[T1]{fontenc} \usepackage{geometry} \geometry{verbose,tmargin=1.5cm,bmargin=1.5cm,lmargin=2.5cm,rmargin=2.5cm} \setcounter{secnumdepth}{2} \setcounter{tocdepth}{2} \usepackage{url} \usepackage[unicode=true,pdfusetitle, bookmarks=true,bookmarksnumbered=true,bookmarksopen=true,bookmarksopenlevel=2, breaklinks=false,pdfborder={0 0 1},backref=false,colorlinks=false] {hyperref} \hypersetup{ pdfstartview={XYZ null null 1}} \usepackage{breakurl} \begin{document} \section*{Matemática $5^o$ ano} Isso \'{e} um teste para acentuação! acentua\'{c}\˜{a}0, assim est\'{a} correto! \end{document} Obrigado -- Fernando A.B. Colugnati

2 1

Has R-help gotten meaner over time?
by Benilton Carvalho 30 Abr '13

30 Abr '13

vai que alguem quer tentar algo parecido para a r-br... http://badhessian.org/2013/04/has-r-help-gotten-meaner-over-time-and-what-d… b

1 0

Series temporais
by Edmar Caldas 30 Abr '13

30 Abr '13

Boa Tarde! Mais uma vez vou lançar uma pergunta na lista e espero sinceramente que alguém responde pelo menos não sei fazer. Sou especialista em SPSS e MODELER mas tenho que fazer algumas coisa no R. Por exemplo tenho um arquivo que preciso ajustar um arima. Preciso saber qual o comando ou a linha que gera a projeção e o MAPE. Arquivo começa em 01-01-2008 e termina em 31-12-2012 Preciso da projeção do dia 01-01-2013 até 31-01-2013 por dia tirando sabados e domingos. E o MAPE. Alguém pode me ajudar, se precisar do arquivo eu mando. Alguém tem algum código? Grato. Edmar

2 1

Re: [R-br] Manipulação de dados diários
by Marcos Vagner Macedo 29 Abr '13

29 Abr '13

Daniel, Segue um exemplo dos dados para tornar o problema reproduzível. Data CDI Ibovespa ALFA DINAMICO FIA ALFA FIA 01/02/97 100.244198 6956 77.5631 0.33373 01/03/97 100.32139 6998 78.05032 0.34725 01/06/97 100.398636 7160 79.3197 0.34982 01/07/97 100.475946 7297 82.01036 0.34084 01/08/97 100.553312 7427 82.28164 0.34426 01/09/97 100.630733 7454 82.79493 0.34285 01/10/97 100.70889 7478 82.04698 0.34263 01/13/97 100.78711 7412 82.32868 0.34275 01/14/97 100.865386 7621 84.52306 0.34776 01/15/97 100.943725 7668 86.36466 0.34922 01/16/97 101.022127 7656 86.74005 0.3513 01/17/97 101.100925 7816 86.74305 0.35152 01/20/97 101.179779 7835 86.84854 0.35297 01/21/97 101.258363 7905 87.77227 0.35453 01/22/97 101.337011 7818 87.63108 0.35711 01/23/97 101.416054 7891 87.96762 0.3565 01/24/97 101.495494 7760 87.7515 0.3547 01/27/97 101.574664 7759 87.71396 0.35453 01/28/97 101.655242 7850 87.6412 0.35487 01/29/97 101.734531 7737 86.20481 0.35258 01/30/97 101.816938 7789 86.37315 0.35225 01/31/97 101.910268 7965 88.24893 0.35641 02/03/97 102.003692 8096 89.45418 0.35969 02/04/97 102.097876 8102 89.39141 0.3595 02/05/97 102.191799 8171 90.26877 0.36263 02/06/97 102.286157 8148 90.26682 0.3611 02/07/97 102.380603 8317 92.41892 0.36529 02/12/97 102.475135 8671 96.43418 0.3736

2 1

Alocação proporcional no R
by Matheus Monteiro 29 Abr '13

29 Abr '13

Ainda trabalhando com os 50 municipios mais populosos do Brasil, divido a população em dois estratos, onde no primeiro estrato esteja os 10 mais populosos e no segundo os 40 restantes. Preciso obter uma amostra com reposição de 10 municípios usando o alocação proporcional e calcular a estimativa da média e variancia. Alguém pode ajudar? pop = c(11376685,6390290,2710968,2648532,2500532,2395785,1861838,1776761,1555039,1416014,1410430,1333767,1244518,1098630,1039610,1016128,953393,867067,830231,817590,805397,801746,774886,742428,680496,668877,654786,643603,619746,619536,613815,600692,587701,568099,561329,526338,525225,515707,491807,483821,474596,474596,472300,460062,446911,442701,433158,425169,424948,422569)

2 6

R 3.0 Ubuntu 13.04
by Frederico Alvarez 29 Abr '13

29 Abr '13

Boa tarde a todos. Estou tentando instalar o R 3.0.0 no ubuntu 13.10 via apt. Mas acontece que no repositório da UFPR não há o pacote r-base-core para 64 bit. Porém, há dois pacotes para 32 bits. Alguém sabe como posso resolver esse problema? Saudações. -- Frederico Steinmetz Alvarez Engenheiro de Produção Elétrica Mestrando em Energia Eólica Laboratório Didático Pedagógico de Fluidos - LDPFLU Departamento de Engenharia Mecânica - DEMEC Universidade Federal de Pernambuco - UFPE (81) 8885-1105 (Oi) (48) 8448-5994 (Oi) (48) 9143-5141 (Vivo)

9 11

Alocação Proporcional
by Matheus Monteiro 29 Abr '13

29 Abr '13

Ainda trabalhando com os 50 municipios mais populosos do Brasil, divido a população em dois estratos, onde no primeiro estrato esteja os 10 mais populosos e no segundo os 40 restantes. Preciso obter uma amostra com reposição de 10 municípios usando o alocação proporcional e calcular a estimativa da média e variancia. Alguém pode ajudar? Estou com dúvida aqui, pois preciso encontrar uma maneira de alocar proporcionalmente a população inteira e não apenas o estrato 2 pop = c(11376685,6390290,2710968,2648532,2500532,2395785, 1861838,1776761,1555039,1416014,1410430,1333767,1244518,1098630,1039610, 1016128,953393,867067,830231,817590,805397,801746,774886, 742428,680496,668877,654786,643603,619746,619536,613815, 600692,587701,568099,561329,526338,525225,515707,491807, 483821,474596,474596,472300,460062,446911,442701,433158, 425169,424948,422569) # order pop <- sort(pop, T) # Alocação sub1 <- pop[(1:10)] sub2 <- pop[(11:50)] # Alocação proporcional DÚVIDA prop <- sub2/sum(sub2) set.seed(51) sub3 <- sample(sub2, 10, replace=TRUE, prob=prop)

2 1

Reduzir objeto auxiliares...
by andrebvs＠bol.com.br 29 Abr '13

29 Abr '13

2 1

(sem assunto)
by Antonio Silva 29 Abr '13

29 Abr '13

Olá Eu estou tentado utilizar o teste de Qui-quadrado para verificar a dependência entre a distribuição de duas espécies, como indicado no livro do Jean Valetin (Ecologia Numérica 2a ed, pg 49). Eu tenho um dataframe com dados de presença/ausência de duas espécies em um determinado número de unidades amostrais A equação indicada no livro é: X2= p*(|ad-bc|-p/2)^2 / ((a+b)*(c+d)*(a+c)*(b+d)) Onde a é o número de dupla presença (1-1), b o número de 1-0, c o número de 0-1 e d o número de 0,0 (dupla ausência) p é a+b+c+d Há alguma função no R para este cálculo? Eu não consegui utilizar a função chisq.test para este cálculo. Obrigado! Antonio Olinto

4 4

Estratificação no R
by Matheus Monteiro 29 Abr '13

29 Abr '13

Estou trabalhando com os 50 municipios mais populosos do Brasil, foi pedido para dividir a população em dois estratos, onde no primeiro estrato esteja os 10 mais populosos e no segundo o restante. Preciso obter uma amostra de 5 municípios de cada estrato e calcular a estimativa da media amostral e variancia. Sei que preciso instalar o pacote gtools mas não estou conseguindo. Não estou conseguindo, alguém pode ajudar? Por enquanto está assim: pop = c(11376685,6390290,2710968,2648532,2500532,2395785,1861838,1776761,1555039,1416014,1410430,1333767,1244518,1098630,1039610,1016128,953393,867067,830231,817590,805397,801746,774886,742428,680496,668877,654786,643603,619746,619536,613815,600692,587701,568099,561329,526338,525225,515707,491807,483821,474596,474596,472300,460062,446911,442701,433158,425169,424948,422569) Npop=length(pop) Npop EstI=c(11376685,6390290,2710968,2648532,2500532,2395785,1861838,1776761,1555039,1416014) EstII=c(1410430,1333767,1244518,1098630,1039610,1016128,953393,867067,830231,817590,805397,801746,774886,742428,680496,668877,654786,643603,619746,619536,613815,600692,587701,568099,561329,526338,525225,515707,491807,483821,474596,474596,472300,460062,446911,442701,433158,425169,424948,422569) nEstI=length(EstI) nEstII=length(EstII)

5 8

Sequência no R...
by andrebvs＠bol.com.br 28 Abr '13

28 Abr '13

2 2

Manipulação de dados diários
by Marcos Vagner Macedo 28 Abr '13

28 Abr '13

Boa tarde, Estou com dificuldade na utilização de dados no formato zoo. Estou usando o comando abaixo. z <- read.zoo("dados.csv", header = TRUE, sep = ",", format = "%m-%d-%Y") Está aparecendo o seguinte erro. Erro em read.zoo("dados.csv", header = TRUE, sep = ",", format = "%m-%d-%Y") : index has 4016 bad entries at data rows: 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 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 ... > Vocês poderiam me ajudar? Além disso, gostaria de saber como pegar sub amostras desse banco de dados. A primeira amostra deve ser um período reduzido do banco de dados (OBS:. dados diários - dias uteis) A segunda amostra deve ser algumas séries no período reduzido. Como o banco é grande o ideal é que os filtros seja feitos de uma forma que seja possível automatizar o processo. Exemplo: 5 séries . Periodo 10/01/2000 a 28/12/2012 Analise: 3 séries Periodo 08/01/2012 a 06/08/2012

2 1

Join my network on LinkedIn
by Sérgio Henrique Almeida da Silva Junior via LinkedIn 27 Abr '13

27 Abr '13

LinkedIn ------------ Sérgio Henrique Almeida da Silva Junior requested to add you as a connection on LinkedIn: ------------------------------------------ I'd like to add you to my professional network on LinkedIn. Accept invitation from Sérgio Henrique Almeida da Silva Junior http://www.linkedin.com/e/-33g73l-hg1l8cl7-5x/vhj1xG42iBdUFK5d_HMEGL4iFMC51… View profile of Sérgio Henrique Almeida da Silva Junior http://www.linkedin.com/e/-33g73l-hg1l8cl7-5x/rso/250437145/r-Yw/name/82777… ------------------------------------------ You are receiving Invitation emails. This email was intended for Fábio Corrêa. Learn why this is included: http://www.linkedin.com/e/-33g73l-hg1l8cl7-5x/plh/http%3A%2F%2Fhelp%2Elinke… (c) 2012, LinkedIn Corporation. 2029 Stierlin Ct, Mountain View, CA 94043, USA.

1 0

Ifelse para duas condições
by Fátima Lima Paula 27 Abr '13

27 Abr '13

Preciso fazer um ifelse para duas condições. Exemplo nome=c("Ricardo","Maria","Tereza","3edro","José", "Germano") sobrenome=c("Cunha","Andrade","Silva","Soares","Silva","Lima") idade=c(68,67,89,78,65,67) idade2=c(68,67,90,77,65,70) ano=c(1944,1945,1923,1954,1947,1945) ano2=c(1953,1945,1923,1951,1947,1956) dataint=c(20121004,20120907,20120703,20120301,20120805,20121209) banco1=cbind(nome,sobrenome,idade,idade2,ano,ano2,dataint) banco1=data.frame(banco1) banco1 banco1$ano=as.character(banco1$ano) banco1$ano2=as.character(banco1$ano2) banco1$idade=as.character(banco1$idade) banco1$idade2=as.character(banco1$idade2) Quero criar uma coluna que me diga sim se as idades e os anos estiverem iguais e não caso contrário. Fiz assim, banco$teste=if(banco1$idade==banco1$idade2 & banco1$ano==banco1$ano2)sim else nao Naturalmente, deu erro. Alguém pode me ajudar, por favor? -- "Minha felicidade depende da qualidade dos meus pensamentos"

2 2

Enéias Façanha invited you to check out Dropbox
by Enéias Façanha via Dropbox 27 Abr '13

27 Abr '13

Hi there, Enéias Façanha wants you to try Dropbox! Dropbox lets you bring all your photos, docs and videos with you anywhere and share them easily. Get started here. https://www.dropbox.com/l/QAUjvCIBD3CpWoInvhv6nT1 Thanks! - The Dropbox Team ____________________________________________________ To stop receiving invites from Dropbox, please go to https://www.dropbox.com/l/23TKPtwTdCSwSn42jGZlSY1 Dropbox, Inc., PO Box 77767, San Francisco, CA 94107

1 0

Ajuda em função
by Sérgio Henrique almeida da silva ju 25 Abr '13

25 Abr '13

Prezados Estou fazendo uma função no R que dado os três lados de um triângulo ele me diz se é um triângulo, se o mesmo é isósceles, escaleno ou equilátero. tri <- function(x,y,z){ if ((x +y > z) & (x + z > y) & (y + z > x)){ if ((x==y) | (y==z)){return ("equilátero")} else if ((x = y) | (x=z) | (y=z)){ return ("isosceles")} else{ return ("escaleno")} else{ return ("não é triângulo")}} Porém não estou tendo sucesso, onde posso estar errando? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 / 94429486

4 6

Eliminar linhas cuja variável comece com uma letra
by Fátima Lima Paula 25 Abr '13

25 Abr '13

Para selecionar as linhas cujo diagnóstico principal comece com a letra F, o Gledson me passou o seguinte comando: dados[substr(dados$diag,1,1)=="F",] Agora, preciso eliminar do banco essas linhas. Pelo que vi, não dá para usar o !=. Alguém poderia me ajudar? Obrigada

2 3

Re: [R-br] Digest R-br, volume 26, assunto 37
by Diogo Braga Mendes 25 Abr '13

25 Abr '13

Use a função which que retorna os indices das linhas que esta condição foi satisfeita. exemplo: dados <- dados[ - *which*(*substr*(dados$diag,1,1) == "F"), ] Em 25 de abril de 2013 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: Eliminar linhas cuja variável comece com uma letra > (Fátima Lima Paula) > > > ---------- Mensagem encaminhada ---------- > From: "Fátima Lima Paula" <fatima.lima.paula(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Cc: > Date: Tue, 23 Apr 2013 15:37:05 -0300 > Subject: Re: [R-br] Eliminar linhas cuja variável comece com uma letra > Consegui por meio do != > > > Em 23 de abril de 2013 15:10, Fátima Lima Paula < > fatima.lima.paula(a)gmail.com> escreveu: > >> Para selecionar as linhas cujo diagnóstico principal comece com a letra >> F, o Gledson me passou o seguinte comando: >> dados[substr(dados$diag,1,1)=="F",] >> Agora, preciso eliminar do banco essas linhas. Pelo que vi, não dá para >> usar o !=. Alguém poderia me ajudar? >> Obrigada >> > > > > -- > "Minha felicidade depende da qualidade dos meus pensamentos" > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Diogo Braga Mendes ------------------------------------------------ Estatístico do IPAMB-Belém/PA e Analista em projetos. Mestre em Estatística-PPGME/UFPA. Bacharel em Estatística, CONRE 9558. Programador C++, JAVA, PYTHON, PHP e JS. Fone: (91) 8239-7799 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ “O âmago da ciência, que sempre se apresenta como a mais importante no que refere a resultados práticos, é a pesquisa altamente teórica e abstrata, nascida da infatigável curiosidade, flexibilidade e força da razão humana” Andrei Sakharov In Análise Harmônica de Processos Estocásticos

1 0

Matriz X (Y = XB + e)
by Tiago Souza Marçal 25 Abr '13

25 Abr '13

Prezados, gostaria de saber se existe uma função que simplica a obtensão desta matriz? m<-matrix(c(rep(1,15),rep(1,5),rep(0,10),rep(0,5),rep(1,5),rep(0,5),rep(0,10),rep(1,5)),nr=15) Att. Tiago.

4 4

Plot limites inferiores, superior e outliers
by T Branquinho 25 Abr '13

25 Abr '13

Caros R-users, Estou avaliando alguns modelos feitos em programas diferentes e gostaria de padronizar os gráficos conforme o R, me esbarrei com algumas dificuldades, como acrescentar o gráfico abaixo a informação dos limites inferiores e superiores (igual fiz no abline) e como identificar no gráfico apenas os outliers?. xnames <- c("1","2","3","4","5","6","7","8","9","10")x <- rnorm(10, 0, 2) y <- rnorm(10, 1, 1) plot(y,x,xlab="Residuals",ylab="Fitted values",main="Residuals vs Fitted"); abline(h = 0, lty = 3, col = "gray")

2 2

Listar só os códigos que aparecem mais de uma vez
by Fátima Lima Paula 25 Abr '13

25 Abr '13

Pessoal, estou tentando ver os códigos que aparecem mais de uma vez no banco. O problema é que na hora que faço o table, como têm muitos que aparecem uma vez, não dá espaço para ver o restante. O André me passou o comando > cbind(table(apply(cod, 1, paste, collapse=""))) Não deu certo. Alguém pode me dar uma dica? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

2 1

Utilização de qui-quadrado na análise de distribuição de espécies
by Antonio Silva 25 Abr '13

25 Abr '13

Olá Eu estou tentado utilizar o teste de Qui-quadrado para verificar a dependência entre a distribuição de duas espécies, como indicado no livro do Jean Valetin (Ecologia Numérica 2a ed, pg 49). Eu tenho um dataframe com dados de presença/ausência de duas espécies em um determinado número de unidades amostrais A equação indicada no livro é: X2= p*(|ad-bc|-p/2)^2 / ((a+b)*(c+d)*(a+c)*(b+d)) Onde a é o número de dupla presença (1-1), b o número de 1-0, c o número de 0-1 e d o número de 0,0 (dupla ausência) p é a+b+c+d Há alguma função no R para este cálculo? Eu não consegui utilizar a função chisq.test para este cálculo. Obrigado! Antonio Olinto PS. há pouco enviei a mesma mensagem para lista, porém sem o assunto. Desconsiderem a primeira e desculpem-me por qualquer incoveniente.

2 2

Fwd: Sua mensagem para a lista R-br aguarda aprovação
by Fátima Lima Paula 25 Abr '13

25 Abr '13

Acho que não vou conseguir enviar as colunas ---------- Mensagem encaminhada ---------- De: <r-br-bounces(a)listas.c3sl.ufpr.br> Data: 17 de abril de 2013 11:22 Assunto: Sua mensagem para a lista R-br aguarda aprovação Para: fatima.lima.paula(a)gmail.com Seu email para 'R-br' com o assunto Espaços em branco Esta em espera até que o moderador da lista revise-a para aprovação. A razão de estar em espera é: O corpo da mensagem é muito grande: 9015783 bytes com um limite de 150 KB Ou a mensagem será postada a lista, ou receberá uma notificação da decisão do moderador. Se desejar cancelar esta postagem, visite o seguinte endereço: https://listas.inf.ufpr.br/cgi-bin/mailman/confirm/r-br/0e28e1e8daa488c7845…

7 15

Matrix X (Y = XB + e)
by Tiago Souza Marçal 25 Abr '13

25 Abr '13

Prezados, gostaria de saber se existe uma função que simplica a obtenção desta matriz? m<-matrix(c(rep(1,15),rep(1,5),rep(0,10),rep(0,5),rep(1,5),rep(0,5),rep(0,10),rep(1,5)),nr=15) Att. Tiago.

2 1

Subset começando pela letra
by Fátima Lima Paula 23 Abr '13

23 Abr '13

Pessoal, tenho uma coluna em um data frame com códigos de diagnósticos. Quero fazer um subset dos que comecem pela letra F. Alguém pode me ajudar?

5 6

(sem assunto)
by Helena Turon 23 Abr '13

23 Abr '13

Valeu Thiago, Consegui saber o período que a estação ficou desativada, é ótimo! Como ele reclama do vetor de datas ser um fator, eu faço o comando que você ensinou no começo, quando leio os dados do arquivo. Isto diminuiu muito o tempo de programação, nem precisei dos vetores com as datas! Obrigada, Helena

1 0

Novamente... NA
by Emerson Cotta Bodevan 23 Abr '13

23 Abr '13

Prezados, tenho um data frame em que uma das variáveis (quantitativa) possui NA. Como faço para contar quantos NAs eu tenho? Como faço para a função table contá-los? Outra pergunta... qual a diferença de NA e <NA>? Percebi que NA aparece em variáveis quantitativas e <NA> aparece em variáveis qualitativas. (Carreguei meu Banco de Dados com o comando read.table). Agradeço qualquer ajuda. Att., Emerson

2 2

listar os códigos com quantidade > 1
by Fátima Lima Paula 22 Abr '13

22 Abr '13

Pessoal, tenho um banco com uma coluna com códigos. Pelo summary, verifiquei que existem códigos com mais de uma quantidade. Se faço o table, estoura a tela. Não dá para ver quais os que têm mais de um iguais. Alguém pode me dizer como faço? Obrigada Fátima aih 3308100991261: 5 3308100275557: 4 3308101004440: 4 3308101326992: 4 3308101328829: 4 3308102043906: 4 (Other) :236099

2 4

substituir datas em linhas específicas
by Fátima Lima Paula 22 Abr '13

22 Abr '13

Pessoal, quando fui substituir uns nomes no meu banco, usei levels(banco$nome) <- gsub("FATIMALIMA", "FATIMA LIMA", levels(banco$nome)) Agora o problema é outro, tenho uma linha com nome, data de nascimento, data de internação, motivo da alta etc. Em uma específica linha, quero trocar o motivo da alta e a data da internação. Não posso usar o comando acima pois existe um monte de outros motivos de saida 21, por exemplo, e eu só quero substituir para a FATIMA LIMA. Alguém, por favor, pode me ajudar? Obrigada

4 10

NA
by Emerson Cotta Bodevan 21 Abr '13

21 Abr '13

Tenho um data frame em que uma das colunas possui "sim", "não" e NA. Como faço para criar um novo data frame que só contenha as linhas em que a variável acima seja igual a "sim"? Não quero "não" e NA. Agradeço qualquer ajuda. Emerson

2 2

Re: [R-br] Contar Dados faltantes
by Helena Turon 21 Abr '13

21 Abr '13

Oi Bernardo, Os 2 vetores deveriam ser comparados. Um tem todas as datas (é teórico) e o observado tem dias faltando. Aplicando o applay, separadamente, nenhum retorna dados faltantes. A maneira que eu fiz este loop seria mais fácil diminuir o comprimento das colunas neles. Estou tentando ver se consigo resgatar o dia e hora em que houve a diferença. Senão farei esta diminuição mesmo. Abraço e obrigada. Helena

2 1

ENC: Series temporiai Arima
by edmar Caldas 21 Abr '13

21 Abr '13

Ainda não obtive resposta. De: edmar Caldas [mailto:edmar_caldas@yahoo.com.br] Enviada em: sábado, 20 de abril de 2013 20:01 Para: 'r-br(a)listas.c3sl.ufpr.br' Assunto: Series temporiai Arima dados1<-read.csv("C:\\Users\\Edmar\\Desktop\\dados.csv",head=T,sep=";",dec=" ,") > attach(dados1) > head(dados1) Data V2 1 02/jan/08 5.007.331,28 2 03/jan/08 4.698.059,06 3 04/jan/08 4.337.845,89 4 07/jan/08 4.297.929,19 5 08/jan/08 3.922.327,63 6 09/jan/08 4.346.548,66 > plot.ts(V2) At´q aqui OK. Porque não consigom o rban preciso de algum pacote adicional.? > rban<-diff(log(V2)) Erro em Math.factor(c(483L, 401L, 340L, 333L, 269L, 342L, 301L, 209L, : log não faz sentido para fatores

3 5

Contar Dados faltantes
by Helena Turon 21 Abr '13

21 Abr '13

Prezados, Para contar os horários ou dias sem observações na estação meteorológica, dados faltantes. Fiz um “for” no R: c=comprimento do vetor teórico (sequencia de datas e horários que deveriam existir observações). h= comprimento dos dados das estações (pode haver horários ou dias seguidos sem observações). igual=0 for (m in 1:c){ for (jj in 1:h){ if ((bancof[m,1] == obsf[jj,1]) & (bancof[m,2] == obsf[jj,2]) ) {igual = igual+1} } } faltantes = c-igual coluna 1 é data dd/mm/aaaa e coluna 2 horário (0, 12, 18) Mas ele está demorado: uns 20 minutos por estação (pois elas têm umas 60.000 observações). Eu tenho umas 200 estações. Serão quatro dias, só para os dados do Brasil, têm ainda outros países da América do Sul. Será que existe alguma função no R que faça isto mais depressa? Obrigada, -- Helena Turon Balbino Meteorologista Instituto Nacional de Meteorologia - INMET 7º DISTRITO DE METEOROLOGIA (SP e MS) Tel: (11) 5051-5700 ********************************************************

5 4

Regressao Logística
by ynha.thays＠gmail.com 20 Abr '13

20 Abr '13

Estou usando a regressão logística no R, gostaria de saber como uso o comando glm.fit, não estou coneguindo entender a explicacao no Manual do R. abraços Thays

1 0

Series temporiai Arima
by edmar Caldas 20 Abr '13

20 Abr '13

dados1<-read.csv("C:\\Users\\Edmar\\Desktop\\dados.csv",head=T,sep=";",dec=" ,") > attach(dados1) > head(dados1) Data V2 1 02/jan/08 5.007.331,28 2 03/jan/08 4.698.059,06 3 04/jan/08 4.337.845,89 4 07/jan/08 4.297.929,19 5 08/jan/08 3.922.327,63 6 09/jan/08 4.346.548,66 > plot.ts(V2) At´q aqui OK. Porque não consigom o rban preciso de algum pacote adicional.? > rban<-diff(log(V2)) Erro em Math.factor(c(483L, 401L, 340L, 333L, 269L, 342L, 301L, 209L, : log não faz sentido para fatores

1 0

ENC: Series Temporias
by edmar Caldas 20 Abr '13

20 Abr '13

Pessoal preciso de uma ajuda. fiz uma serie temporal ele mesmo ajustou um ARIMA. Agora preciso achar o MAPE e as estatísticas do modelo como faço. Tem também o ljung.box

1 0

Re: [R-br] mudar de decimal para inteiro
by Helena Turon 20 Abr '13

20 Abr '13

Fátima, Estou trabalhando com conjuntos de dados grandes em txt. Troquei o ponto e vírgula por vírgula. No bloco de notas mesmo (é bem rápido) e abri: liarq <- read.table(paste(x1), header=F, sep=",", na.strings="NA", dec=".", strip.white=TRUE) Abraço, Helena

1 0

Re: [R-br] Contar Dados faltantes
by Helena Turon 20 Abr '13

20 Abr '13

Prezados Edson, Fernando e Manoel, boa noite. data.table Instalei o pacote (e a versão 3.0 do R), estou aprendendo, parece bem interessante! Teoricamente c e h deveriam ter o mesmo tamanho (60.000 linhas). São 2 vetores com datas e horários: um teórico (com todas) e o outro o real, que deve ter algumas faltas. Segue amostra dos vetores. O (%in%) retorna alguns True que seriam False. Se eu comparar o tamanho dos vetores, tenho uma ideia das faltas. Agora vi que era o que eu estava fazendo com aquele loop. aa <- length(bancof[,1])-length(obsf[,1] Também, acho que meu computador é fraco: Intel(R) Core(TM)2 DUO CPU E7500(a)2.93GHz Abçs e muito obrigada. Helena

1 0

Re: [R-br] mudar de decimal para inteiro
by Fátima Lima Paula 20 Abr '13

20 Abr '13

Acabei de descobrir que quando eu leio o banco esse número "perde o final". É um código. Ele tem os primeiros dígitos depois vem vários zeros e depois mais alguns dígitos. Estou em desespero. Perdi horas e horas de trabalho. Vi que é um problema de importação do banco. Nunca havia trabalhado com bancos grandes. Quando eu peguei o banco, ele veio em txt. Deve ter sido na hora que fui abrir no excel. Alguém pode me dar uma luz em relação a isso? Tenho um banco em txt, separado por ponto e vírgula, sem "header". Como é a melhor forma de trabalhar esse banco no R? Obrigada Em 19 de abril de 2013 18:21, Leandro Marino < leandromarino(a)leandromarino.com.br> escreveu: > Em um número deste tamanho consideraria transformar para character. > Fatalmente é um código de município. > > Já tive problemas em números inteiros largos na hora de fazer o merge uma > vez que eles são interpretados como pontos flutuantes.... > > Tive dois números que apesar de serem diferentes na hora do merge eram > juntados.... > Em 19/04/2013 17:24, "Fátima Lima Paula" <fatima.lima.paula(a)gmail.com> > escreveu: > >> Prezados meus números estão aparecendo dessa forma: >> 3.31000e+12 >> Como faço para que fique >> 3310008457687 >> Alguém pode me ajudar, por favor? >> Obrigada >> >> _______________________________________________ >> 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. >

3 7

Re: [R-br] mudar de decimal para inteiro
by David Valentim Dias 19 Abr '13

19 Abr '13

Se realmente são dados inteiros deve especificar usando interger() ou as.integer(), existe um tratamento especial para já que sua representação muda dependendo de hardware/linguagem! Transformar em character resultada em um uso desnecessário de bytes para armazenar o número. Para quem acha que 0.2 + 0.1 deveria dar 0.3 *> sprintf("%0.20f", 0.2 + 0.1) [1] "0.30000000000000004441"* Recomendo ler: http://floating-point-gui.de/basic/ Em 19 de abril de 2013 17:21, Leandro Marino < leandromarino(a)leandromarino.com.br> escreveu: > Em um número deste tamanho consideraria transformar para character. > Fatalmente é um código de município. > > Já tive problemas em números inteiros largos na hora de fazer o merge uma > vez que eles são interpretados como pontos flutuantes.... > > Tive dois números que apesar de serem diferentes na hora do merge eram > juntados.... > Em 19/04/2013 17:24, "Fátima Lima Paula" <fatima.lima.paula(a)gmail.com> > escreveu: > >> Prezados meus números estão aparecendo dessa forma: >> 3.31000e+12 >> Como faço para que fique >> 3310008457687 >> Alguém pode me ajudar, por favor? >> Obrigada >> >> _______________________________________________ >> 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. > -- Biostatistician Consultant and Cogito Programmer's http://www.facebook.com/CogitoConsultoriaEstatistica Curriculum: http://lattes.cnpq.br/7541377569511492

1 0

mudar de decimal para inteiro
by Fátima Lima Paula 19 Abr '13

19 Abr '13

Prezados meus números estão aparecendo dessa forma: 3.31000e+12 Como faço para que fique 3310008457687 Alguém pode me ajudar, por favor? Obrigada

2 2