similar to: error

Hi, i have such a code tau<-0 for (i in 1:10) { x<-rnorm(20,0,1) # t1=onesampb(x,est=tauloc,SEED=F)$conf.interval if(t1[1]>0 ||t1[2]<0)tau=tau+1 } print (tau) this code gives me Error in if (t1[1] > 0 || t1[2] < 0) tau = tau + 1 : missing value where TRUE/FALSE needed what can be done with such a warning message? i tried x<-x[!is.na(x)] but didnt work. Best

Hi, I wanted to make a transformation with "avas" and "ace" but saw a message " couldn't find function "avas" " what are the possibble reasons of this case? It may be a basic question but unfortunately I am very new in R. Thanks for your helps Regards [[alternative HTML version deleted]]

Hi, Does any one know how we can decide on the correct transformation in (avas and ace) after having drawn the graphs y,g(y) x ,s(x) and g(y) ,s(x) . Is it possible by only looking at patterns the graphs follow for example when y ,g(y) shows a logaritmic pattern can we say that log transform on y is suitable? Thanks for your help. Regards [[alternative HTML version deleted]]

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,

Hello, I am wondering about the X11() graphic device on Windows. I try to plot a raster image but nothing gets displayed. I found some pages where it is mentioned that x11() not always supports raster rendering. Is there any add on for x11, any update or any R-package which solves that displaying problem in Windows? What I try to test it is an example from the package {raster}: library(raster)

Hello, I am wondering about the X11() graphic device on Windows. I try to plot a raster image but nothing gets displayed. I found some pages where it is mentioned that x11() not always supports raster rendering. Is there any add on for x11, any update or any R-package which solves that displaying problem in Windows? What I try to test it is an example from the package {raster}: library(raster)

Hi guys, I try to set up two identical installed up to date CentOS6 machines with Glusterfs/CTDB/Samba . I have set up Glusterfs and it works. I have set up CTDB from CentOS and it seems to work too. Samba is AD integrated and works mainly. The main problem is that file locking seem to not work between the machines at all. If two Win7 clients try to open an document from the same Samba server

Hello list, I'm trying to figure out how exactly the specification of nested random effects works in the lmer function of lme4. To give a concrete example, consider the rat-liver dataset from the R book (rats.txt from: http://www.bio.ic.ac.uk/research/mjcraw/therbook/data/ ). Crawley suggests to analyze this data in the following way: library(lme4) attach(rats) Treatment <-

Hello Everyone,   Learning about R and its wonderful array of functions. If it's not obvious, I usually try to find out what a function stands for. I think this helps me remember better.   One function that has me stumped is "sink." Can anyone tell me if this stands for something?   Thanks,   Paul         __________________________________________________ [[alternative HTML

Ronaldo, Further to my previous posting on your Glycogen nested aov model. Having read Douglas Bates' response and Reflected on his lmer analysis output of your aov nested model example as given.The Glycogen treatment has to be a Fixed Effect.If a 'treatment' isn't a Fixed Effect what is ? If Douglas Bates' lmer model is modified to treat Glycogen Treatment as a purely

The text below is a part of, some work I have to do, which is due in 2 days and I am strung up with a lot of other stuff, so I was hoping someone would take 5 mins and help me ?? Here is a part of my data.frame: year country1 country2 contig comlang pop1 gdp1 pop2 gdp2 rta dist avgflow 1 1992 AUS AUT 0 0 17.4950008 321708.281

Dear all, During my pre-R era I tried (yes, tried) to understand mixed models by working through the 'rat example' in Sokal and Rohlfs Biometry (2000) 3ed p 288-292. The same example was later used by Crawley (2002) in his Statistical Computing p 363-373 and I have seen the same data being used elsewhere in the litterature. Because this example is so thoroughly described, I thought

Hi, I''m trying to understand the lme output and procedure. I''m using the Crawley''s book. I''m try to analyse the rats example take from Sokal and Rohlf (1995). I make a nested analysis using aov following the book. > summary(rats) Glycogen Treatment Rat Liver Min. :125.0 Min. :1 Min. :1.0 Min. :1 1st Qu.:135.8

Your response nicely clarifies a question that I've had for a long time, but which I've dealt with by giving each subject a unique label. Unless I'm missing something, both techniques should work as the toy example below gives exactly the same output in all 3 cases below (forgetting about the convergence problem). Would there be a reason to prefer labeling the levels one way or

Fellow R users, I'm using the BCE {BCE} function to run a Bayesian sediment mixing model. The aim is to find the optimum % contribution from each of the 4 source areas that can yield the target geochemistry. I have geochemistry for 4 source areas called Rat: Rat<-read.table(text="CaO MgO Na2O Al2O3 Topsoils 2.511250 0.7445500 0.7085500 14.10375 ChannelBanks

Hi, hackers : Yesterday night and this morning, I cvsup my system to latest STABLE branch source (RELENG_4). All make buildworld, buildkernel, and installkernel procedure was completed with no error. But when I boot into single user mode and try to make installworld, I always got " touch: not found " error. Does anyone encounter this situation ? I have no idea about that. :-( Here is

Steffen, When the new biomaRt tries to load it errors out because I do not have RMySQL installed. There is not a Windows binary for RMySQL and it does contain C code that I do not know how to build. I do not use the MySQL option in biomaRt. Does RMySQL need to be a required dependency? Below is my screen output and sessionINfo. require(biomaRt) Loading required package: biomaRt Loading required

Hi All, I would like to ask a question on fixed and random effecti in lme. I am fiddlying around Mick Crawley dataset "rats" : http://www.bio.ic.ac.uk/research/mjcraw/statcomp/data/ The advantage is that most work is already done in Crawley's book (page 361 onwards) so I can check what I am doing. I am tryg to reproduce the nested analysis on page 368:

Hi! All, I am working on a dataset 'rat' with dimension 20500x363. I have calculated pca of samples (columns). Now I am trying to plot first two principle components with specified columns in different color. I have done following so far: > dim(rat) [1] 20500 363 >#specifying columns to be colored in red >