- R-br - DINF UFPR

Ajuda em modelo não linear
by Leonardo Guimarães de Oliveira 05 Jun '12

05 Jun '12

Como faço para encontrar o ponto na curva obtida na CMR abaixo onde a curva muda de direção. x=structure(list(Vt = c(0, 0.34, 3.68, 10.38, 21.94, 32.08, 48.54, 79.69, 94.82, 107.99), t = c(0L, 3L, 6L, 9L, 12L, 16L, 24L, 36L, 48L, 72L)), .Names = c("Vt", "t"), class = "data.frame", row.names = c(NA, -10L)) m4=nls(Vt~(107.99)/(1+exp(2-4*k*(t-g))), data=x, start=c(k=0.05,g=0.1)) m4 f4=function(x){(107.99)/(1+exp(2-4*0.02853*(x-8.71643)))} curve(f4, 1,100) Leonardo G. Oliveira Médico Veterinário Doutorando em Ciência Animal - EVZ/UFG

2 1

pergunta simples - diagonal de matrix
by Vitor Aguiar 04 Jun '12

04 Jun '12

Boa noite, se eu tenho uma matrix, por exemplo, 7x7 eu posso facilmente extrair a diagonal usando diag(mat). Mas e se eu quiser algo como "uma segunda diagonal", uma diagonal composta de mat[1,2] mat[2, 3] mat[3, 4] mat[4, 5] mat[5, 6] mat[6, 7]? Tem alguma forma simples de fazer? obrigado, Vitor

2 1

Re: [R-br] Alteração em função de pacote (Junior Beleti)
by Junior Beleti 04 Jun '12

04 Jun '12

Olá Benilton, meus estudos se iniciaram após uma conversa com a professora Ana Beatriz, visto que ela tinha problemas de desempenho no trabalho desenvolvido por ela e pelo professor Paulo. Eu fui analisando os pontos onde ela me disse que teria problemas, e identifiquei na volta.cokri um maior tempo de processamento. Vou processar novamente os scripts, mas em minha máquina o tempo de processamento fui bem maior. Volto a postar os resultados. E o professor Paulo está certo pois esse trabalho não visava desempenho, e como ele mesmo citou, vi uma oportunidade para utilizar processamento paralelo e tentar ajudar no desempenho dessa aplicação. Muito obrigado pela ajuda. Junior Beleti.

1 0

borda externa de mapa
by Renato Braga 04 Jun '12

04 Jun '12

Lista, como eu faco para plotar somente a borda externa de um mapa no formato shapefile? Grato, Renato

2 2

Python para R
by Daniel Marcelino 03 Jun '12

03 Jun '12

Caros, alguem da lista conhece programação python? Estou há dias tentando desenvolver uma rotina em R para explorar dados eleitorais, mas não tenho conseguido ir muito longe. Eu tenho um modelo do que mais ou menos eu preciso, mas a rotina está em python. Alguem encara me ajudar nisso? Daniel ----------------------------------------------------- My contacts: Phone 1+514 343 6111 Ext 3799 Skype: d.marcelino 3200 Jean Brillant, Office C5075 Montreal, QC; H3T 1N8 Canada -----------------------------------------------------

3 2

Combinar vetores com tamanhos diferentes...
by andrebvs 03 Jun '12

03 Jun '12

7 13

borda externa
by André Figueiras Rutz 02 Jun '12

02 Jun '12

Caros, como é que faz para inserir num gráfico a moldura externa (só a linha preta em volta do gráfico?) [image: Inline image 3] grato!

3 8

Gráfico
by Paulo Ricardo Silva Moreira 01 Jun '12

01 Jun '12

Boa tarde prezados, Estou precisando fazer um gráfico, plotando duas variáveis com suas respectivas escalas, uma de cada lado do gráfico. Alguém sabe me dizer qual comando eu poderia usar!? Desde já agradeço, Paulo Ricardo.

4 5

Alteração em função de pacote (novo tópico)
by Junior Beleti 01 Jun '12

01 Jun '12

Olá Benilton, tentando explicar melhor meu problema: Faço uso de um pacote chamada geoComp. Tal pacote é utilizado de forma geral para modelar o padrão espacial de dados composicionais. Foi um trabalho da professora Ana Beatriz. O problema está quanto ao tempo de processamento para realizar o trabalho necessário. Analisando mais a fundo, verifiquei que o maior problema quanto a desempenho está na função "volta.cokri" desse pacote. Nessa função existe uma chamada à função "mvrnorm" e é justamente nessa chamada que ocorre o maior problema de desempenho. Verifiquei também que na função "mvrnorm" existe uma chamada para a tal função "eigen", que se trata de um código-objeto fortran. Inicialmente utilizei a função "mclapply" para tentar melhorar tal desempenho, mas não obtive muito sucesso. O que estava pensando quando insisti na ideia de criar a nova função no pacote base, foi o de que se chamando funções "diferentes", ou seja, criando funções idênticas cada uma com o espaço de endereçamento e nomes diferentes, poderia aumentar o desempenho. A seguir, está um script utilizado: require(geoComp) require(MASS) require(statmod) require(geoR) data(pivo) dados <- pivo[,c(6,7,8,1,2)] dados <- as.geoComp(dados) bor <- cbind(c(0,seq(0,200,l=100),0),c(0,sqrt(200^2-seq(0,200,l=100)^2),0)) ## bor é uma matriz 2x102 estima <- mec(dados) gr <- pred_grid(bor, by=4) md.cov.ck <- cokrigagem(estima[[1]]$Estimativas, loc=gr,dados.comp=dados) preditos.gh <- volta.quad(md.cov.ck,n.pontos=7,Variancia=FALSE) preditos.gh <- data.frame(preditos.gh) write.table(preditos.gh,"pred_by4k7ns1000MBM.txt") preditos.simu <- volta.cokri(md.cov.ck,num.simu=1000,int.conf=0.95) preditos.simu <- data.frame(preditos.simu[[1]]) preditos.simu.ic <- data.frame(preditos.simu[[2]]) write.table(preditos.simu,"predsimu_by4k7ns1000MBM.txt") write.table(preditos.simu.ic,"predsimuic_by4k7ns1000MBM.txt") O pacote geoComp pode ser adquirido no endereço: ( http://www.din.uem.br/~pg45178/geocomp/geoComp.tar.gz) Um resumo do trabalho pode ser encontrado em: ( http://www.leg.ufpr.br/lib/exe/fetch.php/pessoais:abtmartins:rbb09_2.pdf) Obrigado pela atenção. Carlos.

3 3

Re: [R-br] Alteração em função de pacote
by Junior Beleti 01 Jun '12

01 Jun '12

Bom dia Fernando Mayer, segui suas orientações, mas obtive o seguinte erro ao tentar instalar o pacote (R CMD INSTALL base_2.15.0.tar.gz). * installing to library ‘/usr/local/lib/R/library’ * installing *source* package ‘base’ ... ** R ** demo ** inst ** preparing package for lazy loading Erro em asNamespace(ns, base.OK = FALSE) : operation not allowed on base namespace * removing ‘/usr/local/lib/R/library/base’ * restoring previous ‘/usr/local/lib/R/library/base’ Obrigado pela atenção e se tiver alguma outra sugestão. Att, -- Carlos

4 3