[R] Ifelse leading to inconsistent result

Neville O'Reilly noreilly at stat.rutgers.edu
Wed Jun 26 19:18:04 CEST 2013


You are right, I inadvertently deleted part of the script - so I ended up with non-mutually exclusive regions. My apologies


Neville O'Reilly, Ph.D 
Associate Director, Financial Statistics & Risk Management, 
Rm 479 Hill Center, Rutgers University, 110 Freylinghuysen Rd., 
Piscataway, NJ 08854. 
848-445-7669 
http://fsrm.rutgers.edu/ 

----- Original Message -----
From: "David Carlson" <dcarlson at tamu.edu>
To: "Neville O'Reilly" <noreilly at stat.rutgers.edu>, r-help at r-project.org
Sent: Wednesday, June 26, 2013 1:15:06 PM
Subject: RE: [R] Ifelse leading to inconsistent result

I don't see that you have set up mutually exclusive ranges. If we
modify your code so save the maxlogP and minlogP values:

niter = 1e5 # number of iterations is 10^5
CountLoss = rep(0,niter)
CountProf = rep (0,niter)
maxlogP <- rep(0, ninter)
minlogP <- rep(0, ninter)
set.seed(2009) # enables reproducibility of result if script run
again"
for (i in 1:niter)
{
  r = rnorm(100,mean=.05/253,
            sd=.23/sqrt(253)) # generate 100 random normal numbers
  logPrice = log(1e6) + cumsum(r) #vector of 100 days log prices
  maxlogP[i] = max(logPrice) # max price over next 100 days
  minlogP[i] = min(logPrice)
  CountLoss[i] <- ifelse (minlogP[i] < log(950000), 1, ifelse
(maxlogP[i] > log (1000000), 0, 1))
  CountProf[i] <- ifelse (maxlogP[i] < log (1100000),0,1)
}
both <- which(CountLoss+CountProf>1)
length(both)
# 18484
head(both)
# [1]  1  5  9 12 15 25
ifelse (minlogP[1] < log(950000), 1, ifelse (maxlogP[1] > log
(1000000), 0, 1))
# [1] 1
ifelse (maxlogP[1] < log(1100000),0,1)
# [1] 1
exp(maxlogP[1])
# [1] 1204589
exp(minlogP[1])
# [1] 932747.5

Your first simulation meets both criteria along 18483 others!

-------------------------------------
David L Carlson
Associate Professor of Anthropology
Texas A&M University
College Station, TX 77840-4352

-----Original Message-----
From: r-help-bounces at r-project.org
[mailto:r-help-bounces at r-project.org] On Behalf Of Neville O'Reilly
Sent: Wednesday, June 26, 2013 10:41 AM
To: r-help at r-project.org
Subject: [R] Ifelse leading to inconsistent result

I have used ifelse in count variables to count the number of times
in a simulation the values of a vector of logprice fall within
mutually exclusive ranges. However, there is a double count in the
result i.e. i am getting output indicating values falling in
mutually exclusive ranges. Here is the code and result
R script
niter = 1e5 # number of iterations is 10^5
CountLoss = rep(0,niter)
CountProf = rep (0,niter)
set.seed(2009) # enables reproducibility of result if script run
again"
for (i in 1:niter)
{
  r = rnorm(100,mean=.05/253,
            sd=.23/sqrt(253)) # generate 100 random normal numbers
  logPrice = log(1e6) + cumsum(r) #vector of 100 days log prices
  maxlogP = max(logPrice) # max price over next 100 days
  minlogP = min(logPrice)
  CountLoss[i] <- ifelse (minlogP < log(950000), 1, ifelse (maxlogP
> log (1000000), 0, 1))
  CountProf[i] <- ifelse (maxlogP < log (1100000),0,1)
}
sum(CountLoss)
mean(CountLoss) # fraction of times out of niter that stock is sold
for a loss in a 100 day period
sum(CountProf)
mean(CountProf) # fraction of times out of niter that stock is sold
for a profit in a 100 day period

Output
sum(CountLoss)
[1] 64246
> mean(CountLoss) # fraction of times out of niter that stock is
sold for a loss in a 100 day period
[1] 0.64246
> sum(CountProf)
[1] 51857
> mean(CountProf) # fraction of times out of niter that stock is
sold for a profit in a 100 day period
[1] 0.51857

CountLoss and CountProf should sum to less than the number of
interations. When I troubleshoot by reducing the number of
iterations and that size of the logprice, I can't reproduce the
contradicion.

______________________________________________
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