Displaying 20 results from an estimated 7000 matches similar to: "Statistical question about logistic regression simulation"
2004 May 01
2
Generating Lognormal Random variables (PR#6843)
Full_Name: Anthony Gichangi
Version: 1.90
OS: Windows XP Pro
Submission from: (NULL) (130.225.131.206)
The function rlnorm generates negative values for lognormal distribution.
x- rlnorm(1000, meanlog = 0.6931472, sdlog = 1)
Regards
Anthony
2008 May 04
1
Is my understanding of rlnorm correct?
rlnorm takes two 'shaping' parameters: meanlog and sdlog.
meanlog would appear from the documentation to be the log of the mean.
eg if the desired mean is 1 then meanlog=0.
So to generate random values that fit a lognormal distribution I would
do this:
rlnorm(N , meanlog = log(mean) , sdlog = log(sd))
But when I check the mean I don't get it when sdlog>0. Interestingly I
2009 Feb 11
3
Generating Numbers With Certain Distribution in R
Dear all,
Is there a way to generate K numbers of integer (K = 10^6).
The maximum value of the integer is 200,000 and minimum is 1.
And the occurrences of this integer follows
a lognormal distribution.
- Gundala Viswanath
Jakarta - Indonesia
2005 Jun 29
2
MLE with optim
Hello,
I tried to fit a lognormal distribution by using optim. But sadly the output
seems to be incorrect.
Who can tell me where the "bug" is?
test = rlnorm(100,5,3)
logL = function(parm, x,...) -sum(log(dlnorm(x,parm,...)))
start = list(meanlog=5, sdlog=3)
optim(start,logL,x=test)$par
Carsten.
[[alternative HTML version deleted]]
2010 Jun 21
2
How to predict the mean and variance of the dependent variable after regression
Hi, folks,
As seen in the following codes:
x1=rlnorm(10)
x2=rlnorm(10,mean=2)
y=rlnorm(10,mean=10)### Fake dataset
linmod=lm(log(y)~log(x1)+log(x2))
After the regression, I would like to know the mean of y. Since log(y) is
normal and y is lognormal, I need to know the mean and variance of log(y)
first. I tried mean (y) and mean(linmod), but either one is what I want.
Any tips?
Thanks in
2009 May 29
1
Mean of lognormal in base-2
Hi, Does anyone know what the mean value of a lognormal distribution in base-2 is? I am simulating stochastic population growth and if I were working in base-e, I would do:lambda <- 1.1 #multiplicative growth rates <- 0.6 #stochasticity (std. dev)lognormal <- rlnorm(100000, log(lambda) - (s^2)/2, s)## or lognormal <- exp( rnorm( 100000, log(lambda) - (s^2)/2,
2002 Dec 10
1
Lognormal distribution
I am trying to fit a lognormal distribution to a set of data and test its
goodness of fit with regard to predicted values.
I managed to get so far:
> y <- c(2,6,2,3,6,7,6,10,11,6,12,9,15,11,15,8,9,12,6,5)
> library(MASS)
> fitdistr(y,"lognormal",start=list(meanlog=0.1,sdlog=0.1))
meanlog sdlog
1.94810515 0.57091032
(0.12765945) (0.09034437)
But I would
2012 Oct 11
2
Help on probability distribution question
Dear All,
I have a questions I would like to ask about and wonder if you have any thoughts to make it work in R.
1. I work in the field of medicine where physiologic variables are often simulated, and they can not have negative values. Most often the assumption is made to simulate this parameters with a normal distribution but in the "log-domain" to avoid from negative values to be
2010 Feb 24
1
Remove missing observations
Hi everyone
I have the following problem: My dataframe has 3 variables: ID, Year and
and an outcome variable. The dataframe contains repeated measurements
because the subjects filled out a questionnaire every year. The time span
covers 2 years.
Now I want to check if there is a significant change in the outcome over
the 2 years with a paired wilcox.test. The problem: Not every subject
2005 Sep 27
1
Producing empirical bayes estimates in disease mapping for lognormal model
I'm trying to produce empirical bayes estimates based on the lognormal
model in disease mapping
Is there a way this can be done in R?
thanks
Oarabile
2007 Oct 03
1
offset in survreg
Hello,
I have a question regarding the use of an offset term with survreg(),
in the Survival library. In particular, I am trying to figure out on
what scale the offset term should be.
Here's a simple example with no censoring and no coefficients:
---------
y = rlnorm(1000, meanlog = 10, sdlog = 2)
delta = rep(1, 1000)
int = rep(1, 1000)
survreg(Surv(y,delta)~offset(10*int), dist =
2006 Nov 18
3
Random sample from log-normal distribution
Dear all R users,
Please forgive me if my question is too trivial.
Suppose I have two variables, (x,y) which is
log-normally distributed with expected value (mu1,
mu2) and some variance-covariance matrix. Now I want
to draw a random sample of size 1000 from this
distribution. Is there any function available to do
this?
Thanks and regards,
Megh
2009 May 20
1
sem with categorical data
I am trying to run a confirmatory factor analysis using the SEM package. My
data are ordinal. I have read
http://socserv.mcmaster.ca/jfox/Misc/sem/SEM-paper.pdf.
When I apply the hetcor function, I receive the following error:
Error in checkmvArgs(lower = lower, upper = upper, mean = mean, corr = corr,
:
at least one element of 'lower' is larger than 'upper'
Example:
2003 Apr 09
3
plotting the lognormal density curve
I am trying to plot a lognormal density curve on top of an existing
histogram. Can anybody suggest a simple way to do this? Even if someone
could just explain how to plot a regular normal density curve on top of an
existing histogram, it would be a big help.
Also, is there some way to search through the R-help archives other than
simple browsing?
Thank you so much. Your help and time is greatly
2005 Jan 07
3
lognorm
Hi!
I 've a problem to have a lognorm distribution with
mean=1 and var (or sigma)=1.
rlnorm(1000,0,0)
rlnorm(1000,1,1)
rlnorm(1000,0,1)
.... ?
Can you help me?
2009 May 31
1
Bug in truncgof package?
Dear R-helpers,
I was testing the truncgof CRAN package, found something that looked
like a bug, and did my job: contacted the maintainer. But he did not
reply, so I am resending my query here.
I installed package truncgof and run the example for function ad.test. I
got the following output:
set.seed(123)
treshold <- 10
xc <- rlnorm(100, 2, 2) # complete sample
xt <- xc[xc >=
2011 Jul 10
1
Chebyshev Inequality — MVUE
Hello,
I was interested in trying to write an R script to calculate a UCL for a lognormal distribution using the Chebyshev Inequality — MVUE Approach (based on EPA’s guidance found in http://www.epa.gov/oswer/riskassessment/pdf/ucl.pdf). This looks like it should be straight forward, but I am need to calculate an MVUE for the population mean and an MVUE for the population variance, which requires
2003 Jul 25
5
named list 'start' in fitdistr
Hi R lovers!
I'd like to know how to use the parameter 'start' in the function
fitdistr()
obviously I have to provide the initial value of the parameter to optimize
except in the case of a certain set of given distribution
Indeed according to the help file for fitdistr
" For the following named distributions, reasonable starting values
will be computed if `start'
2011 Jul 06
2
how to best present concentrated data points/ ggplot2
Hi all,
I am trying to plot a weighted density plot for two different types and want to show the data points on the x axis.
The code is as follows. The data points are very concentrated. Is there a better way to present it( should I set the alpha value or something else)?
Thanks!
YL
library(ggplot2)
x <- rnorm(10000)
a <- rnorm(5000)
b <- rnorm(5000)
weights.x <- abs(a/sum(a))
2010 Dec 09
1
Calculating odds ratios from logistic GAM model
Dear R-helpers
I have a question related to logistic GAM models. Consider the following
example:
# Load package
library(mgcv)
# Simulation of dataset
n <- 1000
set.seed(0)
age <- rnorm(n, 50, 10)
blood.pressure <- rnorm(n, 120, 15)
cholesterol <- rnorm(n, 200, 25)
sex <- factor(sample(c('female','male'), n,TRUE))
L <-