[R] 3D surface plot with wireframe or persp?
Uwe Ligges
ligges at statistik.tu-dortmund.de
Wed Apr 21 19:06:10 CEST 2010
You cannot specify a dfunction for z, but need to compute the values in
the matrix yourself as in:
persp(data_for_time, data_for_s,
outer(data_for_time, data_for_s, plot_R_i_3d))
Uwe Ligges
On 20.04.2010 17:35, Jin wrote:
>
> Hello Dear,
>
> I have a function, like z=f(x,y), and try a surface plot with this function.
> But, on the reference of "wireframe" requires data option, so I generated x
> and y, and computed z with them. But, still I have a problem to draw a
> surface plot. The code and errors are
>
> ##################################################
> mle_beta0=64.43707;
> mle_beta1=-24365.16;
>
> # generating for the requirement of "wireframe"
> data_for_c = runif(1000,30,85); # range is (30,50)
> data_for_s = sort(1/(273+data_for_c));
> data_for_time = sort(runif(1000,0,100000)) # range is (0,100000)
> data_for_R = exp((-exp(mle_beta0+mle_beta1*data_for_s))*data_for_time)
> data_all = cbind(data_for_s,data_for_time,data_for_R)
>
> # function: plot_R_i_3d = f(data_for_time,data_for_s)
> plot_R_i_3d = function(data_for_time,data_for_s)
> {
>
> R_i = exp((-exp(mle_beta0+mle_beta1*data_for_s))*data_for_time);
> return(R_i)
> }
>
> # tried 1) persp or 2) wireframe
>
> persp(data_for_time,data_for_s,plot_R_i_3d)
>
> ==> Error in min(x, na.rm = na.rm) : invalid 'type' (list) of argument
>
> wireframe(formula=data_for_R ~ data_for_s*data_for_time,data=data_all)
>
> ==> Error in nrow(x) :
> element 1 is empty;
> the part of the args list of 'dim' being evaluated was:
> (x)
> ###################################################
>
> I have two questions:
>
> 1) which one is better to plot in this case - "persp" and "wireframe", Is
> there another suggestion?
> 2) Is it possible to plot with only a function?
> 3) If not, what is wrong in the code above?
>
> Many thank you in advance,
>
> Jin
>
>
More information about the R-help
mailing list