Displaying 4 results from an estimated 4 matches for "probdwn".
2017 Nov 02
2
repeat a function
...function(N) {
ProbUP<- function( a, j, dt) 1 / 6 + ((j ^ 2 * Mfactor(a, dt) ^ 2 + j * Mfactor(a, dt)) / 2) ######### Probability X going up
ProbMID<- function(a, j, dt) 2 / 3 - (j ^ 2 * Mfactor(a, dt) ^ 2) ######## Probability X going middle
ProbDWN<-function( a, j, dt) 1 / 6 + ((j ^ 2 * Mfactor(a, dt) ^ 2 - j * Mfactor(a, dt)) / 2) ####### Probability X going down
TOPProbUP<- function( a, j, dt) 7 / 6 + (j ^ 2 * Mfactor(a, dt) ^ 2 + 3 * j * Mfactor(a, dt)) / 2 ####### Top branch probability going up
TOPProbMID<- fun...
2017 Nov 01
3
repeat a function
...function independently in the global environment and then run the loop with the iF statement, as per below.
for (j in 0:2) {
if (j==0) {
prb["0","1"] <- ProbUP(a,j,dt)
prb["0","0"] <- ProbMID(a,j,dt)
prb["0","-1"] <- ProbDWN(a,j,dt)
}
else {
if (j==jmax) {
prb[paste(j,sep = ""),"1"] <- TOPProbUP(a,j,dt)
prb[paste(j,sep = ""),"0"] <- TOPProbMID(a,j,dt)
prb[paste(j,sep = ""),"-1"] <- TOPProbDWN(a,j,dt)
prb[paste(-j,sep =...
2017 Nov 03
0
repeat a function
..., let them call u and v)
> ProbUP<- function( a, j, dt) 1 / 6 + ((j ^ 2 * Mfactor(a, dt) ^ 2 + j * Mfactor(a,
> dt)) / 2) ######### Probability X going up
> ProbMID<- function(a, j, dt) 2 / 3 - (j ^ 2 * Mfactor(a, dt) ^ 2)
> ######## Probability X going middle
> ProbDWN<-function( a, j, dt) 1 / 6 + ((j ^ 2 * Mfactor(a, dt) ^ 2 - j * Mfactor(a,
> dt)) / 2)
Functions are quite similar, you basically need to find how combination of row and column number results in desired j value. However it could be quite tricky and if your approach works and you do not want...
2017 Nov 02
0
repeat a function
...global
> environment and then run the loop with the iF statement, as per below.
> for (j in 0:2) {
> if (j==0) {
> prb["0","1"] <- ProbUP(a,j,dt)
> prb["0","0"] <- ProbMID(a,j,dt)
> prb["0","-1"] <- ProbDWN(a,j,dt)
> }
> else {
> if (j==jmax) {
> prb[paste(j,sep = ""),"1"] <- TOPProbUP(a,j,dt)
> prb[paste(j,sep = ""),"0"] <- TOPProbMID(a,j,dt)
> prb[paste(j,sep = ""),"-1"] <- TOPProbDWN(a,j...