Dezembro 2011 - R-br

soma condicional
by Luiz André Zardo 19 Dez '11

19 Dez '11

Como posso obter uma soma condicional no R?

4 4

Fwd: Duvida para an[alise
by Adriano Carvalho Costa 19 Dez '11

19 Dez '11

---------- Mensagem encaminhada ---------- De: Adriano Carvalho Costa <acarvalhocosta(a)gmail.com> Data: 19 de dezembro de 2011 13:03 Assunto: Duvida para an[alise Para: r-br-request(a)listas.c3sl.ufpr.br Boa tarde, Estou com uma duvida para analisar este experimento. Alguem pode me ajudar. Penso da seguinte forma. Contrastes 1. Pectina e Amido vs Prop.Ca e CaCl 2. Pectina vs Amido 3. Pectina 1,5 vs Pectina 7,5 e Pectina 15 4. Pectina 7,5 vs Pectina 15 5. Amido 15 vs amido 45 e amido 75 5. Amido 45 vs amido 75 6. Prop. Ca vs CaCl 7. Prop. Ca 0 vc Prop. Ca 15 e Prop. Ca 30 8. Prop. Ca 15 vs Prop. Ca 30 9. CaCl 0 vc CaCl 15 e CaCl30 10.CaCl 15 vs CaCl Ca 30 tratamento Concentrações Pectina 1,5 Pectina 7,5 Pectina 15 amido 15 amido 45 amido 75 Prop. Ca 0 Prop. Ca 15 Prop. Ca 30 CaCl 0 CaCl 15 CaCl 30

2 5

soma acumulada
by Luiz André Zardo 19 Dez '11

19 Dez '11

Prezados colegas, Como faço para criar uma variável em um dataframe que seja soma acumulada de outra variável? i Xi Yi 1 a a 2 b a+b 3 c a+b+c Obrigado.

2 1

H0 é a hipótese nula?
by Pedro Emmanuel Alvarenga Americano do Brasil 19 Dez '11

19 Dez '11

Amigos de R, Um pequeno divertimento para desejar boas festas aos amigos virtuais que ajudam em horas oportunas! http://www.phdcomics.com/comics.php?n=1457 Abraço forte e que a força esteja com vocês, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

1 0

Medidas repetidas na mesma unidade experimental!!
by Ivan Bezerra Allaman 19 Dez '11

19 Dez '11

Boa noite senhores! Há alguns dias estou estudando a metodologia de medidas repetidas em uma mesma unidade experimental ao longo do tempo, com o intuito de ajudar uma pesquisadora. Encontrei materiais, livros, artigos, etc que me dão segurança em trabalhar com tal medotologia. No entanto, há uma dúvida que me deixou num impasse, e preciso da ajuda dos colegas mais experientes que trabalham com estes tipos de dados. O objetivo da pesquisadora foi avaliar o uso de dois tratamentos sobre a redução de um tumor cancerígeno em cães do sexo feminino. Inúmeras variáveis (nominais, contínuas, etc.) foram avaliadas (secreção, leucócitos, etc.). As medidas foram tomadas em um mesmo animal ao longo de um período de 30 dias. Até aí tudo bem. Pois posso utilizar a abordagem de modelos mistos considerando o animal como de efeito aleatório, e no modelo incluir uma estrutura de correlação do tipo Arima, sarima, entre outros. Ainda estou avaliando qual melhor estrutura. Mas a dúvida crucial é que além de ter tomadas várias medidas num mesmo animal, este mesmo animal recebeu mais de um tratamento. Por exemplo: Do dia 1 ao dia 15 recebeu o tratamento A e do dia 16 ao dia 30 o tratamento B. Aí que estou encontrando problema, pois creio que há efeitos residuais do tratamento A sobre o B, confundindo os verdadeiros efeitos de cada tratamento. Uma saída, seria tratar os dois tratamentos como apenas um. Alguém já trabalhou com estrutura semelhante? Alguém poderia me recomendar uma literatura para que eu possa ganhar mais este desafio? Desde já, agradeço a atenção de todos! Abraços! (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

8 14

Comportamento estranho nd rStudio com o lme4 + laticce
by Leonard de Assis 19 Dez '11

19 Dez '11

Bom dia, Estava rodando o código R presente no paper: http://www.jstatsoft.org/v20/i02/paper require(lme4) data("lq2002", package = "multilevel") wrk <- lq2002 for (i in 3:16) wrk[[i]] <- ordered(wrk[[i]]) for (i in 17:21) wrk[[i]] <- ordered(5 - wrk[[i]]) lql <- reshape(wrk, varying = list(names(lq2002)[3:21]), v.names = "fivelev", idvar = "subj", timevar = "item", drop = names(lq2002)[c(2, 22:27)], direction = "long") lql$itype <- with(lql, factor(ifelse(item < 12, "Leadership", ifelse(item < 15, "Task Sig.", "Hostility")))) for (i in c(1, 2, 4, 5)) lql[[i]] <- factor(lql[[i]]) lql$dichot <- factor(ifelse(lql$fivelev < 4, 0, 1)) str(lql) summary(lql) (fm1 <- lmer(dichot ~ 0 + itype + (1 | subj) + (1 | COMPID) + (1 | item), lql, binomial)) rr <- ranef(fm1, postVar = TRUE) str(rr) qq <- qqmath(rr) print(qq$subj) O executei no terminal e rodou perfeitamente (R 2.14 dev compilado dos fontes - linux ubuntu) Quando fui executar no RStudio (Mesmo computador, mesmo usuário), encontro a seguinte mens de erro: > print(qq$subj) Error in x$xscale.components(x$x.limits, at = x$x.scales$at, used.at = x$x.used.at, : could not find function "calculateAxisComponents" Isso é normal ou loucura do RStudio? O resultado esperado do comando é o gráfico na pág 12 do Paper (o da esquerda). Isso eu consegui no R (terminal) -- []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com

1 0

[Dúvida] Combinações de strings de um vetor.
by Pedro Rafael 18 Dez '11

18 Dez '11

Pessoal existe algum objeto no R para realizar combinações 2 a 2, 3 a 3. etc? Exemplo: Considere o vetor: vetor <- c("a","b","c","d") Queria obter todas as combinações possíveis com os elementos do vetor, ou seja, uma saída semelhante a: [1] "ab" "ac" "ad" "bc" "bd", "dc" Até poderia pensar em programar isto do zero usando o R mas será que já existe um objeto para fazer isto pronto e eficiente? Obrigado à todos desde já. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

4 11

[Dúvida] Tratamento de Erros no R.
by Pedro Rafael 18 Dez '11

18 Dez '11

Pessoal como tratar erro no R? Veja a função abaixo: dado = "fa\xE7ile" cod = combn(iconvlist(),2) i <- NULL conver <- vector() for(i in 1:dim(cod)[2]){ iconv(dado, from = as.vector(cod[,i])[1], to = as.vector(cod[,i])[2]) } Fiz todas as combinações possíveis de formato de arquivo e gravei no objeto cod. Tenho um outro objeto chamado dado que está em uma codificação até então "estranha", a sting correta seria façile. Queria ver qual codificação correta para converter o objeto dado para o formato correto. Fiz um FOR que usa o comando ivonc para todas as combinações possíveis. Só que a condição para quando ocorre um primeiro erro, ou seja, quando um tipo de conversão não é adequado. Queria saber como desconsiderar estes erros e dizer que pode continuar após a ocorrência de um erro. Mesmo havendo outra forma de resolver meu problema de codificação gostaria de saber como tratar problemas desse tipo, por exemplo dizer que algo pode ser continuado após um erro ou parado quando um erro ocorre ou uma quantidade específica de erros ocorrerem. A outra pergunta diz respeito a como saber a codificação correta do objeto de forma mais eficiente possível. -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

2 2

Usando Sweave e MikTex no Windows [Off-Topic]
by Augusto Ribas 17 Dez '11

17 Dez '11

Ola. Bem eu começei a usar o Sweave no R e é uma excelente ferramenta pra escrever. Mas eu estou apanhando aqui um pouco. Bem, eu estou criando um Arquivo .Rnw num diretorio que concentro todos os meus arquivos referentes a um determinado projeto. Ai uso a Função Sweave() nele para criar o arquivo .tex So que quando vou usar o arquivo .tex com o miktex ele não vai e pede um arquivo, o Sweave.sty. Este que esta em /share/texmf/tex/latex Ai eu tenho que colocar esse arquivo no diretorio senão não da certo compilar(o termo compilar aqui ta certo, ou como eu devo falar ao usar o latex num arquivo tex pra fazer um dvi ou pdf?) o arquivo tex num pdf. Dai a Minha duvida é: Tem como eu dizer no MikTex onde ta esse aquivo Sweave.sty pra ele usar direto o da instalação? E não precisar estar no diretorio que estou trabalhando? Bem é isso. Obrigado :) -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

3 3

Função para retornar uma lista de combinações de colunas de uma matriz
by Emmanuel Arnhold 16 Dez '11

16 Dez '11

Gostaria de alguma dica para operacionalizar a CMR abaixo em uma função. O resultado final é uma lista com combinações das colunas de uma matriz. x=matrix(c(1,2,3,4,5,6,7,8,9), nc=3);x l1=c(x[,1], x[,2]) l2=c(x[,1], x[,3]) l3=c(x[,2], x[,3]) Lista=list(l1,l2,l3);Lista Obrigado. Emmanuel

3 4