similar to: match lists

Hi all, I have fitted a gee model with the gee package and included restricted cubic spline functions. Here is the model: chol.g <- gee(SKIN ~ rcs(CHOLT, 3), id=ID, data=chol, family=binomial(link="logit"), corstr="exchangeable") To extract the log odds I use: predict.glm(chol.g, type = "link") Now I want to compute the logg odds for specific CHOLT values

*Postdoc Position Computational Systems Biology 'Systems Bioinformatics: Computational Modelling Methods' f/m * * VU University Amsterdam, the Netherlands * *Research project * The goal of the project is to develop and implement computer methods for computational modelling of biological systems. This encompasses methods for computational model assembly, execution, analysis and

Dear All, Probably a very basic question, but can't seem to work my way around it. I want to which row has the maximum value. But what if the row names do not correspond with the row numbers. In the example below, you'll see that the max of example is row 4, but the name of row 4 is "9". How do I get R to return "9" as value, instead of 4. example <-

https://bugs.freedesktop.org/show_bug.cgi?id=90276 Bug ID: 90276 Summary: [NVE6] nouveau E[ PFIFO][0000:01:00.0] read fault at 0x000a5c0000 [UNSUPPORTED_KIND] from CE2/GR_CE on channel 0x007f329000 [unknown] Product: xorg Version: git Hardware: x86-64 (AMD64) OS: Linux (All)

Hi, I have the following function: TEST <- cbind(Pox0t30.stone.1990.2000, sqrt(apply(Pox0t30.stone.1990.2000[,2:11], 1, var))) In rhe resulting Object TEST the added column is titled: sqrt(apply(Pox0t30.stone.1990.2000[,2:11], 1, var) I would like to change it to a sensible variable name like PoxMean. How can I incorporate this into the above function? Ulrich --

Dear R users, Can anyone tell me how to get the p value out of the output from summary.manova? I tried all the methods I can think of, but failed. Many thanks Yu-Kang _________________________________________________________________ ¥ß§Y¥Ó½Ð MSN Mobile ªA°È¡G¦b±zªº¤â¾÷¤W¦¬µo MSN Hotmail http://msn.com.tw/msnmobile

Thanks to some help by Doug Bates (and the updated version of the Matrix package), I've refined my question about fitting nested and non-nested factors in lmer(). I can get it to work in linear regression but it crashes in logistic regression. Here's my example: # set up the predictors n.age <- 4 n.edu <- 4 n.rep <- 100 n.state <- 50 n <- n.age*n.edu*n.rep age.id

Dear XpeRts, I prepared a no qoute Character string by the following command s<-noquote(paste (b1, collapse=",")) where, b1 is the vector of 24 intergers. > dput(b1) c(1L, 2L, 6L, 7L, 12L, 16L, 17L, 20L, 21L, 23L, 25L, 34L, 46L, 48L, 58L, 64L, 65L, 68L, 82L, 97L, 98L, 101L, 113L, 115L) > dput(s)

I'm having a problem reading data to set control totals for a dataframe. I want to adjust a dataframe based on a 2-d table of values, which I get by using : > CurrentX1Sums <- as.matrix(xtabs(~tripid_nu+lineon, data=SurveyData)) > CurrentX2Sums <- apply(CurrentX1Sums, 1, sum) I've created a .csv file with new (target) sums that looks like this: tripid_nu Warner

Hello I have two DSL modems witch are connected to my isp, in future my boss want to buy another connection via DSL modem. Then i will have 3 DSL modems. With every DSL modem i get 3x8 IP`s (netmask 248) from my ISP, now i have question how to configure gateway wich tools should i use. Because ip route and next hop via. wich i use now makes his work fine but with new kernels there is an error

Upgraded to Samba 3.5.6 from 3.5.2. Samba is running as a standalone server on a box that was upgraded from Slackware 13.1 to Slackware 13.37. All Win XP Pro and Win 7 Pro clients can connect to the workgroup and browse shares and folders exactly as before; however, many complain of intermittent pauses, slowdowns and "Not Responding" messages when traversing directories. The samba server

Hello frieds, do I need to do 'net rpc join -U Administrator' (to selfjoin the localhost to the domain) if I have samba3.0.4 running in ROLE_DOMAIN_PDC and with ldapsam-Backend and want to use winbindd for ntlm_auth squid authentification? -- Best regards, Malte mailto:malte.woelky@gmx.de _________________ Malte Woelky -=[SkyNet]=-

I guess this has a simple solution: I have matrix 'mat1' which has row and column names, e.g.: A B C row1 0 0 0 row2 0 0 0 .... rown 0 0 0 I have a another matrix 'mat2', essentially a subset of 'mat1' where the rownames are all in 'mat1' e.g.: B row3 5 row8 6 row54 7 I want to insert the values of matrix mat2 for column B (in reality it could be some or

Hi, I have been studying the R manual on lists but cannot seem to create a simple R list in C - I keep on getting "'dimnames' applied to non-array" whenever I attempt to assign names to the list elements: Wanted output a list structure something like [ type="Bid", price=2.0, volume=1000 ] I can get up to a list of [ "Bid2, 2.0, 1000 ] But for the life of

Hi, I'm experimenting with the mail-filter plugin but whenever I add anything to a message, it always results in an corrupted cache: from the logs: Dec 14 16:09:55 ciphermail-webmail dovecot: imap(martijn=djigzo.com at ciphermail.private): Error: Corrupted index cache file /var/vmail/martijn=djigzo.com at ciphermail.private/dovecot.index.cache: Broken physical size for mail UID 76 Dec 14

I have a problem deploying FastCGI for use with RoR. OS: FreeBSD 5.4 Apache: 2.0.55 FastCGI: 2.4.2 ###### I compiled a new mod_fastcgi.so and copied it to /usr/lib/apache/mod_fastcgi.so ###### My httpd.conf looks like this: LoadModule fastcgi_module /usr/lib/apache/mod_fastcgi.so <IfModule mod_fastcgi.c> AddHandler fastcgi-script .fcgi </IfModule> ###### The testscript

I would like to fit a hierarchical regression model from Witte et al. (1994; see reference below). It's a logistic regression of a health outcome on quntities of food intake; the linear predictor has the form, X*beta + W*gamma, where X is a matrix of consumption of 82 foods (i.e., the rows of X represent people in the study, the columns represent different foods, and X_ij is the amount of

hello, i try to plot 3 graphs which have the same x.axis underneath each other. i'd like the plots to be aligned without margings between the boxes and draw a single x axis beneath the lowest plot. i managed to get the alignment by setting par(mar), but the middle box is stretched and i cant't figure out how to get around this. par(pin) was my guess, but this doesn't do the job

Dear R listers, I am trying to produce a simple (for a stata user) barplot with 4 countries on the x axis, each country observed in 2 subsequent years and 3 variables. Basically, I should have three bars for each year for each country. I am attaching the chart I made in Stata, but I am not sure you'll manage to see it! I did the following: #here I create the data-set TUSE2. The vectors mw, st

Hi all, I've been toying with xapian (mostly using the Python bindings) and I think I've hit a bug in the TermIterator::skip_to() method (or maybe in QuartzAllTermsList::skip_to()). I've attached a c++ source file that demonstrates the issue. In short, if you have a WritableDatabase, ask for the all-terms TermIterator with db.allterms_begin(), and then skip_to() a word that is itself