similar to: sampling random groups with all observations in the group

Hi I am an engineer and am running the package mgcv and specifically the command gamm (generalized additive mixed modelling), with random effects. i have a few queries: 1. When I run the command with 1000/2000 observations, it runs ok. However, I would like to see the results as in vis.gam command in the same package, with the 3-d visuals. It appears no such option is available for gamm in the

I have a question to everybody. After log10 transfprmation, I have done ordinary kriging in gstat in R? I need to back trnasform the prediction values to orgiginal scale. How do I do this in gstat in R? Thanks Zia -- Zia Uddin Ahmed 915 Brad Field Hall Department of Crop and Soil Cornell University Ithaca NY 14850 USA

Hi, Apologies if my problem sound silly! I've been modifying the status.xsl , and now stuck with the problem. I want to disply on the top header (h2) total number of sources mounted on the server . I dont know how/where the webserver is working, thus i dont know how/where to get the total count (if its already stored somewhere at all). Otherwise, what is the easiest way to get the count?

Hi, I'm trying to setup tunnel between a Linux VM (running on a KVM) and my Mac 10.8.4. On the Linux box i am able to start tincd and get interface up for my VPN name, but on the Mac, when i run i get this error ?1 tinc/zt ? sudo /opt/local/sbin/tincd -c /opt/local/etc/tinc/zt -D --debug=4 tincd 1.0.21 (Sep 21 2013 00:13:28) starting, debug level 4 Could not open /dev/tun0: No such file or

R helpers I am getting to grips with R but came across a small problem today that I could not fix by myself. I have 3 text files, each with a single column of data. I read them in using: myData1<-scan("C:/Program Files/R/myData1.txt") myData2<-scan("C:/Program Files/R/myData2.txt") myData3<-scan("C:/Program Files/R/myData3.txt") I wanted to produce a

Hi: I want to give an index with all the dates between Sept. to Nov. as 1, and anything else is 0. It doesn't matter which year it is, as long as it is between Sept. to Nov, then set up to 1, otherwise is 0. My data frame looks like below: ID Date 201 1/1/05 6:07 AM 201 3/27/09 9:45 AM 201 9/29/09 8:44 AM 203 10/16/08 10:01 AM 203 10/28/08 9:45 AM 203 10/31/08 11:12 AM 203 11/7/08 11:32 AM

Hi, Does anyone here have experience running Tinc on NixOS? I'm trying to run Tinc on a NixOS machine, using the similar configuration i had for Ubuntu. My home subnet is 192.168.1.0/24 and my work is 10.16.0.0/24. However, unlike ubuntu, when I start tincd on nixos, and try to 'ifconfig $INTERFACE 192.168.1.10 netmask 255.255.255.0' in my tinc-up, I loose network access on the box

Hi, I'm trying to setup tinc tunnel between my work and home linux machines as shown here http://cl.ly/image/1H1R2X2D403X I am able to setup tinc session between A and B and I can ping them. But I can't ping C or D from B, or A respectively. Home network 192.168.1.0 Work network 10.20.50.0 tinc tunnel is 15.0.0.x. Do I need to run tinc on C and D as well? Or Do i need to manually

Hi Bill, Thanks for prompt reply and letting me know a way around it. I have more than 1200 observations and not all the values are the same. However, my data points are quite similar, for example, 0.079275, 0.078867, 0.070716 in millions and etc. I have run the data without converting it to millions and I still get the same error message. As I have kurtosis value, it should be fine for the

Dear all, I am pretty new to R only having an introduction course, so please bare with me. I am doing my PhD at The Max Planck Institute of Immunobiology where I am analyzing some calorimetry data from some mice. I have a spreadsheet consisting of measurements of the respiratory exchange rate at different time points measured every 9 minutes over some days. My goal is "bin"/average the

moments::anscombe.test(x) does give errors when x has too few values or if all the values in x are the same > moments::anscombe.test(c(1,2,6)) Error in if (pval > 1) pval <- 2 - pval : missing value where TRUE/FALSE needed > moments::anscombe.test(c(2,2,2,2,2,2,2,2)) Error in if (pval > 1) pval <- 2 - pval : missing value where TRUE/FALSE needed You can use tryCatch() to

Another bad case is > moments::anscombe.test(rep(c(1,1.1),length=35)) Error in if (pval > 1) pval <- 2 - pval : missing value where TRUE/FALSE needed I haven't checked the formulas carefully, but I suspect the problem is from taking the cube root of a negative number in z <- (1 - 2/(9 * a) - ((1 - 2/a)/(1 + xx * sqrt(2/(a - 4))))^(1/3))/sqrt(2/(9 * a)) In R, the

Hi, I have a dataset including monthly date from 1971-01-01 to 2009-01-01. The dates are character variables. I want to select the subsample from my original data, so I use mydata3<-window(mydata,start="1972-02-01",end="2005-02-01"), or mydata3<-window(mydata, start=as.date("1972-02-01"),end=as.date("2005-02-01")) or mydata3<-window(mydata,

Hi all, While running the anscombe.test in R, I'm getting an error of *Error in if (pval > 1) pval <- 2 - pval : missing value where TRUE/FALSE needed* for a few time series columns whereas for most of the series the function is working fine. I have checked for those specific columns for missing values. However, there is no NA/NAN value in the dataset. I have also run kurtosis for

Dear users, I have a double for loop that does exactly what I want, but is quite slow. It is not so much with this simplified example, but IRL it is slow. Can anyone help me improve it? The data and code for foo_reg() are available at the end of the email; I preferred going directly into the problematic part. Here is the code (I tried to simplify it but I cannot do it too much or else it

As a followup to pi-day, I attempted to make a .gif of a simulation based estimation of pi by plotting points inside a single quadrant of a circle (a la?http://www.drewconway.com/zia/?p=2667 ). ?When rendering the individual x,y pairs with points() I intermittently see points crop up around (2,0.5) but the input values for x and y are bounded between 0 and 1. square<-structure(c(0, 0, 1, 1, 0,

I'm experience echo on outgoing calls: Snom 200 ----> Asterisk ----> T100P ----> PRI ----> called party I am getting echo on the Snom 200 phone. The called party does not hear the echo. Since the PRI is digital, I don't really understand where the echo is coming from. I turned on echotraining, echocancellation=yes (128), and echowhenbridged in zapata.conf, to no avail.

Hi, I have a problem. I have a data frame looking like: ID val A? .3 B? 1.2 C? 3.4 D? 2.2 E? 2.0 I need to CREATE the following TABLE: CASE?? DIFF A-A??? 0 A-B??? -0.9 A-C??? -3.1 A-D??? -1.9 A-E??? -1.7 B-A??? ... B-B??? ... B-C B-D B-E C-A C-B C-C C-D C-E D-A D-B D-C D-D D-E E-A E-B E-C E-D E-E WHERE CASE IS THE COUPLE OF ELEMENTS CONSIDEREDM AND DIFF IS THE computed DIFFERENCE between

I downloaded the iid.test, but I can't run it. I get the following message: Error: could not find function "iid.test" Where am I supposed to save this package in order that it works? Thanks, EZ [[alternative HTML version deleted]]

This is my first post to R-help and I look forward receiving some advice for a novice like me... I?ve got a simple repeated (4 periods so far) 10-question survey data that is very easy to work on Excel. However, I?d like to move the compilation to R but I?m having some trouble operating on count list data in a neat way. The data C > str(C) 'data.frame': 551 obs. of 13