similar to: ts.data plot

Dear List, i am trying to produce a 3d plot using wireframe using the code: wireframe(Residuals_FD ~ Elevation * Temperature, data = data2, scales = list(arrows = FALSE), drape = TRUE, colorkey = TRUE) As you can see when the code (using the data below) is run the plot area is set-up correctly but the actual surface is missing? Any help would be greatly appreciated. Chris #data Elevation

would someone be kind enough to paste the code below into an R session ( ir you can paste it into a file and just source it ) and take a look at it ? I must be doing something wrong but i can't find it. I start out with a zoo object that has 100 elements in it. then, i only want to keep the rows in which the bid and the ask are both grester than zero. ( maybe my logic is wrong there ? ).

R Help - I'm trying to run a correlation matrix with a covariate of "age" and will at some point will also want to covary other variables concurrently. I'm using the "psych" package and have tried other methods such as writing a loop to extract semi-partial correlations, but it does not seem to be working. How can I accomplish this? library(psych) > set.cor(y =

Hi R Gurus! I conducted regression and anova followings : Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 6.07e-01 5.95e-02 10.19 < 2e-16 *** nemp 2.87e-06 1.04e-07 27.63 < 2e-16 *** as.factor(corridor1)A -8.81e-02 2.13e-02 -4.14 3.6e-05 *** as.factor(corridor1)B

Here I am again with question I'll feel foolish for asking, when I see the answer. I'm trying to produce a report and here's where I get stuck: How do I get R2HTML to produce the same number format? Particularly remove the decimal places for Par and Sal. Are there better methods to produce this type of report? Thanks, L.A. R version 2.10.0 XP

This should be easy but it does not work I have 3 vectors*(activeT,inactT, activeR)*, the idea is that if the last value in inactT is higher than the last in activeT this value has to be append in active T and the last value in another vector call activeR has to be repeated. (at the bottom you can find the vectors) I have done this: activeT=round(as.numeric(activeT)) inactT=

Hi, I'm proposing to make vectorizer-maximize-bandwidth on by default for loop vectorizer because it should generally help performance. I've tested the performance impact on Intel sandybridge machine with speccpu benchmarks: Benchmark Base:Reference (1) ------------------------------------------------------- spec/2006/fp/C++/444.namd 26.84

To anyone who can help: I have two brief questions concerning sapply. Following below is the code for my example. The two problems are described at the end of the code: site <- rep(2:6, each = 12) tillage <- rep(c(1,-1), each = 6, times = 5) carbon <- c(18.23, 16.06, 17.81, 16.07, 17.26, 17.08, 14.92, 15.88, 12.11, 14.23, 16.99, 13.57, 20.34, 20.3,

The pretty picture that I saw at: http://chartsgraphs.wordpress.com/2009/02/09/r-panel-chart-beats-excel-chart/#more-1096 inspired me to try something similar. The code that I wrote is: ------snipsnip--------------------------------------------------------------------- M <- structure(list(date = structure(c(13634, 13665, 13695, 13726, 13757, 13787, 13818, 13848, 13879, 13910, 13939, 13970,

When using get.hist.quote(), I find the dates are broken. This is with R 2.1.1 on Mac OS X `panther'. > library(tseries) Loading required package: quadprog 'tseries' version: 0.9-27 'tseries' is a package for time series analysis and computational finance. See 'library(help="tseries")' for details. > x <-

Howdy I apologize for duplicated posting. But I decided to correct my previous posting. I had the regression results using r <- lm(Y ~ nemp + as.factor(devt), data=d). First, there is the result of anova(r). Here I could not find regression coefficients. Response: Y Df Sum Sq Mean Sq F value Pr(>F) nemp 1 58.2 58.2 1233.23 < 2e-16 ***

The aggregate function for some reason will now work for me. The error I'm getting is: "Error in sort.list(y) : 'x' must be atomic for 'sort.list' Have you called 'sort' on a list?" agPriceList=aggregate(PriceList$Size, list(PriceList$bandNum),sum) *Price list dataframe:* dput(PriceList) structure(list(Price = c(0, 8.18, 8.27, 10.42, 10.5, 10.6, 11.13,

Hi..   R Program is shown ARIMA output as below then SARIMA equation is be   (1 - 0.991B^{12})z_t + 43.557 = (1+0.37B)(1-0,915B^{12})a_t    But I try to calculate it by manual . It look like it 's big different from R sofeware,   I am not sure this equation is correct or not . PLS supoort me to confirm it   Arima Model ( 0,0,1)(1,0,1)   No Transformation Constant   >> 43.557 , t = 10.09

Hello, I'm a bit new to the world of R so forgive my ignorance. I'm trying to do a zero-inflated negative binomial regression and have received an error message and i'm not sure what it means. I'm running R 2.5.1 on XP. I have just tried a really simple version of the model to see if it would run before I put all the variables in. I have attached all the variables to the

Dear Rusers, I have met a difficult problem on explaining the differences of principal component analysis(PCA) between R,S-PLUS and SAS/STATA/SPSS, which wasn't met before. Althought they have got the same eigenvalues, their coeffiecients were different. First, I list my results from R,S-PLUS and SAS/STATA/SPSS, and then show the original dataset, hoping sb. to try and explain it.

Hi, I am running a Cox Mixed Effects Hazard model using the library coxme. I am trying to model time to onset (age_sym1) of thought problems (e.g. hearing voices) (sym1). As I have siblings in my dataset, I have decided to account for this by including a random effect for family (famid). My covariate of interest is Mother's diagnosis where a 0 is bipolar, 1 is control, and 2 is major

Hi, I have a problem using the function 'reshape' for the following dataset: input: X Y Z T.U.1 T.U.2 T.V.1 T.V.2 1.00 1.00 A 24.94 24.74 34.65 35.33 1.00 3.00 A 26.03 25.87 25.76 26.12 1.00 1.00 B 24.76 24.63 34.34 34.48 1.00 3.00 B 26.15 25.98 25.96 25.66 2.00 1.00 A 24.40 24.26 34.55

Is there a way to take the effective user output from top and figure out which AD Domain User is indicated by the smbd process? 2834 root 20 0 584928 38292 5784 R 24.6 0.2 26:06.20 samba 6448 3000042 20 0 643128 44424 12368 S 20.3 0.3 15:23.92 smbd 28432 3000047 20 0 635816 38856 6864 S 1.3 0.2 0:00.04 smbd 28433 3000024 20 0 635816 38868 6876 S

hello, i'd appreciate help with my glmer. i have a dependent which is an index (MH.index) ranging from 0-1. this index can also be considered as a propability. as i have a fixed factor (stage) and a nested random factor (site) i tried to model with glmer. i read that it's possible to use a quasibinomial distribution, for this kind of data, which i than actually did - but firstly (1)

Hello there, I am not quite sure how to interpret the output of Rprof (in the following the output I was staring at). I was poking around the web a little bit for documentation but without much success. I guess if I want to figure out what takes so long in my code the 2nd table $by.total and the total.pct column (pct = percent) is the most helpful. What does it mean that [ or [.data.frame is