similar to: Two-way Random Effects with unbalanced data

Hi all, How does R know to regard a variable as a factor and not a character? For example, consider the following table: Observation Gender Dosage Alertness 1 m a 8 2 m a 12 3

Hi, I wonder if there is a mac equivalent to ctrl+r shortcut for running a script line.. (I am using MAC OS X 10.7.2). I saw this post http://r.789695.n4.nabble.com/Shortcut-to-run-script-tp1016619p1016619.html but couldn't find an answer... Thanks a lot, Asaf -- View this message in context: http://r.789695.n4.nabble.com/Mac-shortcut-for-running-a-script-tp4190771p4190771.html Sent from

Hi, I'm using the uniroot function, and would like to detect an error which occurs, for instance, when the values at endpoints are not of opposite signs. For example: uniroot( function(x) x^2+1, lower=1, upper=2 ). I want to say something like: if "error in uniroot(...)" return NA else return uniroot$root Thanks a lot! Asaf -- View this message in context:

Hi there, I am trying to use predict() with an object returned by cv.glmnet(), and get the following error: no applicable method for 'predict' applied to an object of class "cv.glmnet" What's wrong? my code: x=matrix(rnorm(100*20),100,20) y=rnorm(100) cv.fit=cv.glmnet(x,y) predict(cv.fit,newx=x[1:5,]) coef(cv.fit) Thanks so much, Asaf -- View this message in context:

Hi Joe & Michael, In rev. 151382 you have changed the fptoui implementation of the x86 codegen for win32. Before the change fptoui was lowered to flds 16(%esp) fisttpll 8(%esp) movl 8(%esp), %eax After the change fptoui is lowered to flds 40(%esp) calll _ftol2 Please note that the assumption that _ftol2 doesn't modify ECX isn't true on sandybridge platform.

Hello. I'm interested in creating a Hebrew stemmer to use with Xapian. Hebrew is a complicated language to stem, as it uses the semitic "root" system, rather than prefixes and suffixes, and has many irregularities in accidence (morphology). Fortunately, two bright fellows from the Technion University in Israel have already created a Hebrew morphological analyzer as part of their

Dear All, I wonder if there is an efficient way to fit the generalized linear mixed model for multivariate outcomes. More specifically, Suppose that for a given subject i and at a given time j we observe a multivariate outcome Yij = (Y_ij1, Y_ij2, ..., Y_ijK). where Y_ijk is a binomial(n_ijk, p_ijk). One way to jointly model the data is to use the following specification: g(p_ijk) =

[ cc: llvmdev -- please email the list rather than me directly ] 2009/7/6 asaf <asaf at inbox.com> > LLVM seems a good vm, why is there no 68k version ? > LLVM is an open-source project with development driven by contributors (companies and individuals), and if there isn't a particular feature implemented, it's because no one has had the necessity or desire to implement it.

I'm using the lme function in nmle to estimate the variance components of a fully nested two-level model: Y_ijk = mu + a_i + b_j(i) + e_k(j(i)) lme computes estimates of the variances for a, b, and e, call them v_a, v_b, and v_e, and I can use the intervals function to get confidence intervals. My understanding is that these intervals are probably not that robust plus I need intervals on the

Dear R user, I am trying to reproduce the results in Montgomery D.C (2001, chap 13, example 13-1). Briefly, there are three suppliers, four batches nested within suppliers and three determinations of purity (response variable) on each batch. It is a two stage nested design, where suppliers are fixed and batches are random. y_ijk=mu+tau_i+beta_j(nested in tau_i)+epsilon_ijk Here are the

Dear R user, I am trying to reproduce the results in Montgomery D.C (2001, chap 13, example 13-1). Briefly, there are three suppliers, four batches nested within suppliers and three determinations of purity (response variable) on each batch. It is a two stage nested design, where suppliers are fixed and batches are random. y_ijk=mu+tau_i+beta_j(nested in tau_i)+epsilon_ijk Here are the

Dear R list, What is the proper way to specify a nested model so that the F values agree with the expected mean square errors? Specifically, suppose I have a design where "Heads" are nested within "Machines". I would like to model the following Y_ijk = Mu + Machine_i +Head_j(i) +Error_k(ij). Using the commands below, > summary(aov(Strain~Machine + Head%in%Machine ))

Hi al, I have a dataset (see attached), which basically involves 4 treatments for a chemotherapy drug. Samples were taken from 2 biopsy locations, and biopsy were taken at 2 time points. So each subject has 4 data points (from 2 biopsy locations and 2 time points). The objective is to study treatment difference.? I used lme to fit a mixed model that uses "biopsy.site nested within pid"

hi there im trying to to raise the smbd deamon but i can't!!!! and in the log i get an error... User tftp in passdb, but getpwnam() fails! can ypu please give me an answer assaf

Hi, I'm trying to achieve the following feature that's common in Avaya systems: A user page all extensions in a full duplex mode- he can hear all, and all can hear him via their phones' speaker. When one of the extensions picks up the handset, the call is bridged between the pager and the person who picked up. All the rest are disconnected. Does anyone have an idea or ever