Re: [R] nlme questions

and Bates (2000) Mixed-Effects Models in S and S-PLUS (Springer), I suggest you get it; it is excellent for things like this. I'm sorry I couldn't help more. spencer graves Christian Mora wrote: Dear R users; Ive got two questions concerning nlme library 3.1-65 (running on R 2.2.0

Re: [R] nlme questions

and trt2 are defined as numeric (ie 0 1) then it will fit as a numeric. --Matt -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of Christian Mora Sent: Friday, November 18, 2005 4:01 AM To: Spencer Graves Cc: r-help@stat.math.ethz.ch Subject: Re: [R] nlme

[R] nlme questions

that, but I would like to have a different opinion on these two issues. Thanks in advance Christian Mora __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting

[R] error in gnls

convIter not found. This error occurs only with R 2.1.1. Any ideas? Thanks Christian Mora __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] Extract Coeff, Std Error, etc from gnls output

Dear list members; Is there any trick to extract the coefficients along with std errors, t-values and p-values for each beta from a gnls fit model (similar to the results obtained using summary(lm)$coeff for linear models)? Thanks for any hint cm __

[R] custom loss function + nonlinear models

Hi all; I'm trying to fit a reparameterization of the assymptotic regression model as that shown in Ratkowsky (1990) page 96. Y~y1+(((y2-y1)*(1-((y2-y3)/(y3-y1))^(2*(X-x1)/(x2-x1/(1-((y2-y3)/(y3-y1))^2)) where y1,y2,y3 are expected-values for X=x1, X=x2, and X=average(x1,x2), respectively.

RE: [R] covariate selection?

Hi Ian Have you tried help.search(pca)? Christian -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Martínez Ovando Juan Carlos Sent: Tuesday, October 12, 2004 7:56 PM To: Ian Fiske Cc: [EMAIL PROTECTED] Subject: RE: [R] covariate selection? Hello Ian,

[R] intersect two files

Hi all; Im working with two datasets in R, say data1 and data2. Both datasets are composed of several rows and columns (dataframe) and some of the rows are identical in both datasets. Im wondering if there is any way to remove from one set, say data1, the rows that are identical in the other set,

RE: [R] Tutorial for graphics

A fairly good tutorial is “Gráficos Estadísticos con R” by Juan Carlos Correa and Nelfi González (PDF [2110kB]) that you can find in Documentation/Contributed (in Spanish though). -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Phil Sent: Wednesday, June

RE: [R] Non-Linear Regression Problem

Working on the same idea, Ive generated a data grid with 4 vars, two of them with its own sequence and two with fixed values. As Spencer pointed out one option is to get the values from a simple loop. My question is: How can jump from one set of starting values to the next (on the data grid) in

[R] mode

How can I get the mode (most frequent value) from a dataset (continuos variables)? I can obtain it when the data is discrete (by making a table and looking at the higher frequency) but I don't know how obtain it from, for example, a density plot of the data. Does anyone know how to do it? Thanks

[R] Johnson's Sb distribution

Hi all; I'm working with the library SuppDists trying to fit a Johnson's Sb distribution to a dataset. It works fine, but I need to set one of the location parameters (epsilon) to zero. How can I do this using the function JohnsonFit() or any other similar? ...and Is it possible to define the

RE: [R] Sampling

YourData[sample(209,106),] Regards Christian Mora -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Savano Sent: Monday, December 01, 2003 9:38 AM To: Lista R Subject: [R] Sampling UseRs, I imported a table using read.table. It has 209 observations, I

[R] gnls - Step halving....

Hi all, I'm working with a dataset from 10 treatments, each treatment with 30 subjects, each subject measured 5 times. The plot of the dataset suggests that a 3-parameter logistic could be a reasonable function to describe the data. When I try to fit the model using gnls I got the message 'Step

[R] gnls - Step halving....

Hi all, I'm working with a dataset from 10 treatments, each treatment with 30 subjects, each subject measured 5 times. The plot of the dataset suggests that a 3-parameter logistic could be a reasonable function to describe the data. When I try to fit the model using gnls I got the message 'Step

[R] Question about 'NA'

Hi all, Ive got a database with 10 columns (different variables) for 100 subjects, each column with different # of NA's. I'd like to know if it is possible to use a function to exclude the NA's using only a specific column, lets say: Data2 - omit.exclude(Data1$column1) ??, then Data3 -

[R] Lattice graphs and legend

Dear all, I'm working with grouped data using the nlme package in R. By using the function plot() I'm trying to obtain a trellis-like graph but by default a legend is placed on top of the plot. How can I remove this legend? Thanks for any hint Christian Mora

[R] 4-parameter logistic model

Dear R users I'm a new user of R and I have a basic question about the 4-parameter logistic model. According to the information from Pinheiro Bates the model is: y(x)=theta1+(theta2-theta1)/(1+exp((theta3-x)/theta4)) == y(x)=A+(B-A)/(1+exp((xmid-input)/scal)) from the graph in page 518 of the