Re: [R] problem in font

Em Segunda 18 Setembro 2006 10:02, Uwe Ligges escreveu: > Ronaldo Reis-Jr. wrote: > > Hi, > > > > But, I have this warning in the enviromnent that have this problem: > > > > During startup - Warning message: > > using .GlobalEnv instead of 'package:l

[R] problem in font

Hi, But, I have this warning in the enviromnent that have this problem: During startup - Warning message: using .GlobalEnv instead of 'package:lattice' maybe this is the problem? How to fix it? Thanks ROnaldo -- I've been on this lonely road so long, Does anybody know where it goes, I remembe

[R] problem in font

Hi, This is a false problem. This problem is only in a .RData file, not a global R problem. Thanks ROnaldo -- I've been on this lonely road so long, Does anybody know where it goes, I remember last time the signs pointed home, A month ago. -- Carpenters, "Road Ode" -- > Prof. Ro

[R] problem in font

Hi, after a debian upgrade my dev.copy2eps dont work anymore. I have this message: Error in matchFont(postscriptFonts(family)[[1]], old$encoding) : unknown font I try to change the font family in my .Rprofile: setHook(packageEvent("graphics", "onLoad"), function(...) {

[R] Making R script to run in a console

Hi, is possible to make a R script to run under a console without open the R environment? Something like this example.R #!/usr/bin/R function(name="Put here your name") { print(name) } In a console I make ./example.R name="Ronaldo Reis Júnior" then program print my name. It is possible? Tha

[R] Calculation of lags

Hi, If I have the follow situation: A dependent variable (i.e. number of insects) that is affected by an independent variable (i.e. rain). The problem is that the measure of rain affect the population in other moment. So there exit a lag between the rain and the number of insects. Exist in R a

[R] calculating the spacial autocorrelation for poisson data

Hi, I have some count data. I try to calculate if exist spacial auto-correlation. In geoR I make: variog.p50 <- variog(p50,uvec=c(1:10)) but I think that the correct is using geoRglm I try: covariog.p50 <- covariog(p50,uvec=c(1:10)) But it is not the same calculation. I new in geostatistic,

[R] [OFF] The "best" tool for a space-temporal analyses?

Hi, I try to make an analyses to discover what is the time that an area begin to have spacial autocorrelation. And after, what is the number of individuals responsible for this autocorrelation. The main idea is to discover if exist a contamination of a quadrat from others quadrats and how is t

[R] is possible to use update with lmer models

Hi, the question on subject is possible? example: m1 <- lmer(y~1+(1|subject)) m2 <- update(m1,.~.+x) I try but dont work, exist any other method for this? Thanks Ronaldo -- He won a Toyota now, eh? -- palíndromo -- |> // | \\ [***] | ( õ

Re: [R] Problem with Crawley book example

Em Terça 04 Abril 2006 13:52, Dimitris Rizopoulos escreveu: > change the number of EM iterations and try again, this has worked for > me, some times, when I received the same error message as you did, > i.e., > > library(nlme) > model <- lme(height ~ seed, random = ~ time | rep / seed, control = >

[R] Problem with Crawley book example

Hi, I try to run the example of Crawley's Book on the page 661, but it fail, look > repmeasures <- read.table("../Packages/Crawley/data/repmeasures.txt",header=T) > attach(repmeasures) > rep <- as.factor(rep) > library(nlme) > model <- lme(height~seed,random=~time|rep/seed) Erro em lme.formula(h

Re: [R] manipulating weeks dates

Em Quarta 15 Março 2006 10:26, Gabor Grothendieck escreveu: > as.Date(paste(YEAR, WEEK, 0), "%Y %U %w") Hi, it works, but it use a year with 53 weeks, I need to use with 52 weeks, how to change this? Thanks Ronaldo -- Errigal Mountains <--> Tailoring manure -- anagrama -- |>

[R] manipulating weeks dates

Hi, I have these vectors: > WEEK <- rep(c(1:52),2) > YEAR <- rep(c(2000,1999),c(52,52)) How to make a vector of Date with weeks in years? I try as.date from survival package, but it dont work with weeks, just only with days, months etcs. Thanks Ronaldo -- "Realmente minha cidade e muito facul

[R] manipulating weeks dates

Hi, I have these vectors: > WEEK <- rep(c(1:52),2) > YEAR <- rep(c(2000,1999),c(52,52)) How to make a vector of Date with weeks in years? I try as.date from survival package, but it dont work with weeks, just only with days, months etcs. Thanks Ronaldo -- "Realmente minha cidade e muito facul

[R] simple problem

Hi, look this: > summary(família) Anacardiaceae Annonaceae Bombacaceae Cecropiaceae 2413 Chrysobalanaceae ClusiaceaeEuphorbiacaea Fabacea 114

[R] lme X lmer results

Hi, this is not a new doubt, but is a doubt that I cant find a good response. Look this output: > m.lme <- lme(Yvar~Xvar,random=~1|Plot1/Plot2/Plot3) > anova(m.lme) numDF denDF F-value p-value (Intercept) 1 860 210.2457 <.0001 Xvar1 2 1.2352 0.3821 > summa

Re: [R] Difficulties with for() {while(){}}

Em Qua 16 Nov 2005 11:30, jim holtman escreveu: > What do you mean by 'abort'? Does an 'error' occur that you want to catch? > If so, look at 'try'. Otherwise if it is testable, then just test for the > condition and restart. > Hi, this is not a real error, is a situation whitout a resolution. I

[R] Difficulties with for() {while(){}}

Hi, I have the follow function: function() { ## Init of function ... for(i in test) { ... while(j <= test2) { ... } } } The problem is that sometimes, naturally, the while is not possible to be resolved, and so the program abort. In this case I need that program ret

Re: [R] problems with for: warnings and segfault

Em Sex 11 Nov 2005 13:34, Peter Dalgaard escreveu: > > Look at > > xcoord[indice] <- xcoord+i > > The left hand side is a single element, the right hand side is a > vector so I don't think it means what I think you think it means. I'm stupid, of course. The correct is: > xcoord[indice] <- xcoor

[R] problems with for: warnings and segfault

Hi, I have two problem with a for looping using R Version 2.1.1 (2005-06-20) on a Debian Linux Testing. The first problem: warnings messages Look: > xcoord <- 5 > ycoord <- 5 > indice <- 1 > for(i in c(1:5)) {indice <- indice+1;xcoord[indice] <- xcoord+i; ycoord[indice] <- ycoord } Warning m

Re: [R] lme and lmer syntax

Em Ter 25 Out 2005 15:24, Doran, Harold escreveu: > There is an issue with implicit nesting in lmer. In your lme() model you > nest block/irrigation/density/fertilizer. In lmer you need to do something > like (I dind't include all of your variables, but I think the makes the > point) > > lmer(yield

Re: [R] lme and lmer syntax

Em Seg 24 Out 2005 18:08, Doran, Harold escreveu: > Ronaldo > > See the article on lmer pasted below for syntax. It is the only current > source documenting the code. In lmer(), the nesting structure for the > ranmdom effects is handled in a slightly different way. If your > observations are nested

[R] lme and lmer syntax

Hi, I have this: lme(y~x1+x2,random=~1|x1/x2) How to make this random effect using lmer? I try this: lmer(y~x1+x2+(1|x1/x2) But it dont work. Any idea? Thanks Ronaldo -- System halted! -- |> // | \\ [***] | ( õ õ ) [Ronaldo Reis Júnior

Re: [R] simple question

Em Sáb 22 Out 2005 10:07, Korbinian von Blanckenburg escreveu: > Its just a simple question I guess: > > I have a vector with missing information like > x<-c(0,1,31,131,NA,133,NA,310,NA,112,3,1,2,93) > > How can I make a vector like this no missing in it. I used the x[x<0] > commabd and tried some

[R] lmer and grouping fators

Hi, I make this model using lme m.lme <- lme(Glycogen~Treatment,random=~1|rTrt/Liver) How to make this using lmer? I try > m.lmer <- lmer(Glycogen~Treatment+(1|rTrt/Liver)) Erro em lmer(Glycogen ~ Treatment + (1 | rTrt/Liver)) : entry 0 in matrix[0,0] has row 2147483647 and column 214

[R] how to compare a null model with a model with NAs?

Hi, I make a model: m.null <- glm(y~1) m.comp <- glm(y~x1+x2+x3) I try make a anova like this anova(m.null,m.comp,test="F") The result: Erro em anova.glmlist(c(list(object), dotargs), dispersion = dispersion, : models were not all fitted to the same size of dataset The problem is that

[R] Doubt about Sweave

Hi, At this moment, I make my R output using comments and latex tags in R. Something like this: ### make a plot plot(y~x) dev.copy2eps(file="plot.eps") # @\includegraphics[width=\linewidth]{plot}@ ... After I include the file saved with this code in a latex document using the listings package.

[R] rkward

Hi, I'm testing rkward to use in my class. Now I'm using XEmacs, but for many people XEmacs is more complicated than R. Anybody here use this program? It use plugins, what are these plugins? Thanks Ronaldo -- Cada escola que se abre é uma cadeia que se fecha. -- |> // | \\ [**

[R] Change the mirror

Hi, Please, change the brazilian mirror http://www.termix.ufv.br/CRAN/ to http://www.insecta.ufv.br/CRAN/ in R homepage. Thanks ROnaldo -- Se dois homens no mesmo trabalho concordam o tempo todo, um deles é demais. Se discordam sem parar, então os dois são dispensáveis --Darryl F. Zanuck --

Re: [R] best way to fit a model

Em Sex 09 Set 2005 15:25, Achim Zeileis escreveu: > On Fri, 9 Sep 2005 14:55:04 -0300 Ronaldo Reis-Jr. wrote: > > Hi, > > > > I have some data that have this behaviour: > > |*** > > |* > > | * > > |* > > |

[R] best way to fit a model

Hi, I have some data that have this behaviour: | |*** |* | * |* | * | What is the best and simpler way to fit this in R? Thanks Ronaldo -- Ela pilotava um Continenal 2001 com ignição automática Magiclic... -- |> // | \\ [**

Re: [R] Doubt about nested aov output

Hi Spencer, Em Dom 04 Set 2005 20:31, Spencer Graves escreveu: > Others may know the answer to your question, but I don't. However, > since I have not seen a reply, I will offer a few comments: > > 1. What version of R are you using? I just tried superficially > similar things w

[R] Doubt about nested aov output

Hi, I have two doubts about the nested aov output. 1) I have this: > anova.ratos <- aov(Glicogenio~Tratamento+Error(Tratamento/Rato/Figado)) > summary(anova.ratos) Error: Tratamento Df Sum Sq Mean Sq Tratamento 2 1557.56 778.78 Error: Tratamento:Rato Df Sum Sq Mean Sq F

[R] help with a xyplot legend

Hi, I try to put a legend in a xyplot graphic. xyplot(y~x|g,ylim=c(0,80),xlim=c(10,40),as.table=T,layout=c(2,3), ylab="Número de machos capturados",xlab=expression(paste("Temperatura (",degree,"C)")), key=list(corner=c(0,0),x=0, y=0, text=list(legenda),lines=list(col=cor, lwd = espessura, lty=

[R] Rsquare from glmmPQL or another GLMM?

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 -- |> /

[R] How to calculate the correct SE in a nested or spliplot anova?

Hi! How to calculate the correct SE of mean in a nested or spliplot anova? Nested example: - m <- aov(Glycogen~Treatment+Error(Treatment/Rat/Liver)) > m Call: aov(formula = Glycogen ~ Treatment + Error(Treatment/Rat/Liver)) Grand Mean: 142. Stratum 1: Treatment Terms:

[R] How to automatic obtain the p value from a nested aov?

Hi, I try this nested anova: > m <- aov(Glycogen~Treatment+Error(Treatment/Rat/Liver)) > summary(m) Error: Treatment Df Sum Sq Mean Sq Treatment 2 1557.56 778.78 Error: Treatment:Rat Df Sum Sq Mean Sq F value Pr(>F) Residuals 3 797.67 265.89 Error: Treat

Re: [R] Error in F test on version 2.1.0

Em Seg 09 Mai 2005 10:53, Uwe Ligges escreveu: > Have you read the posting guide? > It asks you to specify a reproducible example > Please so so (both, read the psoting guide and specify the example). > > Uwe Ligges > tot <- c(16,20,23,18,10,18,25,20,24,24,29,12,23,23,21,23,11) > trat <-

[R] Error in F test on version 2.1.0

Hi, I make a upgrade to R 2.1.0 and in some analysis I give an error: anova(model,test="F") Analysis of Deviance Table Model: binomial, link: logit Response: landing/total Terms added sequentially (first to last) Df Deviance Resid. Df Resid. DevF Pr(>F) NULL16

[R] Keep R-mirror updated

Hi, I have a R mirror on www.termix.ufv.br/CRAN I try to keep it updated using rsync and cron. rsync -rzuav --progress --delete cran.r-project.org::CRAN /var/www/termix/mirrors/linux/CRAN/ & But all times the rsync abort with error: 95.99kB/s0:00:02 548960 77% 110.50kB/s0:00:01

Re: [R] Tool for update

Em Sex 25 Mar 2005 11:16, abunn escreveu: > I edited Rprofile to update everything on Tuesdays. I've been doing this > since 2.0 and I think I've had R running almost every Tuesday, which begs > the question of what I would be doing if R hadn't come into existence. > > In any case, It works pretty

[R] Anova for GLMM (lme4) is a valid method?

Hi, I try to make a binomial analysis using GLMM in a longitudinal data file. Is correct to use anova(model) to access the significance of the fixed terms? Thanks Ronaldo -- Todos somos iguais perante a lei, mas nao perante os encarregados de faze-las cumprir. -- S. Jerzy Lec -- |> // |

[R] R + MacOSX + Emacs(XEmacs) + ESS

Hi, I try to use Emacs or XEmacs with R in a MacOS X Panter without X11. Anybody can make this work? How? Thanks Ronaldo -- Um mapa-mundi que nao inclua a utopia, nao e digno de consulta, pois deixa de fora as terras a que a Humanidade esta sempre apontando. -- Oscar Wilde -- |> // |

[R] glm.nb

Hi, I make some simulations with rnbinom and try to test with glm.nb. But in some data set the glm.nb fail. Look: pop <- rnbinom(n=1000,size=1,mu=0.05) > glm.nb(pop~1,maxit=1000) Error in while ((it <- it + 1) < limit && abs(del) > eps) { : missing value where TRUE/FALSE needed look some pop

Re: R: [R] glm.nb stop on Error.

Em Tuesday 09 November 2004 12:12, Vito Muggeo escreveu: > if(class(model)=="try-error") Thanks it works fine :) Ronaldo -- O importante não é ganhar. O que importa é competir sem perder nem empatar. -- |> // | \\ [***] | ( õ õ ) [Ronaldo Reis Júnior

[R] glm.nb stop on Error.

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 }

[R] integrate function

Is possible to integrate this diferential equation: dN/dt = Nr(1-(N/K)) in R using the integrate() function? Or any other diferential equation? If yes, how? If no, anybody know any software on linux that make this? Inte Ronaldo -- O problema de ter os dois pés bem firmes no chão é que você

[R] Sciviews

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

[R] gui for MacOS X / Aqua

Hi, anybody know if Rcmdr or any other GUI like this work with R/Aqua? Thanks Ronaldo -- Let's organize this thing and take all the fun out of it. -- |> // | \\ [***] | ( õ õ ) [Ronaldo Reis Júnior] |> V [UFV/DBA-Entomologia

Re: [R] problem with restore and some .RData

Em Thursday 17 June 2004 17:08, Prof Brian Ripley escreveu: > The short answer is to add importFrom(stats, family) to MASS's NAMESPACE > and reinstall. Another way is to use > > library(stats) > > in your .Rprofile in that directory. Thanks, its works. > The real question is why your .RData is l

[R] problem with restore and some .RData

Hi, I have problem with the restore function in some .RData using R 1.9.0 Look the error: [EMAIL PROTECTED] RAnalise]$ R ... Error: object 'family' not found whilst loading namespace 'MASS' Fatal error: unable to restore saved data in .RData But if I load this .RData with the load() function all

[R] [OFF] program to estimate the best fit

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 -- |> // | \\ [

[R] Help on parameters

Hi, I have a follow analysis Trat1 = quantitative variable Trat2 = qualitative variable with 3 levels (A, B, C) Trat3 = qualitative variable with 3 levels (D, E, F) Resp = Response I try to get the parameters to compare with zero, so I make this model: glm(Resp~Trat1*Trat2+Trat1*Trat3-Trat1-1)

Re: [R] Problems with SJava instalation

Em Dom 14 Mar 2004 07:26, Ronaldo Reis Jr. escreveu: > Em Sex 05 Mar 2004 10:29, Jean Coursol escreveu: > > Hi, > > I try to execute ldconfig and the SJava dont work. I try to install with > gcc-3.3 (was gcc-2.95 in my first instalation). Where is the problem? > > Thanks &

Re: [R] Problems with SJava instalation

Em Sex 05 Mar 2004 10:29, Jean Coursol escreveu: > On 5 Mar, Ronaldo Reis Jr. wrote: > > Hi, > > > > I'm try to use the SJava package. The install is OK. In R I have this > > error message: > > > > -- > > > >> l

[R] Questions about spatial data

Hi, I have two questions about spatial data. 1) I have this grid x <- rep(c(1:8),c(rep(6,8))) y <- rep(c(1:6),8) I need to make some others grids with same length and different numbers of points whit a random spatial pattern. How is the best method to make this in R. 2) I need to produce som

[R] Books

Hi, I am seeking some books about spatial statistics and mixed models using R. I need books that approach binomial (proportion and binary) and poisson (count) data. The books will be used basically for ecologist. Would anybody have suggestions? Thank you Ronaldo -- Newton's Little-Known

[R] Problems with SJava instalation

Hi, I'm try to use the SJava package. The install is OK. In R I have this error message: -- > library(SJava) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library "/opt/lib/R/site-library/SJava/libs/SJava.so": libRSNativeJava

Re: [R] [OFF] Dataset for extra Crawley Chapter

Em Ter 14 Out 2003 12:01, Yang, Richard escreveu: > Ronaldo; > > You can download all datasets used in Crawley's book from his web > site. > > -- Richard > In the web site dont have the dataset for extra chapters, only for the printed book's chapters. Thanks Ronaldo -- Todo mundo precis

[R] [OFF] Dataset for extra Crawley Chapter

Hi, anybody have the dataset used in Gamma Errors chapter of the Crawley's books (An Introduction to Data Analysis using S-Plus). specifically the functionalresponse and the Density datasets. Thanks Ronaldo -- For every problem there is one solution which is simple, neat, and wrong. -- H. L. M

[R] different results depending of variable position.

Hi, I make an analysis and depending of the order of the variables, the significance change, look. m1 <- glm((infec/ntot)~idade+sexo+peso,family=binomial,weights=ntot) > anova(m1,test="F") Analysis of Deviance Table Model: binomial, link: logit Response: (infec/ntot) Terms added sequentially

[R] [Off] Ecological Modelling's Book Suggestion.

Hi, Anybody can suggest a good Ecological Modelling's Book. Thanks Ronaldo -- O céu sueco. -- palíndromo -- |> // | \\ [***] | ( õ õ ) [Ronaldo Reis Júnior] |> V [UFV/DBA-Entomologia] |/ \

Re: [R] how to calculate Rsquare

Can anybody send these articles for me? > NagelKerke, N. J. D. (1991) "A note on a general definition of the > coefficient of determination", Biometrika 78: 691-2. > > Cox, D. R. and Wermuth, N. (1992) "A comment on the coefficient of > determination for binary responses", The Ame

[R] R2HTML

Hi, I try to use R2HTML to make output. It is very good for this. I have this HTML output: -- HTMLStart(HTMLframe=F,echo=T) - > as.title("Comment one") Comment one > summary(1) Table of the summary 1 here > as.title("Comment two") Comment two > summary(2) T

[R] Error in par

Hi, I make a function and it work, but I get a error message in the follow step: image.plot(zlim = c(0,216),legend.width=0.2,legend.only=T, horizontal=T,col=palette(),nlevel=nivcor,offset=0.7, legend.shrink=1) The error is: Error in par(old.par) : invalid value speci

[R] how to calculate Rsquare

Hi, I have something like this: > x <- 1:10 > y2 <- 30+5*x+rnorm(x,sd=3) > y <- c(y1,y2) > x <- c(x,x) > plot(x,y) > x <- 1:10 > y1 <- 1+5*x+rnorm(x,sd=2) > y2 <- 30+5*x+rnorm(x,sd=5) > y <- c(y1,y2) > x <- c(x,x) > f <- factor(rep(c("a","b"),c(10,10))) > m <- lm(y~x+f) > anova(m) Analysis of Var

[R] doubt about graphics

Hi, I have an data like this x <- c(1:11) y <- seq(245,445,20) I make a plot plot(x,y) The first y value in y axis is 250 and the last y valeu is 450, the spacing between values are 50, so the y values showed are: 250, 300, 350, 400, 450 I need that values in y axis are: 240, 280, 320, 360,

[R] problem in pot with subset selection

Hi, I try to make a selective plot, but it make an error. The subset only work with length < 4, it is correct? > plot((ocorrencia/isca)~frag,subset=especieama==c("grupo1","grupo2","grupo3","AnoplotermesSp1","NeocapritermesOpacus"),pch=c(1,2,3,4,5),xlab="Área (ha)",ylab="Proporção de iscas ocu

Re: [R] [OFF] stepwise using REML???

Em Dom 22 Jun 2003 11:18, Douglas Bates escreveu: > > For nested design it may be very dangerous due the difference in > > variance structure, mainly in a splitplot design. ML make > > significative variables that REML dont make. > > It would be good to quote an example that shows this. I'm not su

[R] [OFF] stepwise using REML???

Hi, I know that is not possible make a stepwise procedure using REML in R, I can use ML for this. For nested design it may be very dangerous due the difference in variance structure, mainly in a splitplot design. ML make significative variables that REML dont make. I read an article that is m

[R] glmm and overall goodness of fit

Hi, exist in R any glmm function that have any tools for test for overall goodness of fit? Thanks Ronaldo -- O papel da impressora é sempre mais forte na parte picotada. -- | // | \\ [***] |> ( õ õ ) [Ronaldo Reis Júnior] | V [UF

[R] proportional splitplot with glmmPQL

Hi ALL, sorry by this repeated question, but I'm a student and I dont have money to buy any book more for now and I dont find any answer in network. I have a splitplot with binomial errors (proportion data). I try to make a model with glmmPQL. But my doubt are: Exit any method to compare model

Re: [R] Trying to make a nested lme analysis

%in% Treatment, say rTrT > by > > rats$rTrt = getGroups(~ 1 | Treatment/Rat, data = rats, level = 2) > > then fit the lme model as > > lme(Glycogen ~ Treatment, data = rats, random = ~ 1|rTrT/Liver) > > "Ronaldo Reis Jr." <[EMAIL PROTECTED]> writes: > &

[R] Trying to make a nested lme analysis

Hi, I'm trying to understand the lme output and procedure. I'm using the Crawley's book. I'm try to analyse the rats example take from Sokal and Rohlf (1995). I make a nested analysis using aov following the book. > summary(rats) Glycogen Treatment Rat Liver Min. :12

[R] [OFF]- Xemacs, delete key and ESS

Hi, I use R with XEmacs and ESS, when I'm editing a R file in XEmacs, the delete key work like the Backspace key. This behaviour is only in XEmacs + ESS + R file, it dont have this behaviour on a R section (M+R). Anybody know how to make a delete key to work like the delete key?? ESS 5.1.20-2

Re: [R] Problem to install RMySQL

the version of MySQL! > > My guess is that the version of R is the problem. > > On Sun, 30 Mar 2003, Ronaldo Reis Jr. wrote: > > I try to update my RMySQL but it dont work. > > > > Creating a new generic function for "print" in package > > RMySQL >

[R] Problem to install RMySQL

Hi, I try to update my RMySQL but it dont work. Creating a new generic function for "print" in package RMySQL [1] "print" Error in getProperties(ClassDef) : "slots" is not a valid slot for this object (or was mistakenly deleted) Execution halted ERROR: execution of package source for 'RMySQL'

[R] Ajusts in nls

Hi, I make an ajust using NLS. The ajust is: > library(nls) > modelo <- nls(y~a-b*exp(-c*x),start=list(a=100,b=90,c=0.1609)) > summary(modelo) Formula: y ~ a - b * exp(-c * x) Parameters: Estimate Std. Error t value Pr(>|t|) a 115.2527 2.9139 39.55 < 2e-16 *** b 118.6876 7.892

[R] Howto calculate MS from a unbalanced data

Hi, I make a nested lme analysis, it is ok. The estructure is: lme(y~x,random=~1|x/animal) It is no the exactly formula, but the teoretical formula. The real formula have a new variable with x/animal group. Now I try to construct a nested lme table like the traditional anova table to presenta

[R] Trying to make a nested lme analysis

Hi, I'm trying to understand the lme output and procedure. I'm using the Crawley's book. I'm try to analyse the rats example take from Sokal and Rohlf (1995). I make a nested analysis using aov following the book. > summary(rats) Glycogen Treatment Rat Liver Min. :12

[R] Problems with variable types.

Hi all, I have problems in a dataframe variables types. Look: from a loop function: for(...){ ... dados.fin <- rbind(dados.fin, c(L=j, A=j^2, Nsp=nsps, N=length(amosfin$SP), AmT="am",NAm=nam,

[R] Dataframe in loop

Hi, I try to make a dataframe in a loop function, but I dont have succeed. The function is something like this: for(i in c(10,12)) { expr for(j in c(1:2) { total <- c(1,2,3,4,5,6,7) nspf <- length(levels(as.factor(total))) fin <- data.frame(L=i,N=nspf)

[R] Nested analysis with non-normal errors

I have some doubt to make some analysis in R, I remember that in GLIM all analysis use glm. I have using the new Crawley's Book to help-me. How make a nested or splitplot analysis with poisson errors?? with normal errors use aov, and non-normal errors? Thanks fo all Ronaldo -- Whenever someone

[R] Help on a simple function.

Hi, I try to make a function that have two others functions inside. It is simple, but the problem is that functions inside use the same variable, but with different values. I try something like this: Teste <- function(Pdig(nlinhas),Ldig(nlinhas)) { Pdig <- function(nlinhas) { Tdig <- (15.50 +

[R] Receiving Rnews by mail

Hi all, I possible to receive Rnews automatically in my E-mail? Thanks Ronaldo -- O cerebro e um orgao maravilhoso. Comeca a funcionar assim que voce se levanta da cama e nao p ra ate voce chegar ao escritorio. -- Robert Frost -- | //|\\ [**

[R] testing slope

Hi all, I try to test a linear slope using offset. I have: > m2 <- glm(Y~X*V) > summary(m2) Call: glm(formula = Y ~ X * V) Deviance Residuals: Min1QMedian3Q Max -2.01688 -0.56028 0.05224 0.53213 3.60216 Coefficients: Estimate Std. Error t

[R] Help to make a function

Hi all, I have a simple code to produce some populations graphics. 10: DIM X(71) 20: INPUT "r=";R 30: INPUT "k=";K 40: INPUT "P(0)=";P(0) 50: FOR T=0 TO 69 60: P(T+1) = P(T) + R*P(T) * (1-P(T)/K) 70: NEXT T 80: GRAPH 90: COLOR 0 100: LINE (0,0)-(0,220) 110: LINE (0,0)-(220,0) 120: LINE (0,K)-(220

[R] [off] Statistical list

Hi all, Anybody know if exist a good statistical list? A list for discuss about different means to work in the same dataset. The profit and cost of each mean. Thanks for all Inte Ronaldo -- O fracasso e a oportunidade de comecar de novo inteligentemente -- | //|\\ [**

Re: [R] warning in binomial analysis

Em Thomas Lumley, escreveu: > > When you fit logistic regression models to fairly sparse data you can > often have a situation where for some combination of variables the > response variable is either all 0 or all 1. In that case the maximum > likelihood estimates for at least some of the coeffici

[R] [off-topic] The better analysis

Hi all, Firstly excuse-me by this off-topic question. I'm very confused to decide what is the better analysis to use. Whem the experiments is a tradicional design of manipulative experiments, it is easy to decide. But in natural observations is too complicated. I have 6 farms, each farm cultives

[R] doubt on survreg.

Hi all, I make a censored model using survreg from survival package. But the deviance values are too high, it is possible that must be a error in my model? > summary(completo) conc peso dv w temp Min. :0.7 Min. : 5.00 Min. : 1