- R-br - DINF UFPR

Problema numérico no gnm
by Davi 08 Mai '13

08 Mai '13

Boa noite a todos! Estou tentando simular de um modelo não-linear generalizado e depois ajustar esse modelo aos dados gerados na simulação. A variável resposta pode seguir Poisson ou Binomial, não tem problema. A média é modelada por: n*(1-exp(-(b0+b1*x))*(exp(-gx) ), em que x é uma variável quantitativa, n pode ser uma variável offset, no caso da Poisson, ou o número de ensaios Bernoulli, no caso da binomial. Aí vai o exemplo que vem dando erro, que eu suponho ter a ver com o optim, devido aos valores muito pequenos dos parâmetros b0, b1 e g e do valor da offset n ser muito grande. Agradeço a todos pela atenção! Abraços. # no R require(gnm) x <- sort( rep(0:5,3) ) b0 <- 1.5*1e-6 b1 <- 3*1e-6 g <- 3*1e-1 n <- rep(1e+8,length(x)) p <- (1 - exp(-(b0+b1*x)) ) * exp(-g*x) ; p mu <- n*p ## Funções não lineares # Parametrização original margo.betas <- function(x){ list(predictors=list(beta0=1,beta1=1,gamma=1), variables=list(substitute(x)), term=function(predictors,variables){ pred<-paste("(1-exp(",variables[1],"+",variables[2],"*",predictors[1],"))*(e xp(", variables[3],"*",predictors[1],"))",sep="") }) } class(margo.betas) <- "nonlin" # Reparametrização, com b0 = 1/(ln k0), b1=1/(ln k1) e g=1/(ln k2) margo.ks <- function(x){ list(predictors=list(k0=1,k1=1,k2=1), variables=list(substitute(x)), term=function(predictors,variables){ pred<-paste("(1-",variables[1],"*(",variables[2],"^",predictors[1],"))*(", variables[3],"^",predictors[1],")",sep="") }) } class(margo.ks) <- "nonlin" ######## ABORDAGEM POISSON set.seed(0119) y <-rpois(length(x),lambda=mu) eta.margo1 <- formula(y ~ -1 + Mult(offset(n), margo.ks(x) ) ) eta.margo2 <- formula(y ~ -1 + Mult(offset(n), margo.betas(x) ) ) gnm(eta.margo1,family=poisson,tolerance=1e-50) gnm(eta.margo2,family=poisson,tolerance=1e-50) ######## ABORDAGEM BINOMIAL set.seed(0119) y <-rbinom(length(x),size=n,prob=p) eta.margo1 <- formula(cbind(y,n-y) ~ -1+ margo.ks(x) ) eta.margo2 <- formula(cbind(y,n-y) ~ -1+ margo.betas(x) ) gnm(eta.margo1,family=binomial) gnm(eta.margo2,family=binomial) ___________________________________________________ Davi Butturi-Gomes Doutorando em Estatística e Experimentação Agronômica - ESALQ /USP - Piracicaba Mestre em Biometria - IB/UNESP - Botucatu Ecólogo - IB/UNESP - Rio Claro Skype: davibg87 (19)8206.6565

1 0

Krigagem
by Arthur Telles Calegario 08 Mai '13

08 Mai '13

Bom dia Senhores,Estou com um pequeno problemas em relação a Krigagem e talvez os senhores tenham uma resposta.Na verdade o problema esta em preparar os dados. Estou trabalhando com evapotranspiração no Rio Doce e possuo as coordenadas das estações meteorológicas em graus. Contudo, para a construção do semivariograma, necessito que a distância entre as mesmas seja linear.Assim devo projeta-las. É um procedimento simples quando todos as estações estão em um mesmo 'datum'. Contudo, as estações que estão do 'datum' 23S são insuficientes para a construção do semivariograma. Assim necessito saber como faço para linearizar as distâncias entre estações meteorológicas cujos limites extrapolam um 'datum' apenas. Entenderam? rs. Agradeço desde já,Arthur Calegario.

11 33

Krigagem
by Tito Conte 03 Mai '13

03 Mai '13

Bom dia, Alguém poderia me ajudar com krigagem? Estou usando o fields, e tenho uma matriz com 57703, com x, y, z pontos, quando uso o comando Krig ele da vários erros. O comando é K = Krig(xy,z) xy são as posições e z é a altura Tito Conte

4 3

pp-plot
by Julia Maria 02 Mai '13

02 Mai '13

Pessoal, o R tem alguma forma pré definida para se fazer um Probability-Probability plot? Por exemplo, o pacote circular tem a função: pp.plot(x), que verifica o ajuste a uma distribuição von-mises. E se eu quiser fazer o mesmo com outras ditribuições? Grata

1 0

set.seed
by Daniel Marcelino 02 Mai '13

02 Mai '13

Caros, estou realizando um experimento e acabei constatando que a semente que eu distribuir é fundamental para eu obter as amostras binárias. Eu estou interessado em distribuir 25 unidades para tratamento e o restante para controle. Há alguma maneira de saber, sem ser por tentativa e erro, que minha distribuição terá 25 unidades em cada categoria? Por exemplo, > # 50 Bernoulli trials > set.seed(51) > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50)) > table(T) T 0 1 22 28 > set.seed(12345) > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50)) > table(T) T 0 1 23 27 > set.seed(123) > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50)) > table(T) T 0 1 25 25 > set.seed(2) > T <- sample(c(0,1), 50, replace=TRUE, prob=c(25/50, 25/50)) > table(T) T 0 1 25 25 -- \begin{signature} Daniel Marcelino ☁ dm.silva(a)umontreal.ca ☎ (514) 343 6111 #3799 Skype: d.marcelino ✎ 3200 Jean Brillant, Office C5071 Montreal, QC; H3T 1N8 Canada \end{signature} "Small steps toward a much better world"

4 5

excluir linhas que tenham códigos repetidos
by Fátima Lima Paula 01 Mai '13

01 Mai '13

Pessoal, existe algum comando que exclua as linhas que apresentam códigos repetidos? Não seria deixar apenas um não. É excluir todas. Exemplo nome=c("Ricardo","Maria","Tereza","Pedro","José", "Germano","Paulo","Ana") sobrenome=c("Cunha","Andrade","Silva","Soares","Silva","Lima","Pitanga","Amora") aih=c(6867,6798,8903,7865,8903,6745,6867,8903) ano=c(1944,1945,1923,1954,1947,1945,1946,1944) dataint=c(20121004,20120907,20120703,20120301,20120805,20121209,20120907,20120703) sexo=c("M","F","F","M","M","M","M","F") banco1=cbind(nome,sobrenome,aih,sexo,ano,dataint) banco1=data.frame(banco1) banco1 nome sobrenome aih sexo ano dataint 1 Ricardo Cunha 6867 M 1944 20121004 2 Maria Andrade 6798 F 1945 20120907 3 Tereza Silva 8903 F 1923 20120703 4 Pedro Soares 7865 M 1954 20120301 5 José Silva 8903 M 1947 20120805 6 Germano Lima 6745 M 1945 20121209 7 Paulo Pitanga 6867 M 1946 20120907 8 Ana Amora 8903 F 1944 20120703 Quero excluir as linhas que apresentam as aihs repetidas (6867 e 8903), de forma que o banco fique assim: nome sobrenome aih sexo ano dataint 1 Maria Andrade 6798 F 1945 20120907 2 Pedro Soares 7865 M 1954 20120301 3 Germano Lima 6745 M 1945 20121209 Obrigada -- "Minha felicidade depende da qualidade dos meus pensamentos"

5 5

Função maxDrawdown
by Marcos Vagner Macedo 01 Mai '13

01 Mai '13

Boa noite pessoal, Não estou conseguindo utilizar a função maxDrawdown do pacote PerformanceAnalytics. Alguém pode me ajudar? Estou utilizando dias úteis no formato zoo. Quando utilizo a função tenho o seguinte resultado: > t(round(maxDrawdown(fi_97_12[,2:4]),4)) Worst Drawdown Ibovespa NaN ALFA.DINAMICO.FIA NaN ALFA.FIA NaN Warning messages: 1: In merge.zoo(fx, .xts(, .index(x))) : Index vectors are of different classes: integer POSIXct 2: In merge.zoo(fx, .xts(, .index(x))) : Index vectors are of different classes: integer POSIXct 3: In merge.zoo(fx, .xts(, .index(x))) : Index vectors are of different classes: integer POSIXct > Muito Obrigado, Marcos Macedo

1 0

Juntar colunas
by Edson Lira 01 Mai '13

01 Mai '13

Caros amigos, tenho um banco de dados aonde as quatro primeiras colunas, por problema de importação é o nome do indivíduo do banco, qual a rotina que posso juntar as colunas com os respectivos espaços. Exemplo: C1 C2 C3 C4 ADRIANA DE SOUZA MALHEIROS [ ]'S. Edson Lira Estatístico Manaus-Amazonas

4 5

package Babel com knitr [Resolvido]
by Fernando Colugnati 30 Abr '13

30 Abr '13

Oi Benilton, resolveu com o latin1 mesmo. Pouco antes de sua resposta tinha conseguido um exemplo! Abs e obrigado -- Fernando A.B. Colugnati

1 0

package Babel com knitr
by Fernando Colugnati 30 Abr '13

30 Abr '13

Prezados, estou com um problema na questão de internacionalização com o Latex, via R, usando o knitr (não o sweave). Apesar de incluir \usepackage[portuguese]{babel} não esto tendo sucesso na compilação para PDF com o pdfLaTex, os caracteres saem errados. Em tempo, estou rodando pelo Rstudio, em um MacPro (pensei em ser algum problema de codificação...). A instalação que tenho do Tex é a padrão para Mac. Alguém já passou por isso? Segue um pedaço de código, mas com todo o preambulo que estou usando: \documentclass{article} \usepackage[portuguese]{babel} \usepackage[sc]{mathpazo} \usepackage[T1]{fontenc} \usepackage{geometry} \geometry{verbose,tmargin=1.5cm,bmargin=1.5cm,lmargin=2.5cm,rmargin=2.5cm} \setcounter{secnumdepth}{2} \setcounter{tocdepth}{2} \usepackage{url} \usepackage[unicode=true,pdfusetitle, bookmarks=true,bookmarksnumbered=true,bookmarksopen=true,bookmarksopenlevel=2, breaklinks=false,pdfborder={0 0 1},backref=false,colorlinks=false] {hyperref} \hypersetup{ pdfstartview={XYZ null null 1}} \usepackage{breakurl} \begin{document} \section*{Matemática $5^o$ ano} Isso \'{e} um teste para acentuação! acentua\'{c}\˜{a}0, assim est\'{a} correto! \end{document} Obrigado -- Fernando A.B. Colugnati

2 1