[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