Julho 2014 - R-br - DINF UFPR

análise de agrupamento: índice de similaridade
by Aline Lipsky 30 Jul '14

30 Jul '14

Boa tarde, Estou fazendo análise de agrupamento e utilizei nos meus dados o pacote vegan, que faz análise de dissimilaridade. Gostaria de saber se tem um pacote no R que faz análise de agrupamento com matriz de similaridade, pois estou aplicando o índice de similaridade de Jaccard. Grata, Aline

3 3

Como recarregar uma função de pacote alterada?
by Thiago V. dos Santos 30 Jul '14

30 Jul '14

Pessoal, Mudei alguns parâmetros de uma função de um determinado pacote no R, basicamente cores dos rótulos dos eixos de um gráfico (que eram vermelhos e que preciso que fiquem pretos). Para fazer isso eu fui no diretório "scripts" do pacote e editei três linhas do código .R que continha a função. Fechei o RStudio, carreguei o pacote novamente e testei a função alterada, que para minha surpresa continuou produzindo os mesmos resultados de antes de alteração. Preciso reinstalar o pacote com a função alterada? Tem alguma outra coisa que deve ser feita e que estou pulando? Obrigado por qualquer dica, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/ind… Phone: (612) 323 9898

3 2

Re: [R-br] Geração de números aleatórios com amplitude determinada
by Ivan Bezerra Allaman 29 Jul '14

29 Jul '14

Você quer usar uma normal truncada. Veja o pacote "truncnorm". Cordialmente,

5 4

Como fazer cokrigagem para pontos não amostrados
by Jacqueline Cantú 29 Jul '14

29 Jul '14

Olá pessoal,Gostaria de saber como fazer cokrigagem para pontos não amostrais..Alguém sabe me dizer como faz? Att, Jacqueline

2 2

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 29 Jul '14

29 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