search for: alr3

...ying to find datasets that will give me residuals, after applying the lm function, with no normality, non linearity, and heteroscedacity so I can try to exemplify those cases in the linear regression model. Can you give any advice on what datasets would be appropiate? I can?t use the ones in the alr3 package because those have already been seen in class. Thank you very much :-) natorro -- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean.

...bles. I wanted to ask what is the preferred way to find out if predictor and/or response needs to be transformed and if yes how (log-transform?). I found a procedure in "A modern approach to Regressoin in R" (Sheather, 2009): There they suggest an approach with the function bctrans from alr3...but it seems that it is deprecated. So what is the best way (box-cox test) find the best transformation for predictor and response simultaneously? AFAIK boxcox from MASS is used only used for transformation of the predictor? Thank you very much Johannes --

...ata.frames with 10 variables and 10,000 rows each. The command in question is: temp = merge (demo.2,ph,by="seqn",all.x=TRUE) in 2.9.0, this command took 3.3 seconds. in 2.9.1, it took 35.8 seconds. I have reverted back to 2.9.0. Additional packages loaded are: library(Hmisc) library(alr3) library(epicalc) library(ggplot2) library(lattice) library(reshape) library(survey) library(car) thanks very much for all the effort. R is wonderful.

...ny changes to old functions, including a change in naming conventions, using "camel-case" rather than period-separated names: for example, av.plots() is now avPlots(). For the time-being, old function names, such as av.plots(), are available as deprecated aliases. A few functions from the alr3 package have also been rewritten, renamed, and are now part of car. The alr3 package will be updated later in the year. Regards, John -------------------------------- John Fox Senator William McMaster Professor of Social Statistics Department of Sociology McMaster University Hamilton, Ontario...

Dear all, I'm interested in testing 2 nonlinear restrictions on coefficients of a nls object. Is there a package for doing this? Something in the lines of `test(nls object, res=c("res 1","res 2"),...)' I only found the function delta.method in the alr3 library that calculates the se of a singleton nonlinear restriction of a nls object using the delta method. Thanks in advanced for your help and suggestions. David

Dear all, Is there a way to retrieve standard errors from nls models? The help page tells me that arguments such as se.fit are ignored... Many thanks and best wishes Christoph -- Dr. rer.nat. Christoph Scherber University of Goettingen DNPW, Agroecology Waldweg 26 D-37073 Goettingen Germany phone +49 (0)551 39 8807 fax +49 (0)551 39 8806 Homepage http://www.gwdg.de/~cscherb1

...ny changes to old functions, including a change in naming conventions, using "camel-case" rather than period-separated names: for example, av.plots() is now avPlots(). For the time-being, old function names, such as av.plots(), are available as deprecated aliases. A few functions from the alr3 package have also been rewritten, renamed, and are now part of car. The alr3 package will be updated later in the year. Regards, John -------------------------------- John Fox Senator William McMaster Professor of Social Statistics Department of Sociology McMaster University Hamilton, Ontario...

CRAN (and crantastic) updates this week New packages ------------ Updated packages ---------------- New reviews ----------- This email provided as a service for the R community by http://crantastic.org. Like it? Hate it? Please let us know: cranatic at gmail.com.

...- Ignored: > x=c(12,24,13,11,23,10,9,17,11,14,18) > y=c(10,11,9,8,5,12,11,21,12,13,14) then I load the press source code, I run it using the x and y values defined before and I get the following > press(x,y) Loading required package: drc Loading required package: alr3 Loading required package: lattice Loading required package: magic Loading required package: abind Loading required package: MASS Loading required package: nlme Loading required package: plotrix 'drc' has been loaded. Please cite R and 'drc' if used for a publicatio...

...atim (label=lm) 7 : echo term verbatim (label=diag) 8 : echo term verbatim (label=tukey) Error: chunk 8 (label=tukey) Error in parse(text = x) : unexpected symbol in "yield ~ variety + block + preds(yield+variety +block+preds.sq)sq" It failed in: <<tukey>>= require(alr3) tukey.nonadd.test(ov.lm) @ At that point I ran debug() on the call: debug(tukey.nonadd.test) tukey.nonadd.test(ov.lm) debugging in: tukey.nonadd.test(ov.lm1) debug: { envir <- environment(formula(m)) dd <- eval(m$call$data, envir) subs <- eval(m$call$subset, envir) wg...

Dear R users, How can I do a regression analysis in R where there is more than one observation per x value? I tried the example in Sokal&Rohlf (3rd edn., 1995), page 476 ff., but I somehow couldn??t find a way to partition the sums of squares into "linear regression", "deviations from regression", and within-groups SS. I tried

Dear all, Is there a fairly general R implementation of the delta method that uses numerical derivatives? I realise that the delta method has been implemented using symbolic derivatives (e.g. alr3::delta.method, emdbook::deltamethod, msm::deltamethod and survey:::nlcon), however possibly non-linear estimators using the delta method with numerical derivatives can be quite useful (e.g. predictnl in Stata and the estimate and predict statements for proc nlmixed in SAS). I would like something...

...ta can be found following the given links: either http://www.sendspace.com/file/dnugwc or http://rapidshare.com/files/420569660/sel.day.txt I want to use the Box-Cox transformation to normalize the response as well as the predictor variables. The scatterplot looks like this: library(zoo) library(alr3) load("sel.day.txt") sel.p1<-window(sel, start=as.POSIXct("2008-04-05"), end=as.POSIXct("2009-04-01")) pairs(~v.obs+ snow+ HH6.1+ Q.Enz+ pcpt+ qd,data=sel.p1,gap=0.4,cex.labels=1.5) In Sheather: "A Modern Approach to Regression with R" the function...

Hi all! I try to estimate a statistic of the form: (x1-x2)/(y1-y2), where x1,x2,y1,y2 represent variable means, so each has an estimate and standard error associated with it. How is it possible to estimate the mean and the variance of this ratio? Thank you! [[alternative HTML version deleted]]

I apologize if this comes across as confusing. I will try to explain my situation as best I can. I have R bootstrapping my growth data for fish. It's resampling my database of age and length data and then produces several new datasets for me. In this case, it's resampling my data to create three new datasets of age and length data. Here is my code with my original data called

...r alia: ... coefficients (intercept) a 1.2 3.4 Packages installed: acepack ace() and avas() for selecting regression transformations adlift An adaptive lifting scheme algorithm akima Interpolation of irregularly spaced data alr3 Methods and data to accompany Applied Linear Regression 3rd edition car Companion to Applied Regression coda Output analysis and diagnostics for MCMC drc Analysis of dose-response curves EbayesThre...

...=TRUE) BestDF$max.iter<-factor(BestDF$max.iter,ordered=TRUE) BestDF$instance<-factor(BestDF$instance,ordered=FALSE) BestDF$migr.freq<-factor(BestDF$migr.freq,ordered=FALSE) anova.ax3<-aov(best~(comm.strat+migr.freq+max.cpu+instance+max.iter) ^2,data=BestDF) summary(anova.ax3) library(alr3) residual.plots(anova.ax3,ask=FALSE) -----= END R SCRIPT =----- Thanks in advance for any help you could provide. -------------------------------------------------------------------- Max MANFRIN Tel.: +32 (0)2 650 3168 IRIDIA - CoDE, CP 194/6...

...C_CTYPE=Dutch_Belgium.1252;LC_MONETARY=Du tch_Belgium.1252;LC_NUMERIC=C;LC_TIME=Dutch_Belgium.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] drc_1.4-2 plotrix_2.4-2 nlme_3.1-89 MASS_7.2-41 lattice_0.17-6 [6] alr3_1.1.7 loaded via a namespace (and not attached): [1] grid_2.7.0 d<-data.frame(dose=c(2.00e-05, 4.00e-06, 8.00e-07, 1.60e-07, 3.20e-08, 6.40e-09, 1.28e-09, 2.56e-10, 5.10e-11, 1.00e-11, 2.00e-05, 4.00e-06, 8.00e-07, 1.60e-07, 3.20e-08, 6.40e-09, 1.28e-09, 2.56e-10, 5.10e-11, 1.00e-11), resp...

...ot;,"micEcon","aaMI","abind","accuracy","acepack","actuar","adapt","ade4","adehabitat","adlift","agce","agsemisc","akima","AlgDesign","allelic","alr3","amap","AMORE","AnalyzeFMRI","aod","ape","mimR","minpack.l m","misc3d","mitools","mix","mixreg") TEMP <- c(TEMP,"mlbench","mlica","mlmRev","mml...