Setembro 2011 - R-br

Conexão R com Oracle
by nani.design＠bol.com.br 28 Set '11

28 Set '11

2 1

Ponto de corte em cluster
by Marcelo Cardoso mello 27 Set '11

27 Set '11

Como obter no R um ponto de corte em um cluster (UPGMA ... etc). Algum critério como pseudo test F ou t ou algum outro? Obrigado. M

2 1

Re: [R-br] Gráfico WMF
by Ivan Bezerra Allaman 27 Set '11

27 Set '11

Bom, então um caminho extremamente deselegante seria diminuir a fonte do eixo, o que acha? plot(1:40,runif(40),axes=F,xlab='Eixo X',ylab='Eixo Y') axis(1,at=c(0,seq(2,40,3),40), cex.axis=0.8) axis(2,seq(0,1,.1)) box() savePlot('c:/temp/temp.wmf') dev.off() (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}

9 12

Gráficos 3D para exibir Simulações _
by Paulo Nogueira 27 Set '11

27 Set '11

Olá a todos. Acompanho a lista há bastante tempo e mesmo sem mandar mensagens aprendo muito com as discussões e soluções propostas. Acho que quase tudo que sei de R aprendi por aqui. Sou bastante grato. Minha questão é a seguinte: Estou simulando distribuições sendo que cada uma delas é gerada modificando-se de forma controlada o valor de alguns parâmetros da função geradora. Para compará-las graficamente, exibindo a evolução conforme modificam-se os parâmetros, gostaria de criar um gráfico semelhante ao que coloquei no link abaixo. http://goo.gl/17Erk Não conheço muito de gráficos 3D no R, mas gostaria que vocês indicassem funções/pacotes que eu possa pesquisar e chegar em um resultado semelhante. Toda ajuda é bem vinda. Grato, Paulo Nogueira Starzynski

2 4

Erro: pacote não encontrado
by Pedro Emmanuel Alvarenga Americano do Brasil 27 Set '11

27 Set '11

Amigos de R, O seguinte erro vem ocorrendo diversas vezes comigo com diferentes pacotes... > library('epicalc')Carregando pacotes exigidos: MASSCarregando pacotes exigidos: nnet Anexando pacote: 'epicalc' The following object(s) are masked from 'package:rms': lrtest The following object(s) are masked from 'package:Hmisc': fillin, recode > epicalc::followup.plot(c$prontuario,c$dias,c$CD,by=c$iris)Erro em get(search()[2]) : objeto 'package:epicalc' não encontrado> followup.plot(c$prontuario,c$dias,c$CD,by=c$iris)Erro em get(search()[2]) : objeto 'package:epicalc' não encontrado> search() [1] ".GlobalEnv" "package:epicalc" "package:nnet" "package:MASS" "package:foreign" [6] "package:rms" "package:Hmisc" "package:survival" "package:splines" "tools:rstudio" [11] "package:stats" "package:graphics" "package:grDevices" "package:utils" "package:datasets" [16] "package:methods" "Autoloads" "package:base" > Coincidentemente esse erro irritante que acontece com mais de um pacote vem acontecendo depois que eu comecei a usar o RStudio. Eu achei inicialmente porque eu tinha um monte de pacotes carregados, depois eu achei que era porque eu importei os dados com spss.get e não com o read.spss. Mas nada disse fez diferença nas tentativas e erros. Alguem tem uma luz pra eu tentar evitar esse erro? , Abraço forte e que a força esteja com você, 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

2 10

Re: [R-br] Gráfico WMF
by Ivan Bezerra Allaman 26 Set '11

26 Set '11

Leandro, Se você pediu uma sequência de 2 à 40 de 3 em 3, o último número que irá aparecer é claro que será o 38, por isso o número 40 não aparece. Caso queira que ele apareça, das diversas opções, você pode tentar, seq(2,41,3) ou seq(2,40,2). (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}

4 5

Gráfico WMF
by Leandro Marino 26 Set '11

26 Set '11

Caros, Preciso de uma ajuda. Não sei como resolver isto. Nas configurações abaixo, necessito que o número 40 seja exibido no eixo X. Alguma sugestão? plot(1:40,runif(40),axes=F,xlab='Eixo X',ylab='Eixo Y') axis(1,at=c(0,seq(2,40,3),40)) axis(2,seq(0,1,.1)) box() savePlot('c:/temp/temp.wmf') dev.off() Aproveitando, existe outra função para salvar em WMF (obrigatoriamente) mas que não seja necessário abrir uma nova janela? Isto faz com que seja consumido muito mais memória e processamento! Atenciosamente, Leandro Lins Marino Centro de Avaliação Fundação CESGRANRIO Rua Santa Alexandrina, 1011 - 2º andar Rio de Janeiro, RJ - CEP: 20261-903 ( (21) 2103-9600 R.:236 * <mailto:leandro@cesgranrio.org.br> leandro(a)cesgranrio.org.br P Antes de imprimir pense em sua responsabilidade e compromisso com o MEIO AMBIENTE Esta mensagem, incluindo seus anexos, pode conter informacoes privilegiadas e/ou de carater confidencial, nao podendo ser retransmitida sem autorizacao do remetente. Se voce nao e o destinatario ou pessoa autorizada a recebe-la, informamos que o seu uso, divulgacao, copia ou arquivamento sao proibidos. Portanto, se você recebeu esta mensagem por engano, por favor, nos informe respondendo imediatamente a este e-mail e em seguida apague-a.

2 1

Problemas com vetores
by Thiago Touya 26 Set '11

26 Set '11

Senhores, em meu trabalho, criei 3 vetores, sendo dois deles formados por 0's e 1's. De posse dos três criei uma matriz. Agora estou tentando adicionar mais uma coluna à matriz. Coluna esta que servira para identificar a situação dos dados do outro vetor que foi criado. No fim das contar preciso ter uma matriz com 4 colunas sendo que as duas primeiras formadas por 0's e 1's determinam a condição da terceira coluna e essa condição será encontrada na quarta coluna. Só que a quarta coluna é formada por valores to tipo "Censurado", "Curado" e "Falha". Tentei fazer da seguinte forma: # 1º vetor VetCens <- rbinom(500,1,0.6) # 2º vetor VetCura <- rbinom(500,1,0.3) # 3º vetor VetBN <- rnbinom(500,48,mu=13.5) # Criando a matriz MatSim <- cbind(VetCens, VetCura, VetBN) # Tentativa de criação do quarto vetor com as situações VetSituac <- function(param){ VetCens <- param[1] VetCura <- param[2] situacao <- rep(0,500) for(i in 1:500){ if (VetCens[i]==1 and VetCura[i]==0) situacao[i]==3 } if (VetCens==1 and VetCura==1){ situacao="Curado" } if (VetCens==0 and VetCura==1) or (VetCens==0 and VetCura==0) situacao="Falha" } } Estou tendo problemas com essa ultima função. Pra um individuo ser considerado CURADO devemos ter para os dois primeiros vetores os valores 1,1 Pra um indivíduo ser considerado CENSURADO devemos ter para os dois primeiros vetores os valores 1,0 Pra um indivíduo ser considerado FALHA devemos ter para os dois primeiros vetores os valores 0,1 ou 0,0. Alguém tem uma sugestão? Obrigado -- Att, Thiago Morais de Carvalho (61 - 8169-5094)

6 6

Off toppic - R Windows
by Leonard de Assis 26 Set '11

26 Set '11

Obrigado Norton 360 por detectar o r.dll como Virus, deletar da máquina sem me avisar e fazer eu perder um fim de semana tentando descobrir o que tinha de errado no meu R Só pra desabafar mesmo :( -- []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com

2 1

Compilando e utilizando códigos em C no R
by Gustavo Marcatti 26 Set '11

26 Set '11

Pessoal esse email é para aqueles que tem dificuldades em compilar um código em C e utilizá-lo no R (principalmente usuários Windows). Tem um pacote (inline) que facilita bastante esse processo e pode ser interessante até para usuários unix. Esse pacote compila o código gerando uma .dll (windows) temporária que depois é carregada no R. O código em C pode ser escrito no próprio script do R (normalmente eu deixo em um arquivo separado .r e quando quero utilizar a função utilizo o source(".r"), então o código é compilado/carregado e a função já fica pronta para uso. Depois de encerrar o R o dll é excluída do sistema. Segue um exemplo. # Exemplo: # Algoritmo que retorna o valor mínimo de um vetor # 1º: utilizando o loop for do R min2 <- function(vetor){ n <- length(vetor) Min <- vetor[1] for(i in 1:(n-1)){ if(Min > vetor[i+1]) Min <- vetor[i+1] } Min } # 2º: utilizando um codigo escrito em C #install.packages("inline", dep = T) require(inline) #--- Compilando a função que encotra o menor valor em um vetor ------------ #--- Só é necessário compilar uma vez -------------------------------------- # Argumentos da função em C sigMin <- signature( vetor="double", n="integer", Min = "double") # Código propriamente dito codeMin <- " *Min = vetor[0]; for (int i=0; i < *n - 1; i++) { if(*Min > vetor[i+1]) *Min = vetor[i+1]; }" # Para compilar a função funcMin <- cfunction(sigMin, codeMin, convention=".C") #--- Fim # apenas para facilitar o uso da função min3 <- function(vetor){ n <- length(vetor) Min <- 0 funcMin(vetor,n,Min)$Min } # Teste: set.seed(1) vec <- runif(1000000) system.time(minimo1 <- min2(vec)) minimo1 system.time(minimo2 <- min3(vec)) minimo2 # Mais exemplos podem ser encotrados no pacote inline. Gustavo Marcatti Eng. Florestal- UFV

1 0