simulação de mutinomiais
Senhores, estou interessado em simular uma amostra de multinomiais em que suas probabilidades sao geradas por uma regressão logistica, isto é Y ~ Multinomial(1, Pk0i ,Pk1i ,Pk2i , Pk3i) em que i = 1, ..., 100 ( tamanho da amostra) e k= 1,2,3. Assim, k assume 3 valores para cada i. A dúvida é: como faço para informar essa variação das probabilidades na função rmultinom? ou teria outra maneira de fazer isso? Grato, Márcio
Então, de que pacote é essa função? Têm algum CMR para trabalharmos em cima? Em geral, as funções que geram variáveis aleatórias permitem você passar o valor dos parâmetros, por exemplo, com a rnorm() eu passo rnorm(10, mean=3) para simular de média 3, não tem o equivalente nessa sua função para especificar as probabilidades? Caso não tenha você pode simular de uma uniforme e fazer os cortes de forma que o tamanho do intervalo corresponda às probabilidades desejadas. Leia o guia de postagem. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Geralmente eu simulo multinomial com o sample. dê uma olhada com ?sample []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com Em 16/05/2012 16:49, Walmes Zeviani escreveu:
Então, de que pacote é essa função? Têm algum CMR para trabalharmos em cima? Em geral, as funções que geram variáveis aleatórias permitem você passar o valor dos parâmetros, por exemplo, com a rnorm() eu passo rnorm(10, mean=3) para simular de média 3, não tem o equivalente nessa sua função para especificar as probabilidades? Caso não tenha você pode simular de uma uniforme e fazer os cortes de forma que o tamanho do intervalo corresponda às probabilidades desejadas. Leia o guia de postagem.
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br <mailto:walmes@ufpr.br> twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes <http://www.leg.ufpr.br/%7Ewalmes> linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Sim na função rmultinom() tem como especificar as probabilidades, mas no meu caso não tenho um valor especifico desses probabilidades. As probabilidade são geradas por uma regressão logistica e a parte delas eu gero a multinomial, alem disso para cada i amostrado eu terei 4 vetores de probabilidades. Assim, não sei como faço informar essa variaçao do vetor de probabilidades na função. Desde já grato pela atenção Márcio ________________________________ De: Walmes Zeviani <walmeszeviani@gmail.com> Para: r-br@listas.c3sl.ufpr.br; Márcio Augusto Ferreira Rodrigues <marcioaugusto28@yahoo.com.br> Enviadas: Quarta-feira, 16 de Maio de 2012 16:49 Assunto: Re: [R-br] simulação de mutinomiais Então, de que pacote é essa função? Têm algum CMR para trabalharmos em cima? Em geral, as funções que geram variáveis aleatórias permitem você passar o valor dos parâmetros, por exemplo, com a rnorm() eu passo rnorm(10, mean=3) para simular de média 3, não tem o equivalente nessa sua função para especificar as probabilidades? Caso não tenha você pode simular de uma uniforme e fazer os cortes de forma que o tamanho do intervalo corresponda às probabilidades desejadas. Leia o guia de postagem. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Dúvida: Não seria possível você passar como argumento prob da função rmultinon uma função que gere o seu valor da logit? att, FH 2012/5/17 Márcio Augusto Ferreira Rodrigues <marcioaugusto28@yahoo.com.br>:
Sim na função rmultinom() tem como especificar as probabilidades, mas no meu caso não tenho um valor especifico desses probabilidades. As probabilidade são geradas por uma regressão logistica e a parte delas eu gero a multinomial, alem disso para cada i amostrado eu terei 4 vetores de probabilidades. Assim, não sei como faço informar essa variaçao do vetor de probabilidades na função.
Desde já grato pela atenção
Márcio
________________________________ De: Walmes Zeviani <walmeszeviani@gmail.com> Para: r-br@listas.c3sl.ufpr.br; Márcio Augusto Ferreira Rodrigues <marcioaugusto28@yahoo.com.br> Enviadas: Quarta-feira, 16 de Maio de 2012 16:49 Assunto: Re: [R-br] simulação de mutinomiais
Então, de que pacote é essa função? Têm algum CMR para trabalharmos em cima? Em geral, as funções que geram variáveis aleatórias permitem você passar o valor dos parâmetros, por exemplo, com a rnorm() eu passo rnorm(10, mean=3) para simular de média 3, não tem o equivalente nessa sua função para especificar as probabilidades? Caso não tenha você pode simular de uma uniforme e fazer os cortes de forma que o tamanho do intervalo corresponda às probabilidades desejadas. Leia o guia de postagem.
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Olá, gostaria de saber se há a possibilidade de eu não receber esses e-mail na minha caixa de entrada, contudo, sempre que for de meu interesse consultá-los em algum site. Atenciosamente, 2012/5/17 FHRB Toledo <fernandohtoledo@gmail.com>
Dúvida:
Não seria possível você passar como argumento prob da função rmultinon uma função que gere o seu valor da logit?
att, FH
2012/5/17 Márcio Augusto Ferreira Rodrigues <marcioaugusto28@yahoo.com.br
:
Sim na função rmultinom() tem como especificar as probabilidades, mas no meu caso não tenho um valor especifico desses probabilidades. As probabilidade são geradas por uma regressão logistica e a parte delas eu gero a multinomial, alem disso para cada i amostrado eu terei 4 vetores de probabilidades. Assim, não sei como faço informar essa variaçao do vetor de probabilidades na função.
Desde já grato pela atenção
Márcio
________________________________ De: Walmes Zeviani <walmeszeviani@gmail.com> Para: r-br@listas.c3sl.ufpr.br; Márcio Augusto Ferreira Rodrigues <marcioaugusto28@yahoo.com.br> Enviadas: Quarta-feira, 16 de Maio de 2012 16:49 Assunto: Re: [R-br] simulação de mutinomiais
Então, de que pacote é essa função? Têm algum CMR para trabalharmos em cima? Em geral, as funções que geram variáveis aleatórias permitem você passar o valor dos parâmetros, por exemplo, com a rnorm() eu passo rnorm(10, mean=3) para simular de média 3, não tem o equivalente nessa sua função para especificar as probabilidades? Caso não tenha você pode simular de uma uniforme e fazer os cortes de forma que o tamanho do intervalo corresponda às probabilidades desejadas. Leia o guia de postagem.
À disposição. Walmes.
==========================================================================
Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218
==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
código
mínimo reproduzível.
R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Renato Marciano Maciel* * *Departamento de Engenharia de Biossistemas - DEPEB* *Universidade Federal de São João del Rei - UFSJ *Praça Dom Helvécio, 74, Fábricas São João del-Rei, Minas Gerais / CEP: 36301-160 Telefone/fax: (32) 3379-2541 (35) 9808-1270
Sim, isso é possível. Edite as suas opções de membro na lista pedindo para receber não mais mensagens individuais e sim o digest (conteúdo de um dia). Acesse os links no rodapé da mensagem. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Então é só escrever a função logística, conforme se queira, e determinar essas probabilidades e simular. Um ponto de partida pode ser esse logis <- function(x, A, B) 1/(1+exp(-(A+B*x))) curve(logis(x, 1, 1), -9, 4) curve(logis(x, 2.3, 1), add=TRUE) curve(logis(x, 3.6, 1), add=TRUE) x <- -2 px <- logis(-2, c(1,2.3,3.6), 1) abline(v=x, h=px, lty=3) probs <- diff(c(0,px,1)) probs sum(probs) amostr <- sample(1:4, size=1000, prob=probs, replace=TRUE) prop.table(table(amostr)) À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (5)
-
FHRB Toledo -
Leonard de Assis -
Márcio Augusto Ferreira Rodrigues -
Renato Maciel -
Walmes Zeviani