similar to: saturation binding in nlme

Dear list, Has anyone created a version of SSmicmen that allows testing for group differences? The basic Michaelis-Menten equation is: (Bmax * X) / (Kd + X). The nlme package allows modeling of random effects for Bmax and Kd as needed, but I curious how I can build in group differences? I have receptor binding data for strains of mice, and following Pinheiro and Bates' lead in their

?sink On May 18, 2018 9:47:25 AM PDT, Ed Siefker <ebs15242 at gmail.com> wrote: >I have dose response data analyzed with the package 'drc'. >'summary(mymodel)' prints my kinetic parameters. I want >that text in an ASCII text file. I want to get exactly what I >would get if I copied and pasted from the terminal window. > >I've read the documentation on

I have dose response data analyzed with the package 'drc'. 'summary(mymodel)' prints my kinetic parameters. I want that text in an ASCII text file. I want to get exactly what I would get if I copied and pasted from the terminal window. I've read the documentation on data export to text files here:

Dear everyone, I am doing some work about species richness estimation. Nonparametric estimation (such as Chao1, Jacknife1) can be done just using function "specpool()" and "estimateR()" in package "vegan". The problem is that I can not found any functions for parametric estimation (such as MMMeans, MMruns, Michaelis-Menten). Do you know any function for doing this?

Dear R users: I encountered difficulties in michaelis-menten equation. I found that when I use right model definiens, I got wrong Km vlaue, and I got right Km value when i use wrong model definiens. The value of Vd and Vmax are correct in these two models. #-----right model definiens-------- PKindex<-data.frame(time=c(0,1,2,4,6,8,10,12,16,20,24),

Hello everyone, I have biological data from a competition experiment where a free ligand is titrated against the binding of a protein. Now, I would like to fit a standard on-site binding curve to this data in order to obtain the IC50 and Kd values. Unfortunately I have not been able to find a package/function which allows such a fitting and calculates the results. Does anyone know if there is

this may be a better question for r-devel, but ... Is there a particular reason (and if so, what is it) that the inverse link is not in the list of allowable link functions for the binomial family? I initially thought this might have something to do with the properties of canonical vs non-canonical link functions, but since other link functions (probit, cloglog, cauchit, log) are allowed, I

I'm attempting to calculate a regression in R that I normally use Prism for, because the formula isn't pretty by any means. Prism presents the formula (which is in the Prism equation library as Heterologous competition with depletion, if anyone is curious) in these segments: KdCPM = KdnM*SpAct*Vol*1000 R=NS+1 S=(1+10^(X-LogKi))*KdCPM+Hot a=-1*R b=R*S+NS*Hot+BMax c = -1*Hot*(S*MS+BMax) Y

I'm just getting to grips with using ode function and have used the examples and vignettes to produce a small model of a one-pool, michaelis-menten, enzyme kinetic reaction. The rate of flux of substrate into pool A is constant (fluxoa) however the rate of flux out of pool A is controlled by the HMM equation (v = Vmax/ ( 1 + (Km / Concentration A )) ). This function works fine and

Hello, I have two tab-delimited files which I would like to combine. In the first one I have gene IDs (Unique) on column 1 and than various experimental results from microarray analysis (see attached files list1 ) the second arrays have the same genes IDs (more and in a different order, some are double) (see attached files list2 ) What I would like to do is to search in the second list for gene

Dear R users I recently migrated from Statistica/SigmaPlot (Windows) to R (Linux), so please excuse if this may sound 'basic'. When running a nonlinear regression (V = Vmax * conc / (Ks + conc), i.e. Michaelis-Menten) on SigmaPlot, I get the output listed below: >>>Begin SigmaPlot Output<<< R = 0.94860969 Rsqr = 0.89986035 Adj Rsqr = 0.89458984 Standard Error of

Hi there, I have a data frame as listed below: > Ca.P.Biomass.A P Biomass 1 334.5567 0.2870000 2 737.5400 0.5713333 3 894.5300 0.6393333 4 782.3800 0.5836667 5 857.5900 0.6003333 6 829.2700 0.5883333 I have fit the data using logistic, Michaelis?Menten, and linear model, they all give significance. > fm1 <- nls(Biomass~SSlogis(P, phi1, phi2, phi3), data=Ca.P.Biomass.A)

Hello! I work with : R : Copyright 2006, The R Foundation for Statistical Computing Version 2.3.1 (2006-06-01) On Windows XP Professional (Version 2002) SP2. At this moment I use the function "nls" combined with a selfStar model (SSmicmen, related to Michaelis-Menten equation, and provided by the "stats" package). When I realise the following operation (cf. p 59 of the

Hi, All: The following identifies apparent inconsistencies and proposed fixes between the documentation and the behavior of 'nls'. Specifically, the help file says 'data' 'Can also be a list'. However, I got an error when I tried it with a list that could not be coerced to a data.frame. This is documented in the first example below. Also, the

Dear R experts! I'd to fit data by 'nls' with me-supplied function 'fcn'. 1) I'd like 'fcn' to accept arbitrary arguments, i.e. I defined it as f(...) {<body>}. (Ok, that's not actually impotant). 2) Second, I would NOT like to supply every parameter in the formula. To illustrate this, let's look at the last example of 'nls' help

> I've never looked at the MachO code as I do not have such a platform nor do > I know the file format. > > Could we concentrate on the ELF backend, please. I don't mind using the ELF backend as our test case, it just seems that the ELFWriter/ELFCodeEmitter don't even use the BufferBegin/BufferEnd/CurBufferPtr system exposed by the base MachineCodeEmitter. There is a big

I have tried multiple different methods to figure out how to get a date axis of my preference (start date of each month). Any assistance would be appreciated. The below section is not producing a date axis: plot(totaldays$totaldays,totaldays$y,type="n",ylim=c(0,Emax),xaxt="n") *xlabels<-(strptime(totaldays$totaldays,"%Y-%m-%d",tz=""))

Dear R-users, I have some trouble to perform a 2-samples KS test. Apparently my 2 samples are numerical (see below) but R complains that the "y" term is not.... ------------------------------------------------------------------------------------- > str(subset(mydata,Identified=="NO",select=KD)) 'data.frame': 2889 obs. of 1 variable: $ KD: num -0.272 -0.080

Fit NLE - was: [R] computer algebra in R Original post: https://stat.ethz.ch/pipermail/r-help/2023-November/478619.html Dear Kornad, I think I have started to understand what you try to achieve. The problem is to fit a NLE and compute the parameters of the NL-Eq. I have included the R Help-list back in the loop, as I am not an expert in optimization. Goal: y ~ I0 + IHD * hd + ID * d; where: y

Hello All, Thank you for taking my question. I am looking for information on how R handles interaction terms in a multiple regression using the ?lm? command. I originally noticed something was unusual when my R output did not match the output from JMP for an identical test run previously. Both programs give identical results for the main test and if the models do not contain the interaction