Re: Evaluate, Replace, and Hold
- To: mathgroup at smc.vnet.net
- Subject: [mg77178] Re: Evaluate, Replace, and Hold
- From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
- Date: Tue, 5 Jun 2007 06:44:08 -0400 (EDT)
- Organization: The Open University, Milton Keynes, UK
- References: <f40gib$65a$1@smc.vnet.net>
grenander at gmail.com wrote:
> This is a newbie question. I have a parametric form of a model and I
> would like to evaluate it with different parameters on a list of
> values.
>
> model = a x + b;
>
> Different sets of parameters:
>
> p1 = {a->1, b->0};
> p2 = {a->1, b->1};
>
> x for which I'd like to evaluate the model:
>
> pts = Range[10];
>
> I would like to define a function to evaluate the model at pts using
> any of the sets of paramters:
>
> f[pts, p1]
> f[pts, p2]
>
> I would also like it to work for single points:
>
> f[1, p1]
> f[2, p1]
>
> I've tried various forms for a function definition that dont work. I
> can do the following for a particular set of parameters:
>
> Function[{x}, Evaluate[model /. p1]] /@ pts
>
> but how can I turn it into a function of also the parameters. The
> below doesn't work, for example:
>
> f[pts_, p_]:=Function[{x}, Evaluate[model /. p] ] /@ pts
>
> What I know is from descending the documentation tree (v6.0) starting
> from the examples in FindFit. I have a fundamental misunderstanding of
> replacement rules and delayed evaluation.
>
> Thanks.
You may want to try something along the lines,
In[1]:=
model[x_] = a*x + b;
p1 = {a -> 1, b -> 0};
p2 = {a -> 1, b -> 1};
pts = Range[10];
f[data_List /; VectorQ[data, NumberQ],
params_List] := (model[#1] /. params & ) /@ data;
f[(data_)?NumberQ, params_List] :=
(model[#1] /. params & )[data]
f[pts, p1]
f[pts, p2]
f[1, p1]
f[2, p2]
Out[7]=
{1, 2, 3, 4, 5, 6, 7, 8, 9, 10}
Out[8]=
{2, 3, 4, 5, 6, 7, 8, 9, 10, 11}
Out[9]=
1
Out[10]=
3
Hope this helps,
Jean-Marc