[R] a repetition of simulation

Julian Burgos jmburgos at u.washington.edu
Thu Nov 15 21:44:39 CET 2007


summary(log_v)

Julian

sigalit mangut-leiba wrote:
> Hello,
> In addition to my question a few days ago,
> Now I have a matrix of the coefficients,
> how can I see all the P.Values (Pr(>|z|)) of the covariates from the 1000
> iterations?
> I tried names(log_v) and couldn'n find it.
> Thank you,
> Sigalit.
> 
> 
> On 11/13/07, Julian Burgos <jmburgos at u.washington.edu> wrote:
>> Well, the obvious (but perhaps not the most elegant) solution is put
>> everything in a loop and run it 600 times.
>>
>> coefficients=matrix(NA,ncol=3,nrow=600)
>>
>> for (loop in 1:600){
>>
>> [all your code here]
>>
>> coefficients[loop,]=coef(log_v)
>>
>> }
>>
>> That will give you a matrix with the coefficients of each model run in
>> each row.
>>
>> Julian
>>
>>
>> sigalit mangut-leiba wrote:
>>> I want to repeat the simulation 600 times and to get a vector of 600
>>> coefficients for every covariate: aps and tiss.
>>> Sigalit.
>>>
>>>
>>> On 11/13/07, Julian Burgos <jmburgos at u.washington.edu> wrote:
>>>> And what is your question?
>>>>
>>>> Julian
>>>>
>>>> sigalit mangut-leiba wrote:
>>>>> Hello,
>>>>> I have a simple (?) simulation problem.
>>>>> I'm doing a simulation with logistic model and I want to reapet it 600
>>>>> times.
>>>>> The simulation looks like this:
>>>>>
>>>>> z <- 0
>>>>> x <- 0
>>>>> y <- 0
>>>>> aps <- 0
>>>>> tiss <- 0
>>>>> for (i in 1:500){
>>>>> z[i] <- rbinom(1, 1, .6)
>>>>> x[i] <- rbinom(1, 1, .95)
>>>>> y[i] <- z[i]*x[i]
>>>>> if (y[i]==1) aps[i] <- rnorm(1,mean=13.4, sd=7.09) else aps[i] <-
>>>>> rnorm(1,mean=12.67, sd=6.82)
>>>>> if (y[i]==1) tiss[i] <- rnorm(1,mean=20.731,sd=9.751) else  tiss[i] <-
>>>>> rnorm(1,mean=18.531,sd=9.499)
>>>>> }
>>>>> v <- data.frame(y, aps, tiss)
>>>>> log_v <- glm(y~., family=binomial, data=v)
>>>>> summary(log_v)
>>>>>
>>>>> I want to do a repetition of this 600 times (I want to have 600
>> logistic
>>>>> models), and see all the coefficients of the covariates aps & tiss.
>>>>> Thanks in advance,
>>>>> Sigalit.
>>>>>
>>>>>       [[alternative HTML version deleted]]
>>>>>
>>>>> ______________________________________________
>>>>> R-help at r-project.org mailing list
>>>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>>>> PLEASE do read the posting guide
>>>> http://www.R-project.org/posting-guide.html
>>>>> and provide commented, minimal, self-contained, reproducible code.
>>>       [[alternative HTML version deleted]]
>>>
>>> ______________________________________________
>>> R-help at r-project.org mailing list
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



More information about the R-help mailing list