similar to: Bootstrap direction

Hi All, Two weeks ago, I received instructions from the R list about creating pseudodatasets by bootstrapping an existing dataset 1000 times. I received code that bootstrapped my dataset and ran an ANOVA on each pseudodataset, producing a histogram of F-values. Here is the code I received that produces this F-histogram. X <- data.frame(b=rnorm(2162),c =factor(rep(1:47,each=46))) ##

Dear R-users: How does one calculate in R the odds ratios for a CATEGORICAL predictor variable that has 4 levels. I see r-help inquiries regarding odds ratios for what looked like a continuous predictor variable. I was wondering how to get the pairwise odds ratios for comparisons of levels of a categorical predictor variable. I can't seem to get the correct output using: >

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.

Hi, I would like to know how to capture the console output from running an algorithm for further analysis. I can capture this using capture.output() but that yields a character vector. I would like to extract the actual numeric values. Here is an example of what I am trying to do. fr <- function(x) { ## Rosenbrock Banana function on.exit(print(f)) x1 <- x[1] x2 <- x[2]

Hello James, Parth, Following our discussion on IRC and on code review, the way FeatureVector class works needs some discussion. Presently, the FeatureVector class is defined as follows, with a fixed number of feature count (19): class FeatureVector::Internal : public Xapian::Internal::intrusive_base{ friend class FeatureVector; double label; double score;

> > > > The approach I was thinking would look something like this: > > * instead of Features, which is really a namespace implemented as a > class, we separate out the calculation of the different features > into distinct subclasses of Feature, whose only job is to calculate > a single feature. Currently the FeatureManager calls these (via >

Dear All, I've just installed Samba version 1.9.18p10 on one of our HP-UX boxes. The unix box is advertising itself with two WINS names. netbios name = ukswi0103 netbios aliases = ukswi0104 in smb.conf I've been playing with the idea of changing the behaviour of the server based on what the client calls it. What I'd like to do is have samba do security = user if the client

Is there a painless way to find the names of all packages on CRAN that "Depend" on a specified package? url: www.econ.uiuc.edu/~roger Roger Koenker email rkoenker at uiuc.edu Department of Economics vox: 217-333-4558 University of Illinois fax: 217-244-6678 Champaign, IL 61820

Hello I have 4 NT server : so 4 NT Domain and 1 Linux Box with samba 1.9.18p4-50.4 in smb.conf, I configure netbios aliases = TEST1 TEST2 TEST3 TEST4 sercurity = server include = /etc/samba/server.%L ... and in /etc/samba/server.TEST1 password server = NTSERV1 and in /etc/samba/server.TEST2 password server = NTSERV2 but when I read the nmbd and smbd logs , %L is replaced by a null string

Hello, I am trying to preform a comparative analysis using the partial mantel method. However, i am having trouble creating compatible matrices. I can create matrices from my numerical data using the 'distance(x, method)' function from ecodist, but the row/column names are out put as numbers. While when i create a phyo-distance matrix using 'cophenetic(x)' from APE, the row/column

Hello, I am running the program COMPAR.GEE within the package APE. My dependent variable is binomial, while my independent variable is a multi-state categorical variable. The output reports an estimate for each state of the independent variable except the first one. For example, for the variable X with 3 states, the output is: intercept (estimate) X2 (estimate) X3 (estimate)

I am writing regarding serious problems that a number of researchers at a workshop held by the Smithsonian Tropical Research Institute are encountering with the most recent R version for Macs (the patched version of 2.1.0 (2.1.0a) that you indicate should be used over 2.1.0). The following problems are not specific to my computer but are being experienced by ALL Mac users at this workshop

Hello, I am trying to run the APE program COMPAR.GEE with a model containing a categorical response variable and a mixture of continuous and categorical independent variables. The model runs when I have categorical (binary) response and two continuous independent variables (VAR1 and VAR2), but when I include a categorical (binary) independent variable (VAR3), I receive the following output with

ESTE EMAIL ? AUTOM?TICO. POR FAVOR, N?O ENVIE UMA REPOSTA, POIS ELA N?O SER? RECEBIDA. N?o foi poss?vel entregar a sua mensagem. O destinat?rio n?o recebeu a sua mensagem. Veja a seguir a descri??o do motivo. <steve@ig.com.br>: A caixa postal do destinat?rio excedeu o limite de capacidade de armazenamento. maildrop: maildir over quota. --- Below this line is a copy of the message.

Hi, I have some very noisy, relatively sparse data; a biological response of roughly ~8 subjects at ~8 times points). I've been following the data trend using a lowess line, over-plotted with several values of bandwidth, 'f <- seq(0.3, 0.9, by=0.1)'. At this point, we have no models for these data. I wonder if there is any way under R to assign some sort of confidence interval to

Using Ulrik?s example data (and assuming I understand what is wanted), here is what I would do: ex.dat <- c("FName: fname1", "Fval: Fval1.name1", "Fval: ", "FName: fname2", "Fval: Fval2.name2", "FName: fname3") tst <- data.frame(x = ex.dat, stringsAsFactors=FALSE) sp <- strsplit(tst$x, ':', fixed=TRUE) chk <-

@Don your solution does not solve Vijayan's scenario 2. I used spread and gather for that. An alternative solution to insert mising Fval - picking up with Don's newtst - is newtst <- c("FName: fname1", "Fval: Fval1.name1", "FName: fname2", "Fval: Fval2.name2", "FName: fname3", "FName: fname4", "Fval: fval4.fname4")

Hi, I am trying to find a simple way to numerically solve a system of two equations equal to zero with two unknowns (x_loc and y_loc). Here is a mock data set and below it, the equations I need to solve. theta<-c(180,135,90)/(2*pi) x<-c(0,0,15) y<-c(20,0,0) 0 = -sum((y_loc-y)*(sin(theta)*(x_loc-x)-cos(theta)*(y_loc-y))/(((x_loc-x)^2+(y_loc-y)^2)^0.5)^3) 0 =