Abril 2015 - R-br - DINF UFPR

Re: [R-br] Regex em uma string de busca
by Éder Comunello 16 Abr '15

16 Abr '15

Talvez algo assim... ### <code r> { set.seed(123) df <- data.frame( municipio=sample(c("Limeira", "Manaus", "Dourados"), 30, repl=T), diasemanaincident=sample(c("terca", "quarta", "sexta"), 30, repl=T), trechoklmetros=sample(0:400000, 30), turnoincidente=sample(c("manha", "tarde"), 30, repl=T), stringsAsFactors=F) } which(df$municipio == 'Limeira') # [1] 1 6 15 17 18 19 29 30 which(df$municipio == 'Limeira' & df$diasemanaincident == 'terca') # [1] 15 17 19 vec <- c("df$municipio == 'Limeira'", "df$municipio == 'Limeira' & df$diasemanaincident == 'terca'"); cbind(vec) for (a in vec) print(with(df, which(eval(parse(text=a))))) # [1] 1 6 15 17 18 19 29 30 # [1] 15 17 19 ### </code> Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W] Em 16 de abril de 2015 14:32, Fernando Gama <f.fabiogama88(a)gmail.com> escreveu: > Obrigado pela resposta Eder. :) > > mas os dados não estão em um dataframe e sim eu tenho um vetor de > caracteres na forma de: > > [1] "municipio == Limeira" > > [2] "municipio == Limeira & diasemanaincident == terca" > > [3] "municipio == Limeira & diasemanaincident == terca & *trechoklmetros > == ]12.200, 245.000]*" [1] [1] "municipio == Manaus" > > > [4]"municipio == Manaus & diasemanaincident == sexta" > > [5] "municipio == Manaus & diasemanaincident == sexta & turnoincidente == > manha" > > Meu objetivo é ler essas strings de busca de maneira automática. Como elas > seguem um padrão, posso aplicar o mesmo script em cima de proximas strings > de busca que surgirem... > > > >

1 0

Regex em uma string de busca
by Fernando Gama 16 Abr '15

16 Abr '15

Bom dia amigos, estou tentando criar um padrão com as strings de busca do meu dataset. Essa é minha string: "municipio == Limeira" "municipio == Limeira & diasemanaincident == terca" "municipio == Limeira & diasemanaincident == terca & *trechoklmetros == ]12.200, 245.000]*" "municipio == Manaus" "municipio == Manaus & diasemanaincident == sexta" "municipio == Manaus & diasemanaincident == sexta & turnoincidente == manha" o resultado que pretendo alcançar é: municipio == ' Limeira ' municipio == ' Limeira ' & diasemanaincident == ' terca ' municipio == ' Limeira ' & diasemanaincident == ' terca ' & ** *trechoklmetros > 12.200 & * *trechoklmetros <= **245.000* ) municipio == 'Manaus' municipio == 'Manaus' & diasemanaincident == 'sexta' municipio == 'Manaus' & diasemanaincident == 'sexta' & turnoincidente =='manha' Ou seja, em termos gerais, inserir *aspas simples (ou duplas*) depois do sinal de atribuição == e a segunda condicao é *ajustar o intervalo ]12.200, 245.000]* (em negrito). Tentei gsub, sub, package stringr, porém sem sucesso, nao sou muito familiarizado com regex :/ Quem puder me ajudar... -- Att, Fernando Gama da Mata

2 2

Importar arquivos de imagem para o R
by Pedro Emmanuel Alvarenga Americano do Brasil 15 Abr '15

15 Abr '15

Amigos R, Eu estou fazendo uma analise com um pacote cujas funçoes graficas utilizam o pacote grid. O que eu gostaria de fazer é de salvar dois graficos no mesmo arquivo lado a lado. Eu não tenho familiaridade com grid então estou batendo cabeça. POstei uma mensagem há dois dias e ninguem respondeu então acho que estou num beco sem saida. Me ocorreu que os amigos que trabalham com mapas podem ter a solução. Seria possivel salvar os graficos separadamente, e depois importar ambos os graficos em diferentes janelas do mesmo dispositivo? Se isso for possível, alguem pode me indicar algum exemplo que eu pudesse seguir? Abraço forte, Pedro Brasil

2 1

Re: [R-br] [Regressao logística] - comparacao de coeficientes beta
by Tiago Souza Marçal 15 Abr '15

15 Abr '15

Acho que você deseja fazer a analise de identidade de modelos. O procedimento foi descrito por Regazzi (2004, 2010) para o delineamento inteiramente ao acaso e blocos ao acaso respectivamente. Caso tenha interesse nesta metodologia em R é só entrar em contato. Att. Tiago. Manuella Lech Cantuaria <manuellalech(a)gmail.com> escreveu: Olá a todos, Alguém de vcs conhece algum teste estatístico para comparacao de coeficientes beta obtidos por regressao logistica? Estou analisando diferentes modelos nos quais vario uma das variáveis independentes e minhas amostras sao dependentes. Desde já, muito obrigada :-) Manuella _______________________________________________ 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 2

Análise de Fatores PCA e Dendograma
by Andre Oliveira 15 Abr '15

15 Abr '15

Pessoal boa tarde,estou com dois bancos de dados com a estrutura do USArrests que vem no R só que repetido no tempo. USArrests2007, USArrests2009, USArrests2010. Gostaria de orientação de como anexar estes dados para em um banco de dados único e rodar PCA, Plotar, dendograma e a análise de fatores. Pelo que percebi em meus devo fazer o seguinte. Montar um banco de dados com os três bancos e tirar a soma ou média. Não sei se fere princípios da análise. dados= data.frame(USArrests2007, USArrests2009, USArrests2010) # Não sei se é data.frame aqui mesmo! V1=tapply(V1,grupos,mean) .....Vn=tapply(Vn,grupos,mean) dadosmedio=cbind(V1, ........................,Vn)hclust(dist(dadosmedio)prcomp(dadosmédio) Ou usar a estrutura com Vn=tapply(Vn,grupos,sum) em vez de Vn=tapply(Vn,grupos,mean). Ou usar o próprio banco formado por USArrests2007, USArrests2009, USArrests2010. Caso falei besteira, desconsidere. Ficarei agradecido por sugestões. obrigado pela atenção 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

5 11

Coluna categorizando variável por intervalo
by Fátima Lima Paula 15 Abr '15

15 Abr '15

Prezados, tenho uma coluna com variáveis com uma letra e três números, conforme exemplo abaixo. a<-c("B095","B096","B784","B561","N590","N980","N067","C678","C234","C009","B095", "B561","N980","N980","C234","C009") b<-c("Maria","Paola","José","Pedro","Kleber","Henrique","Marcelo","João","Lia", "Ingrid","Renato","Paulo","Mariana","Mario","Luiza","Luiz") banco<-data.frame(a,b) banco Quero criar uma coluna de forma que as linhas que tiverem de C000 a C300 recebam a palavra"verde", as que tiverem de C301 a C700, recebam "amarelo". As que tiverem de B000 a B550 recebam vermelho e de B 551 a B999 recebam "branco". O restante ficaria NA. Alguém poderia me ajudar? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

3 4

[Regressao logística] - comparacao de coeficientes beta
by Manuella Lech Cantuaria 14 Abr '15

14 Abr '15

Olá a todos, Alguém de vcs conhece algum teste estatístico para comparacao de coeficientes beta obtidos por regressao logistica? Estou analisando diferentes modelos nos quais vario uma das variáveis independentes e minhas amostras sao dependentes. Desde já, muito obrigada :-) Manuella

1 0

Re: [R-br] Gráfico/CL50
by Tiago Souza Marçal 13 Abr '15

13 Abr '15

Use a função points() Att. Tiago. Flavia Batista Gomes - Embrapa Amazonia Ocidental - CPAA <flavia.b.gomes(a)embrapa.br> escreveu: Olá, Tenho um gráfico dose-resposta e preciso marcar um ponto específico, o ponto que representa a CL50. Como insiro essa informação no gráfico? Obrigada, Flávia ----- Mensagem original ----- De: "walmes ." <walmeszeviani(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Segunda-feira, 19 de janeiro de 2015 22:49:32 Assunto: Re: [R-br] Probit/CL50 O pacote MASS tem funções para cálculo de dl50. Como a DL50 é uma função (não linear) dos parâmetros do preditor linear de um modelo glm, sua estimação intervalar geralmente se dá por meio da aplicação do método delta. O pacote car pode ser considerado par isso. Se não estiver enganado, a função chama se deltaMethod(). Consulte a documentação para exemplos de uso. À disposição. Walmes. -- ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== _______________________________________________ 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. ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. Confidentiality note This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure. _______________________________________________ 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 1

Residuos quantile randomizados (qresiduals{stamod})
by Fernando Antonio de souza 13 Abr '15

13 Abr '15

Alguém pode me explicar os residuo quantile aleatorizados? Como posso usá-lo para avaliar ajuste do modelo de distibuição de meu modelo através de quantil-quantil plot, uma vez que o mesmo utiliza aleatorizações para gerar residuos de distribuição contínua, ou seja a cada rodada ele dará valores diferentes. Tem sentido utilizá-lo com a função shapiro.test()? -- ======================================================================= Fernando Souza Zootecnista, DSc. Produção Animal celular: (+55) 82 8113-8781 e-mail:nandodesouza@gmail.com https://producaoanimalcomr.wordpress.com/ ========================================================================

2 3

Sobre R Markdow
by Marcio Leles 12 Abr '15

12 Abr '15

Boa tarde. Estou fazendo um script referente a uma prova e queria que somente saísse os resultados e não os comandos. Estou usando o Rmarkdow no Rstudio. Gostaria de saber se alguém sabe o comando para que ficasse uma questão a cada página, ou seja o comando para quebrar página. Como abaixo eu gostaria que a quatro 6 ficasse em uma página e a 8 na outra. O pedaço do script é. **Resolução Para a Questão 6** ```{r echo=FALSE} dic(DIC1$trat, DIC1$htponto, quali = TRUE, mcomp = "tukey", sigT = 0.05, sigF = 0.05) ``` ```{r include=FALSE} DBC1<-read.table("/Users/Marcio/Documents/_R/rdbc_prova1_2015.txt",header=TRUE, dec = ",") ``` **Resolução Para a Questão 8** ```{r echo=FALSE} dbc(DBC1$trat, DBC1$bloco,DBC1$tabponto, quali = TRUE, mcomp = "tukey", sigT = 0.05, sigF = 0.05) ``` -- Atenciosamente; Prof. Marcio Leles R. de Oliveira http://marcioromarco.webnode.com/ Mensuração e Manejo Florestal Universidade Federal dos Vales do Jequitinhonha e Mucuri - UFVJM Departamento de Engenharia Florestal Campus JK - Rodovia MGT367-KM 583, 5000 - Alto da Jacuba Diamantina - MG, CEP: 39100-000 Tel: (38) 9158-1435

5 4