在R中的值数组上返回函数输出

我是一个初学R程序员，在多元数组问题中挣扎。在R中的值数组上返回函数输出

我试图输入的4个参数值的数组，说a=1:10, b=1:10, p=1:10, q=1:10，成一个函数y=f(x|a, b, p, q)，基于我的数据集，X计算的y值，和给定的4个参数的每一个可能的组合[(a=1,b=1,p=1,q=1),(a=2,b=1,p=1,q=1),...,(a=10,b=1,p=1,q=1),...,(a=10,b=10,p=10,q=10)] = 10^4 = 10,000个可能的组合，因此有10,000个y值。

理想情况下，我希望输出为数组格式，然后我可以在R中绘制图形，从而允许将每个参数绘制为单独的轴。

如果有人能指出我在正确的方向，将不胜感激！

感谢，

罗伯特

来源

2011-04-06 Robert

罗伯特，没有不敬的意图，但你并不是真的在问一个问题。需要做很多步骤来完成你所要求的事情，从你的问题中不可能知道你挂在哪里。您现在可以抓取的最好的初学者R技能是将问题分解成多个步骤。然后尝试做每一步。如果你不能做到这一步，就把它分解成一个问题，然后展示一个你想要什么和尝试过什么的简单例子。 – 2011-04-06 21:23:40

同意JD。在试图编写一个答案时，我意识到，我将最终编码每一步。 – Maiasaura 2011-04-06 22:32:26

是的，对不起，应该更清楚。我试图找出四参数概率分布（第二类的广义Beta）的负对数似然，以帮助我找到我的保险索赔数据的最佳拟合（=最小负对数似然值）。在SAS和R中有各种迭代算法，但是我在改变初始参数值时一直有问题产生稳定的结果。所以我要用“眼球”方法来找出最佳的初始参数。 – Robert 2011-04-12 13:18:46

我同意JD长，该请求不明确，允许一个最终的答案，但有一个问题的第一部分：

all.comb.dfrm <- expand.grid(a=1:10, b=1:10, p=1:10, q=1:10) 
all.comb.dfrm$Y <- with(all.comb.dfrm, f(a,b,p,q))

来源

2011-04-06 21:47:22

我从字面上用不同的对象名称输入了相同的答案。揍我一拳！ +1 – Chase 2011-04-06 22:06:49

谢谢，JD，非常漂亮！ – Robert 2011-04-12 13:14:53

在R中的值数组上返回函数输出

回答

相关问题