Guest User

Untitled

a guest
Apr 27th, 2018
159
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
QBasic 3.37 KB | None | 0 0
  1. prompt1<-menu(c("PDF","CDF"), graphics = TRUE, title = "CDF or PDF?")
  2. # Distribution TYPE
  3. disType<-menu(c("Binomial","Geometric", "Hypergeometric","Poisson"), graphics = TRUE, title = "Which
  4. distribution?")
  5.  
  6. IF(prompt1==1&disType==4){
  7.     lam<-winDialogString("Average"," ")
  8.     lam<-AS.numeric(lam)
  9.    
  10.     Exact<-winDialogString("amount?", " ")
  11.     Exact<-AS.numeric(Exact)
  12.     pPoisson<- round(dpois(Exact,lam),2)
  13.     PRINT(pPoisson)
  14.  
  15. }ELSE IF(prompt1==2&disType==4){
  16.  
  17.     lam<-winDialogString("Average"," ")
  18.     lam<-AS.numeric(lam)
  19.  
  20.     Max<-winDialogString("Amount?", " ")
  21.     Max<-AS.numeric(Max)
  22.    
  23.     pPoisson<-round(ppois(Max,lam),2)
  24.     PRINT(pPoisson)
  25.  
  26. }ELSE{
  27. lower<-winDialogString("lower?", " ") #enters a STRING value FOR "p"
  28. lower<-AS.numeric(lower)
  29.  
  30. higher<-winDialogString("higher?", " ") #enters a STRING value FOR "p"
  31. higher<-AS.numeric(higher)
  32.  
  33.  
  34.  
  35. range <- lower:higher
  36.  
  37. # PDF
  38. IF(prompt1==1){
  39.     IF(disType==1){ #binomial
  40.         p<-winDialogString("Probability of success?", " ") #enters a STRING value FOR "p"
  41.         p<-AS.numeric(p)
  42.         n<-winDialogString("Number of trials?", " ") #enters a STRING value FOR "p"
  43.         n<-AS.numeric(n)
  44.  
  45.         pbinomial<-round(dbinom(range, n, p),2)
  46.         PRINT(pbinomial)
  47.     }
  48.     ELSE IF(disType==2) { #geometric
  49.         p<-winDialogString("Probability of success?", " ") #enters a STRING value FOR "p"
  50.         p<-AS.numeric(p)
  51.         pGeometric<-round(dgeom(range,p),2)
  52.         PRINT(pGeometric)
  53.         # Graph
  54.         plot(range+1, pGeometric, xlab = "X = Number of Trials", ylab = "P(X=x)",
  55.         TYPE = "h", main = "First Ready Terminal, p = 0.95")
  56.     }
  57.     ELSE IF(disType==3) { #Hypergeometric
  58.         fav<-winDialogString("No. of favourable outcomes?", " ") #enters a STRING value FOR "p"
  59.         fav<-AS.numeric(fav)
  60.        
  61.         nonFav<-winDialogString("No. of non-favourable outcomes?", " ") #enters a STRING value FOR "p"
  62.         nonFav<-AS.numeric(nonFav)
  63.  
  64.         chosen<-winDialogString("No. of chosen?", " ") #enters a STRING value FOR "p"
  65.         chosen<-AS.numeric(chosen)
  66.  
  67.         phyp<-round(dhyper(range,fav,nonFav,chosen),2)
  68.         PRINT(phyp)
  69.         # Graph
  70.         plot(range+1, phyp, xlab = "X = Number of Trials", ylab = "P(X=x)",
  71.         TYPE = "s", main = "First Ready Terminal, p = 0.95")
  72.     }
  73.    
  74. # CDF
  75. } ELSE {
  76.     IF(disType==1){ #binomial
  77.         p<-winDialogString("Probability of success?", " ") #enters a STRING value FOR "p"
  78.         p<-AS.numeric(p)
  79.         n<-winDialogString("Number of trials?", " ") #enters a STRING value FOR "p"
  80.         n<-AS.numeric(n)
  81.  
  82.         cbinomial<-round(pbinom(range, n, p),2)
  83.         PRINT(cbinomial)
  84.     }
  85.     ELSE IF(disType==2) { #geometric
  86.         p<-winDialogString("Probability of success?", " ") #enters a STRING value FOR "p"
  87.         p<-AS.numeric(p)
  88.         cGeometric<-round(pgeom(range,p),2)
  89.         PRINT(cGeometric)
  90.         # Graph
  91.         plot(range+1, cGeometric, xlab = "X = Number of Trials", ylab = "P(X=x)",
  92.         TYPE = "s", main = "First Ready Terminal, p = 0.95")
  93.     }
  94.     ELSE IF(disType==3) { #Hypergeometric
  95.         fav<-winDialogString("No. of favourable outcomes?", " ") #enters a STRING value FOR "p"
  96.         fav<-AS.numeric(fav)
  97.        
  98.         nonFav<-winDialogString("No. of non-favourable outcomes?", " ") #enters a STRING value FOR "p"
  99.         nonFav<-AS.numeric(nonFav)
  100.  
  101.         chosen<-winDialogString("No. of chosen?", " ") #enters a STRING value FOR "p"
  102.         chosen<-AS.numeric(chosen)
  103.  
  104.         chyp<-round(phyper(range,fav,nonFav,chosen),2)
  105.         PRINT(chyp)
  106.         # Graph
  107.         plot(range+1, chyp, xlab = "X = Number of Trials", ylab = "P(X=x)",
  108.         TYPE = "s", main = "First Ready Terminal, p = 0.95")
  109.     }
  110. }}
Add Comment
Please, Sign In to add comment