本文整理汇总了Python中random.lognormvariate函数的典型用法代码示例。如果您正苦于以下问题:Python lognormvariate函数的具体用法?Python lognormvariate怎么用?Python lognormvariate使用的例子?那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。
在下文中一共展示了lognormvariate函数的16个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于我们的系统推荐出更棒的Python代码示例。
示例1: random_year
def random_year(self):
s = random.lognormvariate(self.stocks_mean, self.stocks_stddev) - 1
s = Decimal(s)
b = random.lognormvariate(self.bonds_mean, self.bonds_stddev) - 1
b = Decimal(b)
return AnnualChange(year=self.year, stocks=s, bonds=b, inflation=Decimal(0))
开发者ID:hoostus,项目名称:prime-harvesting,代码行数:7,代码来源:montecarlo.py
示例2: sample
def sample(self):
result = self._scale * random.lognormvariate(self._mu, self._sigma)
if self._maximum is not None and result > self._maximum:
for _ignore in range(10):
result = self._scale * random.lognormvariate(self._mu, self._sigma)
if result <= self._maximum:
break
else:
raise ValueError("Unable to generate LogNormalDistribution sample within required range")
return result
开发者ID:anemitz,项目名称:calendarserver,代码行数:10,代码来源:stats.py
示例3: produceReads
def produceReads(self,readDensity = 1,read_length = 50):
"""
Produces uniform reads by walking iteratively over self.regions
"""
mRNA = self.getTotalmRNA()
spanning_length = self.getRegionSpanningLength()
n = spanning_length['total'] - read_length + 1
j = 0
k = 0
for i in range(n):
# "alpha is playing the role of k and beta is playing the role of theta"
dd = max(0,int(round(random.lognormvariate(math.log(readDensity),0.5))))# Notice this is NOT a binomial distribution!!
for d in range(dd):
sequence = mRNA[i:i+read_length]
if(random.randint(0,1) == 0):
strand = 0
else:
strand = 16
flag = strand + 0
print "read_"+str(j)+"."+str(i)+"."+str(d)+"\t"+str(flag)+"\t"+self.chromosome+"\t"+str(self.regions[j].start + k)+"\t60\t"+self.getMappingString(read_length,j,k)+"\t*\t0\t0\t"+str(sequence.upper())+"\t*"
spanning_length['iter'][j] -= 1
if(k >= self.regions[j].getSpanningLength()-1):
j += 1
k = 0
else:
k += 1
开发者ID:ErasmusMC-Bioinformatics,项目名称:samtools_parallel_mpileup_galaxy_wrapper,代码行数:33,代码来源:generate_reads.py
示例4: create_server
async def create_server(self, worker, protocol_factory, address=None,
sockets=None, idx=0):
'''Create the Server which will listen for requests.
:return: a :class:`.TcpServer`.
'''
cfg = self.cfg
max_requests = cfg.max_requests
if max_requests:
max_requests = int(lognormvariate(log(max_requests), 0.2))
server = self.server_factory(
protocol_factory,
loop=worker._loop,
max_requests=max_requests,
keep_alive=cfg.keep_alive,
name=self.name,
logger=self.logger,
server_software=cfg.server_software,
cfg=cfg,
idx=idx
)
for event in ('connection_made', 'pre_request', 'post_request',
'connection_lost'):
callback = getattr(cfg, event)
if callback != pass_through:
server.event(event).bind(callback)
await server.start_serving(
sockets=sockets,
address=address,
backlog=cfg.backlog,
sslcontext=self.sslcontext()
)
return server
开发者ID:quantmind,项目名称:pulsar,代码行数:33,代码来源:__init__.py
示例5: test_lognormal
def test_lognormal(self):
"""Statistics.lognormal test"""
import random
import Biskit.gnuplot as gnuplot
import Biskit.hist as H
cr = []
for i in range( 10000 ):
## Some random values drawn from the same lognormal distribution
alpha = 1.5
beta = .7
x = 10.
R = [ random.lognormvariate( alpha, beta ) for j in range( 10 ) ]
cr += [ logConfidence( x, R )[0] ]
ca = logArea( x, alpha, beta )
if self.local:
gnuplot.plot( H.density( N.array(cr) - ca, 100 ) )
globals().update( locals() )
self.assertAlmostEqual( ca, 0.86877651432955771, 7)
开发者ID:ostrokach,项目名称:biskit,代码行数:27,代码来源:lognormal.py
示例6: find_distributions
def find_distributions(self, n, c, source):
"""
Finds distribution functions
"""
if source[c][n] == 'NoArrivals':
return lambda : 'Inf'
if source[c][n][0] == 'Uniform':
return lambda : uniform(source[c][n][1], source[c][n][2])
if source[c][n][0] == 'Deterministic':
return lambda : source[c][n][1]
if source[c][n][0] == 'Triangular':
return lambda : triangular(source[c][n][1], source[c][n][2], source[c][n][3])
if source[c][n][0] == 'Exponential':
return lambda : expovariate(source[c][n][1])
if source[c][n][0] == 'Gamma':
return lambda : gammavariate(source[c][n][1], source[c][n][2])
if source[c][n][0] == 'Lognormal':
return lambda : lognormvariate(source[c][n][1], source[c][n][2])
if source[c][n][0] == 'Weibull':
return lambda : weibullvariate(source[c][n][1], source[c][n][2])
if source[c][n][0] == 'Custom':
P, V = zip(*self.parameters[source[c][n][1]])
probs = list(P)
cum_probs = [sum(probs[0:i+1]) for i in range(len(probs))]
values = list(V)
return lambda : self.custom_pdf(cum_probs, values)
if source[c][n][0] == 'Empirical':
if isinstance(source[c][n][1], str):
empirical_dist = self.import_empirical_dist(source[c][n][1])
return lambda : choice(empirical_dist)
return lambda : choice(source[c][n][1])
return False
开发者ID:gitter-badger,项目名称:Ciw,代码行数:32,代码来源:simulation.py
示例7: _trigger_func
def _trigger_func(self, stream_id):
while True:
length = random.randint(1,20)
c = [random.uniform(0.0,75.0) for i in xrange(length)]
t = [random.uniform(-1.7, 21.0) for i in xrange(length)]
p = [random.lognormvariate(1,2) for i in xrange(length)]
lat = [random.uniform(-90.0, 90.0) for i in xrange(length)]
lon = [random.uniform(0.0, 360.0) for i in xrange(length)]
tvar = [self.last_time + i for i in xrange(1,length+1)]
self.last_time = max(tvar)
ctd_packet = ctd_stream_packet(stream_id=stream_id,
c=c, t=t, p=p, lat=lat, lon=lon, time=tvar)
log.warn('SimpleCtdPublisher sending %d values!' % length)
self.publisher.publish(ctd_packet)
time.sleep(2.0)
开发者ID:daf,项目名称:coi-services,代码行数:27,代码来源:ctd_stream_publisher.py
示例8: generateLogNormalRVs
def generateLogNormalRVs(mu, sigma, count, histDelta):
rVec = numpy.zeros(count)
rMax = 0.
for ii in range(count):
rVec[ii] = random.lognormvariate(mu, sigma)
if (rMax < rVec[ii]):
rMax = rVec[ii]
# build the histogram ...
deltaR = histDelta
rMinHist = 0.
rMaxHist = rMax + 2. * deltaR
numBins = (rMaxHist - rMinHist) / deltaR
numBins = int(numBins) + 2
rHist = numpy.zeros(numBins)
# print ' making histogram ... ', deltaR, rMaxHist
for ii in range(count):
iBin = int((rVec[ii] - rMinHist) / deltaR + 0.0001)
rHist[iBin] += 1
for iBin in range(numBins):
rHist[iBin] /= float(count)
return (rVec, rHist)
开发者ID:cancerregulome,项目名称:gidget,代码行数:27,代码来源:miscMath.py
示例9: find_service_time
def find_service_time(self, n, c):
"""
Finds the service time function
"""
if self.mu[c][n][0] == 'Uniform':
return lambda : uniform(self.mu[c][n][1], self.mu[c][n][2])
if self.mu[c][n][0] == 'Deterministic':
return lambda : self.mu[c][n][1]
if self.mu[c][n][0] == 'Triangular':
return lambda : triangular(self.mu[c][n][1], self.mu[c][n][2], self.mu[c][n][3])
if self.mu[c][n][0] == 'Exponential':
return lambda : expovariate(self.mu[c][n][1])
if self.mu[c][n][0] == 'Gamma':
return lambda : gammavariate(self.mu[c][n][1], self.mu[c][n][2])
if self.mu[c][n][0] == 'Normal':
return lambda : gauss(self.mu[c][n][1], self.mu[c][n][2])
if self.mu[c][n][0] == 'Lognormal':
return lambda : lognormvariate(self.mu[c][n][1], self.mu[c][n][2])
if self.mu[c][n][0] == 'Weibull':
return lambda : weibullvariate(self.mu[c][n][1], self.mu[c][n][2])
if self.mu[c][n][0] == 'Custom':
P, V = zip(*self.parameters[self.mu[c][n][1]])
probs = list(P)
cum_probs = [sum(probs[0:i+1]) for i in range(len(probs))]
values = list(V)
return lambda : self.custom_pdf(cum_probs, values)
return False
开发者ID:drvinceknight,项目名称:Ciw,代码行数:28,代码来源:simulation.py
示例10: random
def random(self, start, minVal, maxVal, candidate = None):
for z in xrange(config.rand_retries):
v=int(start*random.lognormvariate(0,1)+.5)
#logging.debug("lognorm: start=%d, v=%d", start, v)
if v>=minVal and v<=maxVal and start!=v:
return v
raise MutateFailed("lognorm random gen failed")
开发者ID:AikidoGuy,项目名称:petabricks,代码行数:7,代码来源:mutators.py
示例11: visit_medic
def visit_medic(self):
"""
Lognormal with :
logarithmic mean: 1.024
logarithmic std dev: 0.788
"""
name = 'medic'
if self.debug:
print self.name, "at ",name, simpy.now()
arrive = simpy.now()
yield simpy.request, self, self.resources[name]
wait = simpy.now() - arrive
self.monitors[name].observe(wait)
time = random.lognormvariate(mu=1.024, sigma=0.788)
tib = self.numberOfForms * time
yield simpy.hold,self,tib
yield simpy.release, self, self.resources[name]
p = random.random()
if p < 0.99:
for i in self.visit_dispenser():
yield i
else:
for i in self.visit_exit():
yield i
开发者ID:ivihernandez,项目名称:staff-allocation,代码行数:26,代码来源:Customer.py
示例12: next
def next(self):
change = random.lognormvariate(0, self.sigma)
if self.current + change > 0 and random.choice((True, False)):
change = -change
self.current += change
delay()
yield self.current
开发者ID:goose42,项目名称:something,代码行数:7,代码来源:price_source.py
示例13: _get_new_ctd_packet
def _get_new_ctd_packet(self, length):
rdt = RecordDictionaryTool(stream_definition_id=self.stream_def._id)
#Explicitly make these numpy arrays...
c = numpy.array([random.uniform(0.0,75.0) for i in xrange(length)])
t = numpy.array([random.uniform(-1.7, 21.0) for i in xrange(length)])
p = numpy.array([random.lognormvariate(1,2) for i in xrange(length)])
lat = numpy.array([random.uniform(-90.0, 90.0) for i in xrange(length)])
lon = numpy.array([random.uniform(0.0, 360.0) for i in xrange(length)])
h = numpy.array([random.uniform(0.0, 360.0) for i in xrange(length)])
start_time = ntplib.system_to_ntp_time(time.time()) - (length + 1)
tvar = numpy.array([start_time + i for i in xrange(1,length+1)])
rdt['time'] = tvar
rdt['lat'] = lat
rdt['lon'] = lon
rdt['temp'] = t
rdt['conductivity'] = c
rdt['pressure'] = p
# rdt['coordinates'] = rdt0
# rdt['data'] = rdt1
g = rdt.to_granule(data_producer_id=self.id)
return g
开发者ID:Bobfrat,项目名称:coi-services,代码行数:28,代码来源:ctd_stream_publisher.py
示例14: _get_new_ctd_packet
def _get_new_ctd_packet(self, stream_id, length):
rdt = RecordDictionaryTool(taxonomy=tx)
#Explicitly make these numpy arrays...
c = numpy.array([random.uniform(0.0,75.0) for i in xrange(length)])
t = numpy.array([random.uniform(-1.7, 21.0) for i in xrange(length)])
p = numpy.array([random.lognormvariate(1,2) for i in xrange(length)])
lat = numpy.array([random.uniform(-90.0, 90.0) for i in xrange(length)])
lon = numpy.array([random.uniform(0.0, 360.0) for i in xrange(length)])
h = numpy.array([random.uniform(0.0, 360.0) for i in xrange(length)])
tvar = numpy.array([self.last_time + i for i in xrange(1,length+1)])
self.last_time = max(tvar)
log.warn('Got time: %s' % str(tvar))
log.warn('Got t: %s' % str(t))
rdt['time'] = tvar
rdt['lat'] = lat
rdt['lon'] = lon
rdt['height'] = h
rdt['temp'] = t
rdt['cond'] = c
rdt['pres'] = p
# rdt['coordinates'] = rdt0
# rdt['data'] = rdt1
g = build_granule(data_producer_id=stream_id, taxonomy=tx, record_dictionary=rdt)
return g
开发者ID:pombredanne,项目名称:coi-services,代码行数:31,代码来源:ctd_stream_publisher.py
示例15: initRandomParticle
def initRandomParticle(self, systemRadius, maxSize, systemCenter):
position = [random.randrange(systemCenter[0],systemRadius),random.randrange(systemCenter[1],systemRadius)] #this creates a square system
velocity = [random.randrange(-50,50),random.randrange(-50,50)]
mass = random.lognormvariate(1,6)
radius = math.pow(mass,0.25)
color = (random.randrange(0,255),random.randrange(0,255),random.randrange(0,255),0)
return self(position,velocity,radius,mass,color)
开发者ID:lalaithion,项目名称:nBodySim,代码行数:7,代码来源:ParticleClass.py
示例16: lognormal
def lognormal(mean=None, sigma=None, n=1, minimum=None, maximum=None):
"""
Generate a list with a lognormal distribution.
Option - setting mininum and maximum rejects and regenerates
any value outside the truncation range.
Setting minimum and maximum to None or 'None' makes the range
unbounded.
"""
if minimum is None or minimum is 'None':
minimum = None
else:
minimum = float(minimum)
if maximum is None or maximum is 'None':
maximum = None
else:
maximum = float(maximum)
mean = float(mean)
sigma = float(sigma)
answer = []
while len(answer) < n:
x = random.lognormvariate(mu=mean, sigma=sigma)
if x > minimum and (x < maximum or (maximum is None)):
answer.append(x)
return answer
开发者ID:dynaryu,项目名称:eqrm,代码行数:27,代码来源:distributions.py
注:本文中的random.lognormvariate函数示例由纯净天空整理自Github/MSDocs等源码及文档管理平台,相关代码片段筛选自各路编程大神贡献的开源项目,源码版权归原作者所有,传播和使用请参考对应项目的License;未经允许,请勿转载。 |
请发表评论