search for: formulate

Dear All: I need to have some data frame objects. First aa object: pH Formulation time Subject [1]1.2 F 0 1 [2]7.4 S 1 2 [3] MF 2 3 [4] 3 4 [5] n i Then, I need to produce 2*3(pH*formulation) different tables. This table includes column of (pH, Formulation, time S1 S2 S3 ?K..Si) and S1=

...erstand the use of the constraints in matrix form. I use them like: constr_mat<- -diag(2) constr_vec<- rep(-0.05,2) constr_mat<- rbind(constr_mat, diag(2)) constr_vec<- c(constr_vec, rep(-1, 2)) To get parameters in the interval [-1, 0.05]. (And this works so far) Now I would like to formulate a constraint like Sum(Abs(x_i))<=limit But I have no idea how I could get such a condition by using the matrix constraint formulation. I have already searched the help, web and all searchable R resources I could find, but got no helping result. Is constrOptim the right function to use? Mayb...

Dear all: I am a new R-user and I have 2 questions about it. 1) I have a dataframe. Based on ?formulation? and ?subject?, a dataframe is split into 4 dataframes. The example is as follows. Moreover, I want to calculate ?test? value for these 4 dataframes. My question is that the ?test? values not correct and I do not know where the problem is. 2) There are 12 ?test? (y) values from

HI, sorry but i don't understand how to make a function with as.function() formula<-"2+3*x" formu<-as.symbol(formula) > formu 2+3*x formul<-as.function(alist(x=,formu)) curve(formul,1,5,col="blue") Error in xy.coords(x, y, xlabel, ylabel, log) : x and y lengths differ > typeof(formul) [1] "closure" and not plot the curve function, Why?

...action effects at site b correlation between site a and b study*treatment interaction effects residual (between-subject) effects at site a residual (between-subject) effects at site b correlation between site a and b residuals (between-subject) effects I'm having trouble seeing how to formulate this correctly in lme() Hope someone can help ... cheers Bob Kinley > Toy Study Treatment Subject Site Result 1 T1 1 a 13.901820 1 T1 1 b 19.158889 1 T1 2 a 16.026299 1 T1 2 b...

I'm trying to fit a nested mixed model using lmer and have some questions about the output and my model formulations. I have replicate measures on Lines which are strictly nested within Populations. (a) So if I want to fit a model where Line is a random effect and Populations are fixed and the random Line effect is constant across Populations, I have: measure_ijk = mu + P_i + L_ij +

...00 to 1700 nm. Instead of fitting with octane[i] = a[0] * NIR[0,i] + a[1] * NIR[1,i] + ... I want to fit the data with: octane[i] = a[0] * NIR[0,i] + a[1] * NIR[1,i] + ... + b[0]*NIR[0,i]*NIR[0,i] + b[1] * NIR[0,i]*NIR[1,i] + ... i.e. quadratic with interaction terms. But I don't know how to formulate this. May I have some help please? Thanks, Kelvin [[alternative HTML version deleted]]

I need a little piece of advice concerning passing data frames into the functions. As I do a lot of similar fits at a time, I'd like to write a small function doing the fits for all relevant variables automatically. However, I usually get error messages of the following kind: (I present here a part of a test code). ##################################################### # Data set:

Hi all - R2.0.1 on OSX;MASS library;nlme library I am trying to emulate the solution to a problem set that has normally been run in Genstat, using R. The problem that I am having at the moment is with the following glmm question (using glmmPQL from the MASS library): "We have two different forest habitats (first rotation thicket, and high forest) which we want to survey for the

Hi, I would like to see if a matching variable is an effect-modifier in a conditional logistic regression. Naturally, the matching variable can't enter directly in the model but as an interaction with terms that are in. However, I have problems in formulating the correct model the term that's already in the model is a factor. I am using treatment contrasts and the problem is that if I

Hello, I'm trying to analyse data from an incomplete design with four factor : - fr : number of the batch - op : ID of operator - meth : method used - mat : nature of the material used and one variable - mv : mesure - trmv : transformed mesure str(matvol) `data.frame': 120 obs. of 6 variables: $ fr : Factor w/ 30 levels "1","2","3","4",..:

...he problem is a simple balanced parentheses problem in CFL-reachability, and it can be computed efficiently. The paper you mentioned is a very nice paper that proposed an efficient algorithm to solve Dyck-CFL-reachability problem in bidirectional graphs. However, precise pointer analysis cannot be formulated as a Dyck-CFL-reachability problem. The alias analysis presented in that paper is an over approximation. Hence the time bound from our paper still applies. Regards, Lian On Fri, Oct 18, 2013 at 4:39 PM, Daniel Berlin <dberlin at dberlin.org> wrote: > I notice you guys formulate your...

...exactly fits into the bidirectional graph problem, unless i'm missing something. > > The paper you mentioned is a very nice paper that proposed an > efficient algorithm to solve Dyck-CFL-reachability problem in > bidirectional graphs. Right, and your matching call sites can be formulated as one. However, precise pointer analysis cannot be > formulated as a Dyck-CFL-reachability problem. FWIW: This remains to be seen, and certainly also depends on your definitions of "precise". I agree you cannot currently formulate precise andersen style field sensitive pointer an...

Dear list, we are a statistical/epidemiological departement that - after a few years of rapid growth - finally is getting around to formulate a general data storage and retention policy - mainly to ensure that we can reproduce results from published papers/theses easier in the future, but also with the hope that we get more synergy between related projects. We have formulated what we feel is a reasonable draft, requiring basic...

Hi, I'm trying to make a regression of the form : formula <- y ~ Asym_inf + Asym_sup * ( (1 / (1 + (n1 * (exp( (tmid1-x) / scal1) )^(1/n1) ) ) ) - (1 / (1 + (n2 * (exp( (tmid2-x) / scal2) )^(1/n2) ) ) ) ) which is a sum of the generalized logistic model proposed by richards. with data such as these: x <- c(88,113,128,143,157,172,184,198,210,226,240,249,263,284,302,340) y <-

Dear ?, I'm sure that there are many ways to do what you want; here's one: > cbind(concept_df, category= + ifelse(apply( + sapply(chemical_df$chemical, + function(x) grepl(x, concept_df$concept)), + 1, any), + "chemical", "")) concept category 1 butan

In the almost current development version (2009-05-22 r48594) and also in https://svn.r-project.org/R/trunk/src/library/base/man/factor.Rd ?factor contains (compare the formulations marked by ^^^^^^) \section{Warning}{ The interpretation of a factor depends on both the codes and the \code{"levels"} attribute. Be careful only to compare factors with the same set of levels (in

On 06/10/2010 08:16 PM, Chris Pearce wrote: > I have no plans to change the 'index' packet format further. I have proposed an alternative formulation of the index packet at http://github.com/bemasc/OggIndex/blob/master/Proposed-modified-spec.txt That repository also contains a working implementation of the alternative formulation. I have reviewed the details with Chris extensively.

...ane[i] = a[0] * NIR[0,i] + a[1] * NIR[1,i] + ... >> I want to fit the data with: >> octane[i] = a[0] * NIR[0,i] + a[1] * NIR[1,i] + ... + >> b[0]*NIR[0,i]*NIR[0,i] + b[1] * NIR[0,i]*NIR[1,i] + ... >> i.e. quadratic with interaction terms. >> But I don't know how to formulate this. >> May I have some help please? >> Thanks, >> Kelvin > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/li...

...of fitting with octane[i] = a[0] * NIR[0,i] + a[1] * NIR[1,i] + ... > I want to fit the data with: > octane[i] = a[0] * NIR[0,i] + a[1] * NIR[1,i] + ... + > b[0]*NIR[0,i]*NIR[0,i] + b[1] * NIR[0,i]*NIR[1,i] + ... > i.e. quadratic with interaction terms. > But I don't know how to formulate this. > May I have some help please? > Thanks, > Kelvin [[alternative HTML version deleted]]