- R-br - DINF UFPR

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 22 Nov '16

22 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

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

20 Nov '16

3 3

Residuos do Teste T de Student
by Graciliano Galdino 19 Nov '16

19 Nov '16

Olá a todos, Quando faço uma análise qualquer usando teste T (t.test) não aparecem os resíduos na lista dos resultados da análise, como aparece na ANOVA. Como faço para extraí-los? Abraço -- Graciliano Galdino A. dos Santos Biólogo Doutorando em Ciências Florestais - PPGCF Universidade Federal Rural da Amazônia - UFRA

4 4

experimento com delineamento diferente
by Ze Henrique 18 Nov '16

18 Nov '16

Boa tarde a todos, gostaria de saber se alguém poderia me ajudar em um problema. Foi feito a coleta de solo em 2 altitudes, dentro de cada altitude em 4 sistemas diferentes (mata, pastagem, lavoura degradada e lavoura conservada) e dentro de cada um dos sistemas em duas diferentes profundidades (0-20 cm e 20-40 cm). Sendo que o que tenho de repetição é apenas três amostras coletadas por profundidade na mesma sistema e altitude. Seria uma parcela sub-subdividida em DIC ou qual o delineamento e alguém conhece algum lugar onde eu possa ter estas informações ou algum pacote no R que faça este tipo de análise. Gostaria de verificar as interações se existe e fazer um teste de média. Desde já agradeço a todos pela contribuição.

2 2

Erro checando o pacote: No Visivle global function
by Elias Carvalho 18 Nov '16

18 Nov '16

Pessoal, Desculpe, acabei enviando sem assunto. Estou construindo um pacote onde uso o pacote Rgraphviz do site " http://bioconductor.org/biocLite.R", porém quando vou checar o pacote gera o seguinte alerta: check.package: no visible global function definition for ‘biocLite’Undefined global functions or variables: biocLite Alguém sabe como resolver isso ? Embora seja um Warning eu gostaria de resolver -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

2 1