similar to: predict.lm() question

Dear all, Is there a less cumbersome way to rename a column by name (as opposed to index) than -- names( X)[ names[ X] == "bob"]<-"sue" ? A semi-related question: how does one get the index of a column by name, something along the lines of col.index( X, "sue") ? Chip Barnaby --------------------------------------------------------- Chip Barnaby

Hello, I would like to create a subscriptable collection (presumably a list) of lm() models. I have a data frame DX containing 6 groups of data. The general idea is (NOT RUN) ... for (i in 1:6) { DXS = subset( DX, <whatever>); LMX[ i] = lm( <formula>, data = DXS); } Now access model results by subscript ... e.g. coefficients( LMX[ 2]). Or would it be [[ 2]]? I

Greetings, Q #1 ------ How does one combine data frames by row ... no cleverness a la merge(), just add rows. For example, given A with 20 rows and B with 30 rows, I want C = combine( A, B) having 50 rows. Columns having matching names should be filled from both (all) sources with suitable coercion, unmatched would get NA in rows whose source does not have that column. Q #2 ------ Given

Dear all -- I have a data frame containing data related to heat gain through windows. The general form is ... Key ProfA IAC <many other numeric columns> AAA 0 .7 AAA 10 .6 AAA 0 .66 AAA 20 .45 (more AAA rows) (then AAB rows) 'Key' identifies the physical configuration ... rows with a given Key contain data for same window under various

Dear all -- I want to compute weighted.mean() for grouped rows. Data frame extract is just below. For each Key, I want the mean of IAC weighted by Wt. DP0[1:20,] Key IAC Wt 2 C3-PD030020050.PD030020050.3.12.3.0 0.765 0.8590000 3 C3-PD030020050.PD030020050.3.12.3.0 0.764 0.8449651 4 C3-PD030020050.PD030020050.3.12.3.0

Dear R-ers, My aggregation saga continues. Using the following sequence, I can calculate any statistic for row groups and merge the result back to all associated rows ... > WM = by( D60, D60[ "KeyProfA"], FUN=function(x) weighted.mean( x$IAC, x$Wt)) > D60$IAC.WM = as.numeric( WM[ D60$KeyProfA]) > class( WM) [1] "by" Questions ... 1) Is this a reasonable way

Hi, I'm a complete linux newbie, let alone wine, easy-peasy and netbooks. Having said that, I'm looking at buying the eee pc 701 sd and putting easy-peasy on it. However the reason I need a netbook is to run a cricket scoring program that is only available in windows. I didn't want to fork out for the netbook until I knew the cricket software will run, so I tried the following: 1.

Hi, We have a problem with Syslinux 2.07 and APM on an SMP system. When either network booting with pxelinux.0 or CDROM booting with isolinux we get a kernel error message: Jan 7 09:59:51 ddcserver1 kernel: kernel BUG at apm.c:1756! Jan 7 09:59:51 ddcserver1 kernel: invalid operand: 0000 Jan 7 09:59:51 ddcserver1 kernel: Jan 7 09:59:51 ddcserver1 kernel: CPU: 1 Jan 7 09:59:51 ddcserver1

Can anyone here help me out with a client side Outlook issue? I know outlook is a horrible imap client :( We intermittently have a problem where mail gets "stuck" in Outlook users' outboxes. Might have anything to do with dovecot? Is there anyway I can troubleshoot this server side? Are there any settings I can tweak? We have a new and an old dovecot server (see below) and both

Is there a limit to the number of columns that a data.frame can have? For example, can I read.csv() a file that has 1000 columns and 10,000 rows, will it break or is it limited by available memory. ................................................................................... Philip Sobolik 781-862-8719 x111 Wrightsoft Corporation 781-861-2058 fax 394 Lowell

Hi, i just want to know whether there are any alternative file systems for HP SFS. I heard that there is Cluster Gateway from Polyserve. Can anybody plz help me in finding more abt this Cluster Gateway. Thanks and Regards, Ashok Bharat -----Original Message----- From: lustre-discuss-bounces at lists.lustre.org on behalf of lustre-discuss-request at lists.lustre.org Sent: Tue 2/12/2008 3:18 AM

The PDF file on the web, which is an appendix on nonlinear regression associated with the CAR book, is very nice. When I ran through the code presented there, I found something odd. The code does a certain model in 3 ways: Vanilla NLS (using numerical differentation), Analytical derivatives (where the user supplies the derivatives) and analytical derivatives (using automatic differentiation). The

Hi, I am having a hard time getting what the summary.lm-output for an ANOVA / ANCOVA means. Examples I find always seem to deal with simpler cases than what I meet in my data. My main problem is understanding the output when getting significant INTERACTION TERMS (what never occurs in examples :(). The following is the output after summary.lm(ancova) where "week" is continuous,

Dear R users, I have been having problems getting believable estimates from anova on a model fit from lmer. I get the impression that F is being greatly underestimated, as can be seen by running the example I have given below. First an explanation of what I'm trying to do. I am trying to fit a glmm with binomial errors to some data. The experiment involves 10 shadehouses, divided between

Dear List, i am trying to produce a 3d plot using wireframe using the code: wireframe(Residuals_FD ~ Elevation * Temperature, data = data2, scales = list(arrows = FALSE), drape = TRUE, colorkey = TRUE) As you can see when the code (using the data below) is run the plot area is set-up correctly but the actual surface is missing? Any help would be greatly appreciated. Chris #data Elevation

estimados, ajusté una superficie de respuesta con el paquete "rsm" de R y hago un grafico 2D de la superficie de respuesta con contour(), pero como el argumento de esta función es la superficie de respuesta ajustada y ésta se ha ajustado codificando las variables, de acuerdo a la recomendación de la documentación del paquete rsm, obtengo un gráfico con los ejes en unidades codificadas

Hi, I have a question about estimating IRT models by using nlme, not just rasch model, but also other models. Behavior Research Methods <http://www.springerlink.com/content/1554-351x/> Volume 37, Number 2 <http://www.springerlink.com/content/1554-351x/37/2/>, 202-218, DOI: 10.3758/BF03192688 Using SAS PROC NLMIXED to fit item response theory models (2005). Ching-Fan

flight 17901 linux-linus real [real] http://www.chiark.greenend.org.uk/~xensrcts/logs/17901/ Regressions :-( Tests which did not succeed and are blocking, including tests which could not be run: test-amd64-i386-qemuu-rhel6hvm-intel 4 xen-install fail REGR. vs. 12557 test-amd64-i386-pv 4 xen-install fail REGR. vs. 12557 test-amd64-amd64-pv 4

flight 17916 linux-linus real [real] http://www.chiark.greenend.org.uk/~xensrcts/logs/17916/ Regressions :-( Tests which did not succeed and are blocking, including tests which could not be run: test-amd64-amd64-pv 4 xen-install fail REGR. vs. 12557 test-amd64-amd64-xl 4 xen-install fail REGR. vs. 12557 test-amd64-amd64-xl-win7-amd64 4

flight 18150 linux-linus real [real] http://www.chiark.greenend.org.uk/~xensrcts/logs/18150/ Regressions :-( Tests which did not succeed and are blocking, including tests which could not be run: test-amd64-i386-rhel6hvm-amd 7 redhat-install fail REGR. vs. 12557 test-amd64-i386-pair 17 guest-migrate/src_host/dst_host fail REGR. vs. 12557 Regressions which are regarded as