similar to: plot logist

Hi, exist in R or any other program for linux that estimate the best fit for data using severals functions? Somethink like tablecurve Thanks Ronaldo -- Um menino-prodigio e uma crianca cujos pais tem muita imaginacao. -- Jean Costeau -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior ] |> V [UFV/DBA-Entomologia

Hi all, In my continuous transition of GLIM to R I try to make a glm with binomial errors. The data file have 3 vectors: h -> the factor that is ajusted (have 3 levels) d -> number of animais alive (the response) n -> total number of animals To test proportion of alive, make d/n. In GLIM: $yvar d$ $error binomial n$ $fit +h$ scale deviance = 25.730 (change = -9.138) at cycle 4

Per http://rwiki.sciviews.org/doku.php?id=developers:projects:gsoc2010 -- and http://rwiki.sciviews.org/doku.php?id=developers:projects:gsoc2010:syrfr -- I am applying to mentor the "Symbolic Regression for R" (syrfr) package for the Google Summer of Code 2010. I propose the following test which an applicant would have to pass in order to qualify for the topic: 1. Describe each of the

Hi, I make an analysis sequence on R. In some cases the function glm.nb fail to ajust the model. Its Ok. The problem is that this error stop the program. I need treat this error and not stop the program. Something like this: ... model <- glm.nb(y~x,maxit=1000) if(glm.nb fail) { teste[i] <- 0 } else { teste[i] <- anova(modelo)$"P(>|Chi|)"[2] } ... I try this: ...

Hello, I have a simple theorical question about regresion... Let's suppose I have this: Model 1: Y = B0 + B1*X1 + B2*X2 + B3*X3 and Model 2: Y = B0 + B2*X2 + B3*X3 I.E. Model1 = lm(Y~X1+X2+X3) Model2 = lm(Y~X2+X3) The Ajusted R-Square for Model1 is 0.9 and the Ajusted R-Square for Model2 is 0.99, among many other significant improvements. And I want to do the anova test to choose the best

Package: xen-utils-3.2-1 Version: 3.2.0-5 Severity: important virtinst Debian package is pointing the pygrub path to /usr/bin/ and hvmloader to /usr/lib/xen/boot/. I suggest ajust the path creating a symbolic link of /usr/lib/xen-3.2-1/bin/pygrub to /usr/bin/pygrub and /usr/lib/xen-3.2-1/boot/hvmloader to /usr/lib/xen/boot/hvmloader. In the last case, need ajust the directory

Hi all, I've got a problem with the nls function. I have an adjustment which works when I fix one of the argument of my function (Xo=150) : *Xo*=150 f<- function (tt*,Xo*,a,b) ifelse(tt<*Xo*,a*exp(-b**Xo*),a*exp(-b*tt)) ajust<-nls(RER~f(tt,*Xo*,a,b),data=data.frame(tt=Ph2[,2*k],RER=Ph2[,2*k+1]),start=list(a=0.5,b=0.014)) But, when I use it as a "normal" parameter (and

Anybody use Sciviews with Linux or MacOS X? It work with these systems? Thanks Ronaldo -- A hacker does for love what others would not do for money. -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior ] |> V [UFV/DBA-Entomologia ] | / \ [36571-000 Vi??osa - MG ] |> /(.''`.)\

Hi, I know that Rsquare in glm or in non-linear models is "wrong", but some people like this. How I make to estimate the Rsquare from a model ajusted with glmmPQL or another GLMM? Thanks for all Ronaldo -- A simplicidade ?? o ??ltimo degrau da sabedoria. -- Victor Hugo -- |> // | \\ [***********************************] | ( ?? ?? ) [Ronaldo Reis J??nior

Hi, If I have the follow example: y = count data x1 = continuous data x2 = qualitative data (2 factors levels) I make a model glm(y~x1,family=poisson) The model has only one curve, to calculate r-square, I do "x1 Deviance"/"total Deviance" But in a model with x2 glm(y~x1+x2,family=poisson) I have two curves, but I dont have, using summary() or anova(), the deviance

Dear all, I'm trying to estimate the parameters of a lognormal distribution fitted from some data. The tricky thing is that my data represent the time at which I recorded certain events. However, in many cases I don't really know when the event happened. I' only know the time at which I recorded it as already happened. Therefore I want to fit the lognormal from the cumulative

Hi, I have submitted a package to CRAN in ftp://ftp.tuwien.ac.at/incoming/pastecs_1.0.tar.gz. The R extension manual advises to send a mail to cran at ci.tuwien.ac.at about it, but it fails. I have also tried WWWadmin at ci.tuwien.ac.at. Are these email adresses still valid? Best, Philippe Grosjean ...........]<(({?<...............<?}))><............................... ) ) ) )

Hi all, I have a bit of a newbie question here that I hope y'all can help with. I've run a cluster analysis using hclust on about 500 objects (using R1.3 under Win 2000). The problem is that the tips of the dendrogram are so close together on the plot that the labels overlap and are unreadable. I've used "cex" to reduce the label sizes but this isn't sufficient with so

I stumbled across this link and was wondering if this has any relation to the R project. http://www.fas.umontreal.ca/BIOL/Casgrain/en/labo/R/index.html -jennifer ----------------------------------- J. Steinbachs, Ph.D. Computational Biologist http://compbiology.org ----------------------------------- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing

Hello all, Since I got no answer to my first mail, I suspect I did not formulate my question adequately, so I try again. In R, one can debug or trace a function 'foo' with debug(foo) or trace(foo), respectively. This leads to a special treatment of the function 'foo' until one enter undebug(foo) or untrace(foo). I would like to know if there is a convenient way to know at any

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Is it possible to create a time series from existing data frame? One column of that frame contains dates (like "01/01/01") and the second one the values. The period is non-regular (it is not strictly daily nor working-days-daily). Thank you. lukas - -- Lukas Kubin lukas.kubin at permonik.com phone: 00420603836180 -----BEGIN PGP

I'm doing a non linear regression with 8 parameters to be fitted: J.Tl.nls<-nls(Gw~(a1/(1+exp(-a2*Tl+a3))+a4)*(b1/(1+exp(b2*Tl-b3))+b4),data=Enveloppe, start=list(a1=0.88957,a2=0.36298,a3=10.59241,a4=0.26308, b1=0.391268,b2=1.041856,b3=0.391268,b4=0.03439)) First, I fitted my curve on my data by guessing the parameters'

Hi, I have a little doubt about contrib packages. I make some little functions in a R package format, this functions are for my personal use (envelopes for glm models and a plot bar routine). Some friends ask me about the liberation of this functions in a web site, R contrib etc. what is the correct procedure for send a package to R contrib? The functions is not create for me, is adaptaded to R

Buenas noches: Les consulto por lo siguiente, tengo este gráfico Estos son el resultado de modelos mixtos, los datos son almacenados en un data.frame, resumiendo el código, quitando el x <- data.frame() como el plot() El primero es fitted: goatMilk$dim, fitted(fitWilm_goat) plot la línea de arriba El segundo es predict: goatMilk$dim, predict(fitWilm_goat) El tercero es plot + lines

Dear R and S users, I've been working on fitting finite mixture of negative exponential distributions using maximum likelihood based on the example given in MASS. So far I had much success in fitting two components. The problem started when I tried to extend the procedure to fit three components. More specifically, likelihood = sum( ln(c1*exp(-x/lambda1)/lambda1 + c2*exp(-x/lambda2)/lambda2