Dear R helpers,
I have following input files. (Actually they are more than 10 rates but here i
am considering only 2 rates to write my problem)
rate1.csv
min1 max1 min2 max2 min3 max3
1.05 1.30 1.30 1.65 1.65 1.99
rate2.csv
min1 max1 min2 max2 min3 max3
2.05 2.30 2.30 2.65 2.65 2.99
The simple R code I had used to read these files was as given below -
# OLD CODE
rates1 = read.csv('rate1.csv)
rates2 = read.csv('rate2.csv')
rate1_min1 = rates1$min1
rate1_max1 = rates1$max1
rate1_min2 = rates1$min2
rate1_max2 = rates1$max2
rate1_min3 = rates1$min3
rate1_max3 = rates1$max3
rate2_min1 = rates2$min1
rate2_max1 = rates2$max1
rate2_min2 = rates2$min2
rate2_max2 = rates2$max2
rate2_min3 = rates2$min3
rate2_max3 = rates2$max3
Since I am not aware how many rates I will be dealing with (it can be 2, 4, or
10), hence I can't hardcode my input like I had defined above.
So with the guidance of R - help, I had rerwitten the input code as given below.
# NEW CODE
n = 2 # It gives me no of rates
X = paste('rate', 1:n, '.csv', sep = ' ' )
# the output is "rate1.csv" " rate2.csv"
newrate = lapply(X, read.csv)
# the output is as given below
newrate
[ [1] ]
min1 max1 min2 max2 min3 max3
1 1.05 1.30 1.30 1.65 1.65 1.99
[ [2] ]
min1 max1 min2 max2 min3 max3
1 2.05 2.30 2.30 2.65 2.65 2.99
## _________________________________________________
## PROBLEM
# PROBLEM - A
If I apply the command (As given in OLD CODE above)
rate1_min1 = rates1$min1
# The output which I get is
min1
1 1.05
On similar lines I define
rate1_min1 = newrate[[1]] [1]
# The output which I get is
min1
1 1.05
rate1_max1 = newrate[[1]] [2]
will give me output
max1
1 1.30
and so on.
So I will have to define it 12 times to assign the respective min and max
values.
My problem is instead of hard coding like I have done in OLD CODE, I need to
assign these respective input values using a loop or some other way (as it
becomes cumbersome if no of rates exceed say 10 and also since I am not aware
how many rate files i will be dealing with) so that in the end I should be able
to assign (say) the values like
rate1_min1 = 1.05
rate1_max1 = 1.30
rate1_min2 = 1.30
rate1_max2 = 1.65
rate1_min3 = 1.65
rate1_max3 = 1.99
rate2_min1 = 2.05
rate2_max1 = 2.30
rate2_min2 = 2.30
rate2_max2 = 2.65
rate2_min3 = 2.65
rate2_max3 = 2.99
If there are say 10 rates, then this will continue till
rate10_min1 = ......
rate10_max1 = ......
.........
so on.
## ________________________________________________________
# PROBLEM - B
# Suppose Rij = ith Rate and jth range. (There are 3 ranges i.e. j= 3).
data_label = expand.grid(c("R11", "R12", "R13"), c("R21", "R23", "R23"))
# gives the output like
data_label
Var1 Var2
1 R11 R21
2 R12 R21
3 R13 R21
4 R11 R22
5 R12 R22
6 R13 R22
7 R11 R23
8 R12 R23
9 R13 R23
If instead of two rates, suppose there are three rates, I will have to modify
my code as
data_label = expand.grid(c("R11", "R12", "R13"), c("R21", "R23", "R23"),
c("R31", "R32", "R33"))
If I define say
n = no_of_rates # I will be taking this figure from some otehr csv file.
My PROBLEM = B is how do I write the above code pertaining to data_label in
terms of 'n'.
## _______________________________________________________________________
I understand I am giving the impression as I need some kind of spoon feeding. I
am giving below the original complete OLD CODE (consisting of 4 rates) I had
written in the beginning. Problem with me is to modify it for varaible number
of rates.
I sincerely apologize for the inconvenience I have caused so far by asking
queries. Please guide me.
Regards
Maithili
# ____________________________________________
MY ORIGINAL CODE (Actual HARD CODE)
## ONS - PPA
# INPUT
rate_1 = read.csv('rate1_range.csv')
rate_2 = read.csv('rate2_range.csv')
rate_3 = read.csv('rate3_range.csv')
rate_4 = read.csv('rate4_range.csv')
prob_1 = read.csv('rate1_probs.csv')
prob_2 = read.csv('rate2_probs.csv')
prob_3 = read.csv('rate3_probs.csv')
prob_4 = read.csv('rate4_probs.csv')
rate1_min_1 = rate_1$rate1_range1_min
rate1_max_1 = rate_1$rate1_range1_max
rate1_min_2 = rate_1$rate1_range2_min
rate1_max_2 = rate_1$rate1_range2_max
rate1_min_3 = rate_1$rate1_range3_min
rate1_max_3 = rate_1$rate1_range3_max
rate2_min_1 = rate_2$rate2_range1_min
rate2_max_1 = rate_2$rate2_range1_max
rate2_min_2 = rate_2$rate2_range2_min
rate2_max_2 = rate_2$rate2_range2_max
rate2_min_3 = rate_2$rate2_range3_min
rate2_max_3 = rate_2$rate2_range3_max
rate3_min_1 = rate_3$rate3_range1_min
rate3_max_1 = rate_3$rate3_range1_max
rate3_min_2 = rate_3$rate3_range2_min
rate3_max_2 = rate_3$rate3_range2_max
rate3_min_3 = rate_3$rate3_range3_min
rate3_max_3 = rate_3$rate3_range3_max
rate4_min_1 = rate_4$rate4_range1_min
rate4_max_1 = rate_4$rate4_range1_max
rate4_min_2 = rate_4$rate4_range2_min
rate4_max_2 = rate_4$rate4_range2_max
rate4_min_3 = rate_4$rate4_range3_min
rate4_max_3 = rate_4$rate4_range3_max
## _______________________________________________________________________________________________________________________________________________________________________________________________
## Rij : ith rate jth range i.e. R23 represents 3rd Range of Rate 2
data_lab <- expand.grid(c("R11", "R12", "R13"), c("R21", "R22", "R23"),
c("R31", "R32", "R33"), c("R41", "R42", "R43"))
range_prob <- list()
range_prob[[1]] <- c(prob_1$rate1_prob1,prob_1$rate1_prob2,
prob_1$rate1_prob3)
range_prob[[2]] <- c(prob_2$rate2_prob1,prob_2$rate2_prob2,
prob_2$rate2_prob3)
range_prob[[3]] <- c(prob_3$rate3_prob1,prob_3$rate3_prob2,
prob_3$rate3_prob3)
range_prob[[4]] <- c(prob_4$rate4_prob1,prob_4$rate4_prob2,
prob_4$rate4_prob3)
## COMPUTATIONS
pdf <- expand.grid(range_prob)
data_lab$probs <- apply(pdf, 1, prod)
joint_probs = xtabs(probs ~ Var1 + Var2+Var3+Var4, data = data_lab)
write.csv(data.frame(joint_probs), 'joint_probs.csv', row.names = FALSE)
#________________________________________________________________________________________________________________________________________________________________________________________________________________________
ONS = read.csv('joint_probs.csv')
Names <- apply(ONS[, -5], 1, paste, collapse = ' ')
write.csv(data.frame(rate_combination = Names, Probability = ONS[, 5]),
'prob_table.csv', row.names = FALSE)
#
______________________________________________________________________________________________________________________________________________________________________________________________________________________________
N = 10000 # No of random numbers to be generated for each of the 4 rates
no_of_instances = round(N*(read.csv('prob_table.csv')$Probability))
# # THE ACTUAL PROBLEM - TO GENERATE 10000 random numbers in proportion and
depending on the range of teh given rate in teh combination
write.csv(data.frame(rate1_range = ONS[,1], rate2_range = ONS[,2], rate3_range
= ONS[,3], rate4_range = ONS[,4], no_of_instances), 'Final Table.csv',
row.names = FALSE)
## Random number generation for Rate1
## -----------------------------------
HM = read.csv('Final Table.csv')
rate1_rates=rep(c("R11","R12","R13"),27)
rate1_rates
number = HM$no_of_instances
rate1_combination=rep(rate1_rates,number)
rate1_combination
interval.min=ifelse(rate1_combination=="R11",1.05,ifelse(rate1_combination=="R12",1.3,1.6))
interval.max=ifelse(rate1_combination=="R11",1.3,ifelse(rate1_combination=="R12",1.6,1.99))
rand.nums_rate1=runif(rate1_combination,min=interval.min,max=interval.max)
rate1_series = data.frame(rate1_combination,rand.nums_rate1)
write.csv(data.frame(rate1_series), 'rate1_series.csv', row.names = FALSE)
## Random number generation for Rate2
rate2_rates=rep(c("R21","R21","R21","R22","R22","R22","R23","R23","R23"),9)
rate2_rates
number = HM$no_of_instances
rate2_combination=rep(rate2_rates,number)
rate2_combination
interval.min=ifelse(rate2_combination=="R21",2.05,ifelse(rate2_combination=="R22",2.3,2.6))
interval.max=ifelse(rate2_combination=="R21",2.3,ifelse(rate2_combination=="R22",2.6,2.99))
rand.nums_rate2=runif(rate2_combination,min=interval.min,max=interval.max)
rate2_series = data.frame(rate2_combination,rand.nums_rate2)
write.csv(data.frame(rate2_series), 'rate2_series.csv', row.names = FALSE)
## Random number generation for Rate3
rate3_rates=rep(c("R31","R31","R31","R31","R31","R31","R31","R31","R31","R32","R32","R32","R32","R32","R32","R32","R32","R32","R33","R33","R33","R33","R33","R33","R33","R33","R33"),3)
rate3_rates
number = HM$no_of_instances
rate3_combination=rep(rate3_rates,number)
rate3_combination
interval.min=ifelse(rate3_combination=="R31",3.05,ifelse(rate3_combination=="R32",3.3,3.6))
interval.max=ifelse(rate3_combination=="R21",3.3,ifelse(rate3_combination=="R32",3.6,3.99))
rand.nums_rate3=runif(rate3_combination,min=interval.min,max=interval.max)
rate3_series = data.frame(rate3_combination,rand.nums_rate3)
write.csv(data.frame(rate3_series), 'rate3_series.csv', row.names = FALSE)
## Random number generation for Rate4
rate4_rates=rep(c("R41","R42","R43"), each = 27)
rate4_rates
number = HM$no_of_instances
rate4_combination=rep(rate4_rates,number)
rate4_combination
interval.min=ifelse(rate4_combination=="R41",4.05,ifelse(rate4_combination=="R42",4.3,4.6))
interval.max=ifelse(rate4_combination=="R41",4.3,ifelse(rate4_combination=="R42",4.6,4.99))
rand.nums_rate4=runif(rate4_combination,min=interval.min,max=interval.max)
rate4_series = data.frame(rate4_combination,rand.nums_rate4)
write.csv(data.frame(rate4_series), 'rate4_series.csv', row.names = FALSE)
## Generation of combination names as per the no of random numbers generated
PP = read.csv('prob_table.csv')$rate_combination
PPP = rep(PP, number)
write.csv(data.frame(rate_combination = PPP, rate1 =
round(read.csv('rate1_series.csv')$rand.nums, digits = 4), rate2 =
round(read.csv('rate2_series.csv')$rand.nums, digits = 4), rate3 =
round(read.csv('rate3_series.csv')$rand.nums, digits = 4), rate4 =
round(read.csv('rate4_series.csv')$rand.nums, digits = 4)),
'Main_Random_Numbers.csv', row.names = FALSE)
_________________________________________________________________________________________________________________________
## END of code : ONS-PPA
[[elided Yahoo spam]]
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
