Novembro 2016 - R-br

[OFF -Tópic] Cotas em concursos para negros e índios!
by Andre Oliveira 12 Jun '49

12 Jun '49

Pessoal, tem um colegar meu que vou fazer um concurso ele é descendente de negro ele é pardo. Ele veio com a pergunta como provar que sou negro, pois, tem vagas reservada para negros e índios? Eu sinceramente não sei e não tenho ideia pois nunca tive esta curiosidade. O que de fato comprova sua raça? Algum documento? O que são usados em concursos para provar? obrigado André Oliveira Souza

5 6

Legenda eixo X
by Luciano F. Sgarbi 10 Jun '49

10 Jun '49

Prezados, Estou tentando centralizar a legenda do eixo X em uma figura com vários gráficos. Estou conseguindo centralizar a legenda da linha inferior como no exemplo, mas as linhas do meio não. par(mfrow=c(2,2)) plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 1",outer=TRUE) ### o problema está aqui, pois a legenda desce para a parte inferior da figura. plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 2",outer=TRUE)## Esse dá certo Alguém podia dar uma luz de como eu posso fazer isso? Abraços -- *Luciano F. Sgarbi* Mestrando em Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás, campus II Goiânia-GO 74001-970 Brazil

2 2

Cálculo do coeficiente de assimetria (skewness coefficient) - ecologia de populações
by Juliana Loureiro 01 Dez '16

01 Dez '16

Olá pessoal! Estou com uma dúvida relacionada ao cálculo do coeficiente de assimetria, que ajuda a prever a direção de mudança de uma população (crescimento ou declínio). No caso, eu trabalho com a população de uma espécie de palmeira e dividi os indivíduos em classes de altura. Gostaria de saber como faço pra calcular esse coeficiente utilizando o R. Agradeço muito, Juliana Loureiro

4 3

Dúvida sobre o gls
by Luiz Leal 01 Dez '16

01 Dez '16

Boa dia a todos, tenho algumas dúvidas sobre os mínimos quadrados generalizados (gls) no que se refere aos pressupostos de normalidade. Quando utilizo gls(y ~ x, ...) pressupõe-se que y tem distribuição normal? E os resíduos também? Ou deve me preocupar apenas em analisar o gráficos dos resíduos (bem distribuídos em torno de zero e sem padrão)?. A grande totalidade dos dados que trabalho quase nunca y segue uma distribuição normal. A outra dúvida se refere ao cálculo no R. Conforme script abaixo consegui reproduzir a fórmula de calculo da função gls() nos casos em que weights = varFixed() e weights = varFunc(), mas não entendi a fórmula de cálculo nos casos em que weights = varExp() e weights = varPower(). Li a documentação do R, mas não consegui reproduzir. Alguém sabe?Desde já agradeço.Luiz # Dados x=c(1,2,3,3,4,5,6,7,8,9) y=c(7,5,7,8,8,5,7,9,9,9) # Carregar pacote library(nlme) # Matrix X X=model.matrix(~x) # varFixed / varFunc r=fitted(lm(y~x)) omegafit=matrix(0,length(r),length(r)) ; diag(omegafit)=r solve( t(X) %*% solve(omegafit) %*% X ) %*% t(X) %*% solve(omegafit) %*% y gls(y ~ x, weights = varFixed(~r), method = "ML")$coefficients gls(y ~ x, weights = varFunc(~r) , method = "ML")$coefficients # Dúvida # varExp() r=exp( 2*fitted(lm(y~x)) ) omegafit=matrix(0,length(r),length(r)) ; diag(omegafit)=r solve( t(X) %*% solve(omegafit) %*% X ) %*% t(X) %*% solve(omegafit) %*% y gls(y ~ x, weights = varExp(), method = "ML")$coefficients # varPower r=abs(fitted(lm(y~x)))^2 omegafit=matrix(0,length(r),length(r)) ; diag(omegafit)=r solve( t(X) %*% solve(omegafit) %*% X ) %*% t(X) %*% solve(omegafit) %*% y gls(y ~ x, weights = varPower(), method = "ML")$coefficients

3 4

Extraindo parte de uma formula
by Pedro Emmanuel Alvarenga Americano do Brasil 30 Nov '16

30 Nov '16

Amigos de R, Estou tentando fazer uma função com alguns loops que comparam modelos em diferentes subsets de dados. No caminho, eu precisaria extrair a parte dependente das formulas dos modelos, mas não estou sabendo como. O que me interessa é repetir a parte dependente da fomrula nos diferente subsets. > library(survival) > f <- coxph(Surv(time, status) ~ ph.ecog + age, data=lung) > f$formula Surv(time, status) ~ ph.ecog + age # Eu gostaria de extrair o "Surv(time, status)" para um objeto de tal forma que eu consiga seguir um modelo igual ao abaixo. objeto <- Surv(time, status) # Esse passo que eu não estou conseguindo. E deveria vir do primeiro fit "f" f2 <- survfit(objeto ~ 1, data=lung, subset = lung$age >= 63 , type = "kaplan-meier") f2b <- survfit(objeto~ 1, data=lung, subset = lung$age < 63 , type = "kaplan-meier") Abraço forte, Pedro Brasil

3 3

logica condicional e criar arquivo netCDF
by Mateus Dias Nunes 25 Nov '16

25 Nov '16

Olá, Eu calculei a média e o percentil 95% (P95), da área que desejo trabalhar, porém preciso fazer duas coisas: 1 - fazer a média somente dos casos acima do P95, ou seja, não consegui achar a lógica para que mostre os valores numéricos onde casos <- UV[UV>=resultado] e fazer a média somente destes casos acima do P95 para cada ponto de grade. 2- criar um arquivo .nc que contenha os valores de "casos". Abaixo o link do meu script https://www.dropbox.com/s/ph03e1b9l6ziylg/TESTE.R?dl=0 E aqui abaixo o link do dado .nc que estou utilizando https://www.dropbox.com/s/81k73nffpt1zwq7/310times_fillmissJAN.nc?dl=0 Desde já agradeço, ____________________________________________________________________________ MATEUS DIAS NUNES MESTRANDO DO PROGRAMA DE PÓS-GRADUAÇÃO EM METEOROLOGIA - PPGMET GRUPO DE ESTUDOS EM POLUIÇÃO ATMOSFÉRICA - GEPA UNIVERSIDADE FEDERAL DE PELOTAS - UFPEL TELEFONE: +55 (53) 98112-5154 SKYPE: mateus.diasnunes ____________________________________________________________________________ Sent with Mailtrack <https://mailtrack.io/install?source=signature&lang=en&referral=nunes.mateus…>

2 1

Função escondida - como achar
by Pedro Emmanuel Alvarenga Americano do Brasil 25 Nov '16

25 Nov '16

Amigos de R, Instalei o pacote SurvCorr que faz quase o que eu quero. Entao eu digitei no console >survcorr E a função aparece no console. Eu a copiei para um script e fiz pequenas modificações pra poder se adequar a minha necessidade e renomiei a função para survcorr2 e salvei num arquivo. Quando fui testar percebi que há no final do script uma função chamada "pearson". Inicialmente eu achei que era uma função do pacote original. Então eu carreguei o pacote achando que a survcorr2 iria chamar a pearson do pacote survcorr, mas vem o aviso... >survcorr2(formula1 = Surv(los, desfecho.los) ~ 1, obj = y, obj2 = x2, data = bv) Error in survcorr2(formula1 = Surv(los, desfecho.los) ~ 1, obj = y, obj2 = x2, : could not find function "pearson" Procurei na documentação e não achei em qualquer lugar essa função. Mesmo editando o script para onde essa função aparece como SurvCorr::pearson pra ter certeza que ele procura onde devia, não funciona. Ela tem que estar la em algum lugar porque o exemplo da survcorr original funciona muito bem. Eu já tinha me deparado com o problema de saber que a função existe e não conseguir olhar dentro da função, mas nunca tinha me deparado com não achar a função. Alguma luz em como achar essa função? Pedro Brasil

3 2

sobre a utilização do paste para formar nomes de objetos
by Antonio Silva 24 Nov '16

24 Nov '16

Olá a todos Quero fazer cálculos de íncides que utilizam duas matrizes. No total tenho mais que 100. matriz1 <- matrix(.... matriz2 <- matrix(.... matriz3 <- matrix(.... ls() matriz1 matriz2 matriz3 Tentei fazer um looping for (a in 1 to ...) { for (b in 1 to ...) { resp <- função(paste("matriz",a,sep=""),paste("matriz",b,sep="")) assing(paste("resp_",a,"_",b,sep=""),resp) } } No entanto as funções não reconhecem o resultado de paste("matriz",a,sep="") ou as.name(paste("matriz",a,sep="")) como o objeto matriz1, matriz2, etc. Alguém sabe uma forma de "montar" o nome de um objeto de forma que este seja reconhecido como tal? Obrigado Antônio Olinto

3 3

Re: [R-br] Heterocedasticidade e teste de média
by sznelwar＠uol.com.br 23 Nov '16

23 Nov '16

2 1

Re: [R-br] Loop para ler e trabalhar com diversos arquivos simultaneamente
by sznelwar＠uol.com.br 22 Nov '16

22 Nov '16

1 0