similar to: 3D plot/surface rotation

Is there any such thing as a program file in R? That is, is it possible to compile and save code for performing data management and analytic tasks for a given project into a single/multiple file(s) (i.e., like a script file (S-Plus), a do file (STATA), or a *.sas file (SAS))? Any words of wisdom or a point in the direction of some documentation would be helpful. If it is not possible to work

Is there any existing function for creating contingency tables that will display counts, row, column, and cell percentages in the same tables....anything similar to crosstabs in S? Marc W. Zodet, MS Health Statistician Center for Cost and Financing Studies Division of Statistical Research and Methods 2101 East Jefferson Street, Suite 500 Rockville, Maryland 20852 Phone: 301-594-7072 Fax:

Is there a way to associate text descriptions with variables in a data.frame? For example... Let's say that in my data.frame I have a variable named var1. var1 represents the responses to the question "When was the last time you saw your physician?" When I tabulate the variable var1 I'd like the output to be a bit more descriptive and contain the more descriptive

I need some guidance from someone who is familiar/has some experience with the survey package. The data that I am using is from the Medical Expenditure Panel Survey (www.meps.ahrq.gov <http://www.meps.ahrq.gov/> ). The STRATA and PSU variables are varstr01 and varpsu01 respectively. When I try to specify them with the svydesign function I get an error message. An excerpt of my session

What is the best way to sort a dataframe? For example, how would I go about sorting a dataframe (with variables V1-V5) by ascending V1, V2 and descending V3 while retaining V4 and V5. Also, is there a relatively easy way by which to re-order my columns? Thank you for your time. Marc Marc Zodet [[alternative HTML version deleted]]

Is there a function in R which performs stepwise estimation in ways similar to SAS/STATA (i.e., allows the analyst to specify the significance levels for removal/addition of terms). I've been asked to evaluate two final models: one resulting from a backwards selection in R (stepAIC) and one resulting from a backwards selection using PROC LOGISTIC in SAS. The final terms are slightly

How do I go about generating a WEIGHTED mean (and standard error) of a variable (e.g., expenditures) for each level of a categorical variable (e.g., geographic region)? I'm looking for something comparable to PROC MEANS in SAS with both a class and weight statement. Thanks. Marc [[alternative HTML version deleted]]

I'm trying to create a function which will allow me to subset a data set based on values of various specified variables. I also want to then apply some other function(s) (e.g., summary). This is what I've tried so far.... > test.fx <- function(dta, expvar, expval) { + newdta <- subset(dta, eval(expvar)>expval) + summary(newdta$eval(expvar)) + } > >

I am working on plotting a response surface using wireframe(). The default style/orientation is z | | y | \ | \ | \ | \ | \ | \ | \ | \|________________x 0 Now what I want the orientation of axes is: z | | | |

How can I control de "wrap-around" color behaviour in the persp function ? I am using something like : persp(bb[1:100,2:97], col= rainbow(8,start=0.1, end=0.8))) Depending on the rainbow length value I get several "wrap-around" blocks of the selected color range...something that I wanted to avoid... My idea is to use the color in order to make a separation from a certain

I'm attempting to select a model using stepAIC. I want to use a forward selection procedure. I have specified a "scope" option, but must not be understanding how this works. My results indicate that the procedure begins and ends with the "full" model (i.e., all 17 independent variables)...not what I expected. Could someone please point out what I'm not

Hi, Is there a way to plot more than one trend surface using the functions contourplot() and wireframe(). I have found an add=T in contour(), but no equivalent argument in contourplot() and wireframe()? I have taken the example 11-2 (pages 441-451) from Design and analysis of experiments (Montgomery 2001, 5th edition) to see if this could be done in R. I have managed to plot individual

Hallo! ? I just want to make a 3D plot of a surface of a cone and want to plot some single points around. ? I tried wireframe but cannot find how to plot single points ? I tried scatterplot3d but there the surface is not simple to plot. And: How can I rotate the point of view by the z-axis ? I tried persp3d but how can I add some single points? ? Example: ? library(lattice) library(scatterplot3d)

I'm trying to simulate trend data over a five year period. I want different trend profiles...the simplest being a linear trend. I've been using the following code: patBdta1 <- NULL for(i in 1:100) patBdta1 <- rbind(patBdta1,c(yr1= mean(rbinom(50,1,.50)), yr2 =mean(rbinom(50,1,.51)), yr3 =mean(rbinom(50,1,.52)),

Greetings, A couple of questions: 1. I am using wireframe. It prints 3d plots nicely on screen, but when I use a postscript device the z-axis label merges with the z-axis. Is there any option to control their separation? 2. Again, using wireframe, I rotate the whole plot using the screen parameter. However, the x- and y-axis labels (xlab and ylab) remain horizontal. How could I rotate them?

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

It is a trivial matter to adjust the linewidth of the axes and plot line in a 2D plot: x <- 1:10; y <- 1:10; par(lwd=2); # axis (and others) linewidth plot(x,y,"l",lwd=6); # plot linewidth Is there anything comparable to this for persp or wireframe? I've tried a number of ways, and the framebox/axes seem to be linked to the linewidth used in generating the

Hi, trying this again - any help would be great! Basically, I have two matrices of equal dimension, one of them should produce something similar to a heatmap.. The 2nd matrix should be the "heights" for each value on the heatmap - producing a sort of 3d surface plot. Viewing this seems like a problem too, as I need to scale the x axis so that ~140 labels are visible. I've tried

Hi all, I'd like to plot the Real and Imaginary parts of some f(z) as two different surfaces in wireframe (the row/column axes are the real and imag axes). I know I can do it by, roughly speaking, something like plotz <- expand.grid(x={range of Re(z)}, y={range of Im(z), groups=1:2) plotz$func<-c(Re(f(z),Im(f(z)) wireframe(func~x*y,data=plotz,groups=groups) But that seems like a

I'm attempting to use stepAIC to select a model through a forward procedure. I want to consider up to all 3-way interactions. I've attempted to use the following code: m2.Fwd3way <- stepAIC(m1.Ionly, direction="forward", scope=list(upper=~(var1 + var2 + var3 + var4)^3, lower=~1)) When I submit this the trace indicates that only 2-way