我在编写收入函数时遇到问题。我得到的信息是:“查尔斯接管了一个卖鲜鱼的市场摊位。每天早上,他从批发商那里购买鱼,价格为每公斤 4 英镑,最高 100 公斤,每公斤 2.5 英镑任何超过 100 公斤的东西。他以每公斤 8 英镑的价格出售鱼。如果在一天结束时有任何鱼仍未售出,他以每公斤 1 英镑的价格将其卖给猫粮厂。查尔斯想知道要卖多少鱼每天购买以最大化他的平均每日利润。”
到目前为止,我将过去 150 天的每日销售数据的给定数据建模为伽马函数,并使用矩方法来查找参数,下面是我使用的代码:
First_moment_Xbar = sum(nfsold)/150
Second_moment = sum((nfsold)^2)/150\
estimator_of_lambda = (First_moment_Xbar)/(Second_moment - (First_moment_Xbar)^2)
estimator_of_k = (First_moment_Xbar)^2 / (Second_moment - (First_moment_Xbar)^2 )
estimator_of_k
estimator_of_lambda
我得到了:
> estimator_of_k
[1] 4.992957
> estimator_of_lambda
[1] 0.05733275
我写了一个购买 n 公斤鱼的成本函数:
cost=function(n)
{
total = 0
if(n<=100)
{
total = n*4
}
else{
total = 100*4 + (n-100)*2.5
}
return(total)
}
现在,当有 d 公斤的潜在销售并且查尔斯买了 n 公斤鱼时,我必须编写一个收入函数,收入(n,d)。
我认为潜在的销售将围绕:
round(rgamma(1, estimator_of_k, estimator_of_lambda), digits = 0)
据我所知,在猫粮之前的鱼类销售将围绕我发现的 Gamma 建模。
请问一下如何找到收入函数,谢谢。