Question

(+30) Write a python program per the following specifications:

1. Populate an array(list) of size 50 randomly with only integers 0 and 1
2. Repeat step1 n = 1000 times using either a while loop or a for loop

At this point you should have a total of 50000 observations

3. Display the number of 0’s (use the count() function from prior labs)
4. Display the number of 1’s (use the count() function from prior labs)
5. Using the Binomial distribution formulas
   1. Display the expected mean (average) of the 1’s == n*size*.5
   2. Calculate and display the standard deviation (SD)

NOTE: sd should be > 100

5. 3. Display the range as follows:   mean + SD mean - SD
6. Answer the question: is the total number of 1’s from 4 above. within the range as calculated in 5.c by printing

                                     ‘Yes’ if it is within the range calculated in 5c

                                     ‘No’   if it is not

7. Display the n*size (== 50000) integers 50 integers per line That is, display the integers in myList 50 per line for 1000 lines See sample output
8. Display the average of all n*size == 50000 integers (should be < 1 by adding the answers from 3 and 4 above and divide by n*size

my work:

import random
import math
myArray=[]
a=0
b=1
size= 50
j=0
n=1000
#populate array with size=50 random integers in range 0-1 and repeat for 1000 times
while(j<n):
for k in range(size):
randNum=random.randint(a,b)
myArray.insert(k,randNum)
j=j+1
#display array size
print("size of my array is ",len(myArray))
#print myArray 50 values per line
k=0
sub=[]
while(k<len(myArray)):
sub=myArray[k:k+50]
print("...",k,sub)
k=k+50

#count no. of ones and zeroes
zeros=0
ones=0
i=0
while(i<len(myArray)):
if(myArray[i]==1):
ones=ones+1
else:
zeros=zeros+1
i=i+1

print("number of 1 is ",ones)
print("number of 0 is ",zeros)

#calculate mean and SD of binomial distribution
#probability of occurence of 1 is p = 1/2 and occurence of 0 is q =1/2
#mean of Binomial distribution = np
#SD(standard deviation) of Binomial distribution = (npq)^1/2
#here n=50000 or length of array
n=len(myArray)
p=0.5
q=0.5
mean= n*p #expected mean of no of ones
sd = math.sqrt(n*p*q)
print("===============================")
print("mean of 1s in My Array: ",mean," sd(standard deviation): ",sd)
Lrange=int(round(mean-sd))
Rrange=int(round(mean+sd))
print("range: mean+-sd: ",Rrange," ",Lrange) #round off the range to nearest integer
#check if no of ones is within range
if(ones>= Lrange and ones<= Rrange):
output="Yes"
else:
output="No"

print("Actual number of ones is within the range mean+-sd ??? Ans is: ", output)

getting error line 11
for k in range(size):
^
IndentationError: expected an indented block please help

Answer 1

import random
import math
myArray=[]
a=0
b=1
size= 50
j=0
n=1000
#populate array with size=50 random integers in range 0-1 and repeat for 1000 times
while(j<n):
for k in range(size):
randNum=random.randint(a,b)
myArray.insert(k,randNum)
j=j+1
#display array size
print("size of my array is ",len(myArray))
#print myArray 50 values per line
k=0
sub=[]
while(k<len(myArray)):
sub=myArray[k:k+50]
print("...",k,sub)
k=k+50

#count no. of ones and zeroes
zeros=0
ones=0
i=0
while(i<len(myArray)):
if(myArray[i]==1):
ones=ones+1
else:
zeros=zeros+1
i=i+1

print("number of 1 is ",ones)
print("number of 0 is ",zeros)

#calculate mean and SD of binomial distribution
#probability of occurence of 1 is p = 1/2 and occurence of 0 is q =1/2
#mean of Binomial distribution = np
#SD(standard deviation) of Binomial distribution = (npq)^1/2
#here n=50000 or length of array
n=len(myArray)
p=0.5
q=0.5
mean= n*p #expected mean of no of ones
sd = math.sqrt(n*p*q)
print("===============================")
print("mean of 1s in My Array: ",mean," sd(standard deviation): ",sd)
Lrange=int(round(mean-sd))
Rrange=int(round(mean+sd))
print("range: mean+-sd: ",Rrange," ",Lrange) #round off the range to nearest integer
#check if no of ones is within range
if(ones>= Lrange and ones<= Rrange):
output="Yes"
else:
output="No"

print("Actual number of ones is within the range mean+-sd ??? Ans is: ", output)

code screenshot:

Expected output: