- R-br - DINF UFPR

Geração de números aleatórios com amplitude determinada
by Marcos Bissoli 29 Jul '14

29 Jul '14

Boa noite a todos, Preciso gerar uma sequência de números aleatórios, de uma distribuição normal. No entanto, a função <rnorm> não permite definir limites mínimo e máximo para o vetor de dados. Assim, criei uma função para tal fim: >Gerador.Aleatorio <- function(n,menor,maior,x,s) +{ + dados <- rnorm(n,x,s) + while ((min(dados)<menor)|(max(dados)>maior)) + { + dados <- rnorm(n,x,s) + } + dados <- dados + print(dados) +} Onde: - n é o número de valores que quero gerar - menor é o limite mínimo do vetor - maior é o limite máximo do vetor - x é a média do vetor - s é o desvio padrão do vetor O problema é que, como usei um laço ("loop"), o custo de tempo computacional dessa função tornou-a simplesmente impraticável. Até 100 valores eu consegui gerar. Mas agora tentei com 284 e já jantei, fumei, tratei de assuntos familiares e... Nada! E minha demanda, na verdade, é para gerar mais de 1500 valores, simulando estudo epidemiológico observacional, a partir de valores já descritos em um estudo piloto que realizei (ou seja, eu tenho já estimativas de n, mínimo, máximo, média e desvio padrão). Há alguma função, pacote ou sintaxe que torne isso possível? E qual seriam esses "limites do custo computacional"? Grato desde já por qualquer ajuda. Abraços fraternos, -- MARCOS BISSOLI Faculdade de Nutrição Universidade Federal de Alfenas E-mail: mbissoli(a)gmail.com Twitter: #mbissoli Facebook: https://www.facebook.com/MarcosBissoli Alfenas, Minas Gerais, Brasil *****Pense na Natureza antes de Imprimir***** Divulgue ON-LINE Eu apoio a ENEN "na luta por um Brasil sem fome" "por ĉiu popolo ties propran lingvon, por ĉiuj popoloj la esperantan" (para cada povo sua própria língua, para todos os povos o Esperanto) E nunca votarei no PSDB/DEM!

3 2

Re: [R-br] Inclusão de variáveis numa regressao
by Robert Iquiapaza 28 Jul '14

28 Jul '14

Tente algo similar vars=colnames(seubancodedados) fmla <- as.formula(paste("y ~ ",paste(vars[grep("iso1",vars)],collapse="+"))) lm(fmla,data=seubancodedados) Sds Em 28 de julho de 2014 16:39, Paulo Justiniano <paulojus(a)leg.ufpr.br> escreveu: > > pq ovce nao cria um data-frame apenas com as variáveis que quer usar? > > > > > > On Mon, 28 Jul 2014, Alexandre Loures wrote: > > Obrigado Paulo! >> >> Mas meu data-frame tem outras colunas além de y e iso1. >> >> >> >> >> Em 28/07/2014 15:38, Paulo Justiniano escreveu: >> se o seu data-frame só tiver as colunas y e iso* >> baste fazer: >> >> lm(y ~. , ...) >> >> ond o "." significa todas as demais colunas do data-frame exceto o >> y >> >> >> >> On Mon, 28 Jul 2014, Alexandre Loures wrote: >> >> Boa tarde! >> >> Preciso rodar uma regressão em que há muitas variáveis cujo >> nome inicia da mesma forma. Por exemplo: >> >> iso1ARG >> iso1BRA >> iso1PRY >> iso1URY >> >> Alguém sabe de alguma forma de incluir essas variáveis de uma >> só vez? >> >> Ou seja, ao invés de utilizar o padrão >> >> reg <- lm (y ~ iso1ARG + iso1BRA + iso1PRY + iso1URY + ...) >> >> Haveria alguma maneira em que definindo apenas iso1 todas as >> variáveis que inicia assim seriam incluídas? >> >> Por exemplo: >> >> reg <- lm (y ~ iso1*) >> >> Outra dúvida: nos resultados há como definir para que seja >> mostrado apenas algumas variáveis. >> >> Suponha as seguintes variáveis: rta fta bothin gsp sur1995 >> mercef >> >> Mas quero que apareça apenas os resultados para sur1995 e >> mercef. Há como fazer? >> >> >> Desde já muito obrigado! >> >> >> -- >> Alexandre Rodrigues Loures >> >> >> >> >> _______________________________________________ >> 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. >> >> >> -- >> 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 >> Tel.: +55 83 9804-1000 >> e-Mail: alexandre.loures(a)ymail.com >> Site: www.ccsa.ufpb.br/ppge >> >> >> > _______________________________________________ > 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

Como fazer predição para pontos não amostrados?
by Jacqueline Cantú 28 Jul '14

28 Jul '14

Olá pessoal..Tenho outra dúvida..Construí o semivariograma cruzado, ajustei um modelo e quero fazer a predição para pontos que não foram amostrados para isso construí uma malha fina...Minha dúvida é: como fazer a predição espacial para os pontos não amostrados?Lembrando que são 2 variáveis, sendo que para predizer os valores de uma variável será utilizado os valores da outra variável (cokrigagem) Segue rotina.. require (geoR)library(sp)require(gstat)require(splancs)require(classInt) #lendo os dadosh <- read.table("phca.txt ",head=T)names(h)class(h)coordinates(h) = ~x+y g = gstat(NULL, "ph", ph~1, h)g = gstat(g, "ca", ca~1, h)plot(variogram(g))vm = variogram(g)variogram(g) ##Ajuste do modelo esfericosph <- fit.lmc(vm,g, vgm(0.3, "Sph", 200, 0.2))sph #modelo ajustado plot(vm,sph) ##Predição para pontos não amostrados#construção de uma malha finaborda <- read.table("borda.txt") borda plot(borda)polygon(borda)apply(borda,2,range) gr1 <- expand.grid(x=seq(239182.8, 241027.9, by=5),y=seq(7236744, 7238593, by=5))gi <- polygrid(gr1, bor=borda)plot(gi) Att, Jacqueline

2 1

Como fazer mapa temático
by Jacqueline Cantú 26 Jul '14

26 Jul '14

Olá pessoal,Construí o semivariograma cruzado, ajustei modelo e escolhi o melhor modelo e preciso fazer o mapa temático para cada uma das variáveis. Alguém sabe me dizer como fazer o mapa temático? Muio obrigada pela atenção.Att, Jacqueline

5 6

Pular amostras com singularidade - optim
by jborssoi . 24 Jul '14

24 Jul '14

Caros, gerei J amostras e estou fazendo a estimação dos parâmetros do meu modelo com optim, sendo que a cada nova estimação o chute inicial é atualizado. Porém, para algumas amostras está dando problema de singularidade e interrompe o processo. Gostaria de saber se tem uma forma de pular as amostras que dão o problema? Grato pela ajuda! Joelmir

3 2

Ajuda Urgente - Pooled Incidence Rate
by Ary Serpa Neto 24 Jul '14

24 Jul '14

Olá amigos, Estou precisando construir um forest plot (similar a imagem em anexo) usando "pooled incidence rate". Alguem pode me dar uma luz de como fazer isso no R? Desde ja agradeço a atenção Grande abraço Ary Serpa Neto, MD, MSc

3 2

Teste para comparar distribuições
by Luiz Roberto Martins Pinto 23 Jul '14

23 Jul '14

Caros amigos, Fiz uma pesquisa no Google, mas não consegui uma solução satisfatória. Então peço ajuda. Preciso de encontrar um teste que *quantifique* a semelhança entre dois vetores para checar se eles têm o mesmo tipo de distribuição, por meio de p.value. O interesse principal é comparar os vetores x e ab. Data set: http://www.datafilehost.com/d/f3757310 library(Hmisc) library(grDevices) load('NNC D_d3kk2.RData') # Avaliação qualitativa (gráfica): dev.set(which=1) bpplot(d3) boxplot(d3) # Aparentemente os vetores x, abC e ab são muitos semelhantes # A minha expectativa é encontrar um teste que quantifique esta semelhança por meio de p.value # Então fiz os testes abaixo # Avaliação quantitativa: v2=d3$ab # Comparando o vetor x com o vetor ab: v1=d3$x round(ks.test(v1, v2)$p.value,2) # Out put: p.value=0. Então o vetor x é estatisticamente DIFERENTE do vetor ab. # Comparando o vetor abC com o vetor ab: v1=d3$abC round(ks.test(v1, v2)$p.value,2) # Out put: p.value=0.65. Então o vetor x é estatisticamente IGUAL do vetor abC. # Comparando o vetor abC com o vetor x: v2=d3$x round(ks.test(v1, v2)$p.value,2) # Out put: p.value=0. Então o vetor abC é estatisticamente DIFERENTE do vetor x. ### *Testes alternativos*, para comparar o vetor x com ab: v2=d3$ab v1=d3$x t1=round(ks.test(v1, v2)$p.value,2);t1 # Out put: p.value=0 t2=round(ks.test(v1, v2, alternative = "l")$p.value,2);t2 # Out put: p.value=0 t3=round(ks.test(v1, v2, alternative = "g")$p.value,2);t3 # Out put: p.value=0.85 # Então, para este Data set (NNC D_d3kk2.RData), e para o t3=0.85. Então, por meio deste teste, o vetor x é estatisticamente IGUAL do vetor ab. # *No entanto*, fazendo os testes com conjuntos de dados diferentes (10000 Data set) (mas, sempre com mesmo tamanho de vetor), quando comparo o vetor x com o vetor ab: ## em aproximadamente 60% das comparações (ou em 6000 Data set) t1>0.05 ## em aproximadamente 80% das comparações (ou em 8000 Data set) t2>0.05 ## em aproximadamente 80% das comparações (ou em 8000 Data set) t3>0.05 ## em aproximadamente 85% das comparações (ou em 8500 Data set) t1 *OU* t2 >0.05 ## em aproximadamente 85% das comparações (ou em 8500 Data set) t1 *OU* t3 >0.05 ## em 0% das comparações (ou em 0 Data set) t2 *E* t3 >0.05 ## em 100% das comparações (ou em 10000 Data set) t2 *OU *t3 >0.05 Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC Laboratório de Estatística Computacional Universidade Estadual de Santa Cruz Ilhéus-Bahia luizroberto.uesc(a)gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831

3 4

Re: [R-br] Como fazer mapa temático
by Waldir De Carvalho Junior 22 Jul '14

22 Jul '14

Prezada Jacqueline sua questão está muito genérica. É preciso ser mais específico para te ajudar. No entanto, os pacotes 'sp', 'rgdal', 'plotkml', entre outros, tem a possibilidade de plotar os mapas. verifique Att __________ Dr. Waldir de Carvalho Junior Engº Agroº D.Sc. Ciência do Solo - Mapeamento Digital Embrapa Solos Empresa Brasileira de Pesquisa Agropecuária (Embrapa) Rio de Janeiro/RJ Waldir.Carvalho(a)embrapa.br Telefone: +55 (21) 2179-4614 +55 (21) 2179-4614 | Fax: +55 (21) 2274-5291 www.embrapa.br/solos | twitter.com/embrapa Confira também: www.facebook.com/agrosustentavel Call Send SMS Add to Skype You'll need Skype Credit Free via Skype ----- Mensagem original ----- De: r-br-request(a)listas.c3sl.ufpr.br Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Terça-feira, 22 de julho de 2014 12:00:02 Assunto: Digest R-br, volume 43, assunto 22 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. Como fazer mapa temático (Jacqueline Cantú) ---------------------------------------------------------------------- Message: 1 Date: Tue, 22 Jul 2014 17:16:55 +0300 From: Jacqueline Cantú <jacqueline.cantu(a)hotmail.com> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Como fazer mapa temático Message-ID: <BLU177-W17CD82038F1F8B202CBD5A8BF10(a)phx.gbl> Content-Type: text/plain; charset="iso-8859-1" Olá pessoal,Construí o semivariograma cruzado, ajustei modelo e escolhi o melhor modelo e preciso fazer o mapa temático para cada uma das variáveis. Alguém sabe me dizer como fazer o mapa temático? Muio obrigada pela atenção.Att, Jacqueline -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140722/ad115a93/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 43, assunto 22 ***************************************** ____________________________________________________________________________ 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.

1 0

Identificar pontos num gráfico e adicionar rótulos
by Helder Gramacho 21 Jul '14

21 Jul '14

Bom dia pessoal, Estou tentando usar o identify() para identificar pontos num gráfico e em seguida adicionar um rótulo para os pontos identificados, este assunto já havia sido discutido aqui na lista e foi resolvido com o comando abaixo, no meu caso estou utilizando o RStudio, só que depois que executo o comando parece que o gráfico trava e os rótulos não são adicionados, está faltando algo nos comandos? ### <code r> dados <- data.frame(x=1:10,y=1:10,Nome=paste('ponto',letters[1:10])) plot(dados$x,dados$y) identify(dados$x,dados$y, labels=dados$Nome) ### </code> Obrigado, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

3 4

Selecionar vários níveis dentro de um fator com for()
by ASANTOS 16 Jul '14

16 Jul '14

Boa noite pessoal, Tenho um fator (indiv) no objeto dados_cal[,3], composto de vários níveis sendo soldado1, .... ate soldado 20 e operario1, ... ate operario20, gostaria de selecionar só os operário, para isto eu fiz: #Operários x<- c("operario1","operario2","operario3","operario4","operario5","operario6", "operario7" ,"operario8","operario9","operario10","operario11","operario12", "operario13","operario 14","operario15","operario16","operario17","operario18", "operario19","operario20") RES=NULL for(l in 1:length(x)){# Rayon des cercles teste<-dados_cal[dados_cal[,3]==x[l]] RES=rbind(RES,teste) } # Mas o meu for() aplicado ao objeto x não funcionou, alguém poderia dar um help? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

2 2