similar to: Odds ratios for categorical variable

Hello, Is there an easy function for switching list to matrix. My list is of genetic distances between species pairs: A A 0 A B .5 A C .25 B C .5 and I want a distance matrix such as: A B C A 0 .5 .25 B .5 0 .5 C .25 .5 0 for use in a mantel test. Thank you for the help! cheers, charlie -- ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Charles G. Willis Department of Organismic

Hello, How can one write the output of a GLM/GEE model to a text file, such that the results appear in text file in the same (or similar) format they appear in the R console? For instance, I have the following model: Traitresult <- compar.gee(Y~X, data = data, family = "binomial", phy = tree, scale.fix = TRUE, scale.value = 1) How would I write "Traitresult" to a text

Hello, I am trying to test for non-linearity in a set of non-parametric data. Furthermore, I would like to do so in a multi-variate model. Frankly, I don't even know if this is possible, and if it is possible, I don't know if it is implemented in an R package. Any advice would greatly appreciated. Thank you! Charlie -- ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Charles G.

I have been trying to run a logistic regression using a number of studies. Below is the syntax, error message & data. Any advice regarding what I am doing wrong or solutions are appreciated, regards Bob Green > logreg <- read.csv("c:\\logregtest.csv",header=T) > attach(logreg) > names(logreg) [1] "medyear" "where" "who"

Hi there, I've looked through the very helpful advice about adding fitted lines to plots in the r-help archive, and can't find a post where someone has offered a solution for my specific problem. I need to plot logistic regression fits from three differently-sized data subsets on a plot of the entire dataset. A description and code are below: I have an unbalanced dataset

I really like these Trellis graphics but how do I get this code to skip the missing? logreg<-read.csv("logreg.csv", header=TRUE, sep=",", na.string=" ") attach(logreg) bwplot(yesno~bc_pcb_tot |varlist, data=logreg, main="Box Cox PCB transformation", auto.key=TRUE, fontfamily = "HersheySans" ) Dean Sonneborn M.S. Public Health Sciences *

Hi all, I am trying to bootstrap a small data set into 1000 "pseudodatasets" and then run an ANOVA on each one. Can anyone provide guidance on how I could do this? Thank you. -Dan Janes ************************************************ Dan Janes, Ph.D. Harvard University/OEB 26 Oxford St. Cambridge, MA 02138 Office: 617-496-2375 Fax: 617-495-5667 Email: djanes at oeb.harvard.edu

Ciao, ho aperto con R un file di classe data frame con 15000 righe e 29 colonne. Nella console però sono visualizzate solo la prime e l'ultima colonna e le ultime 8000 righe circa. E' possibile una visualizzazione completa? Grazie Sabrina [[alternative HTML version deleted]]

I have been working on the following code but keep getting an err message. My current thinking is that the problem is on the indexing but do not know how to fix it. Any help please? ungulate <- read.csv("Ungulate.csv",row.names=1) ungulate <- as.matrix(ungulate);colnames(ungulate)<-NULL;rownames(ungulate)<-NULL habitat <- read.csv("Ungulate_vegetation.csv")

Hi, I was trying to plot the logistic regression from a regression "logreg" I just ran. I downloaded the "car" package from the R website and went to Packages -> install package from local zip file I checked in my library file and the package is there. I restarted R. I then ran the command: reg.line(logreg,col=palette()[2], lwd=2, lty=1) And I get the error: Error: could

Dear all, <BR><BR>I was used to work on R1.6 and I have now passed on R1.6.2 but I can''t read my <BR>file (and that is a big problem!!).<BR>I made a data sheet with some <BR>spreadsheet in Excell, and save it as separeted by tab .txt.<BR>I write in R<BR>read.table ("file.txt",h=T,sep="/t",dec=",")<BR>But R

Hi ! I'm runing a linux redhat 7.2 box with samba 2.2.4. I've downloaded the samba source in order to install it with the needed compilation options (--with-smbwrapper --with-automount --with-smbmount --with-pam ...). Samba is now up and running but I would like to build an rpm with this compiled version. Any idea of how to transform this samba compiled version to a rpm package ?

[This email is either empty or too large to be displayed at this time]

Hello R-users, I have code that gives me the important variables from an analysis. I need to input these variables into a different analysis. To do this, I need to modify them slightly... 1) remove all numbers at the end of the variables, 2) remove all periods. I tried to do it with the awkward code below. It works to remove all the numbers, but when I try to remove the period everything

All, I have three columns of data: id, date, hab. I am trying to set up a matrix that has the id as the rows, date as columns, and the hab value as the data values. Each id/date combination can only have one hab value. I would like for it to look something like this" date 1, date 2, date 3. id 0001 2 1 2 0002 3 5 6 I don't know how to set

Hi, I'm working on the codes below however every time I run them when they get to OpenBUGS I keep getting the error message: array index is greater than array upper bound for hab. Any help would be greatly appreciated, Suzie Codes: ungulate <- read.csv(file.choose ()) #ungulate ungulate <- as.matrix(ungulate);colnames(ungulate)<-NULL;rownames(ungulate)<-NULL

Consider this data structure (df1) ... Group Year PctProf FullYr 1 Never RF 2004 87 88 2 Cohort 1 2004 83 84 3 Cohort 2 2004 84 86 4 Cohort 3 2004 87 87 5 Cohort 4 2004 73 74 6 Never RF 2005 85 86 7 Cohort 1 2005 81 82 8 Cohort 2 2005 81 81 9 Cohort 3 2005 78 79 10 Cohort 4 2005 72 74 11

Hi there: recently i try to use LogicReg package for a tree model(logistics fit ) . i list my code and error below: > dim(model.dat) [1] 48000 745 > fit1 <- logreg(resp = model.dat[,745], bin=model.dat[, 9:700], type = 3, select = 3, ntrees = c(1,2), nleaves=c(1,7), ) Insufficient declaration LGCn1MAX in logreg() is 20000 LGCn1MAX should be at least 48000 Please fix and

On 08/05/15 12:49, Mario Pio Russo wrote: > Hi all > > just going back to this, Steve can you please expalin where did you take > this number from? 11644473600 > > that seems not to work correctly on my system so I was lookin g how to > tailor down the right value > > thanks > > >

great thanks Steve ___________________________________________________________________________________________ Mario Pio Russo, System Admin SWG IT Services Dublin, Phone & FAX: +353 1 815 2236, eMail: mariopiorusso at ie.ibm.com IBM Ireland Product Distribution Limited registered in Ireland with number 92815. Registered Office: IBM House, Shelbourne Road, Ballsbridge, Dublin 4 (Embedded