search for: stephanis

Dear R-users, Is there a function in R that classifies data in more than 2 groups using logistic regression/classification? I want to compare the c-indices of earlier research (lrm, binary response variables) with new c-indices obtained from 'multiple' (more response variables) logistic regression. Best regards, Stephanie Delalieux Department Biosystems M?-BIORES Group of Geomatics

I (sometimes) get this error message: WARNING[17191]: app_voicemail.c:2411 leave_voicemail: No entry in voicemail config file for 'irstname.lastname' I can see the value of the argument is "firstname.lastname" when this line executes in the std-exten macro: exten => s-NOANSWER,1,Voicemail(u${ARG1}) ; If unavailable... But the error message drops the first character. It

Hello, I have a general formatting question. I have two columns of data: ColA <- c("m", "m", "m", "m") ColB<- c("d","d","d","d") And I would like to reorder them into a new column that looks like this: ColC<-

When I use the command for PACF, lags 5, 10, 15, and 20 are labeled. I would like to label lag 1. I would greatly appreciate if someone could tell me how to do this. Below is the command that I am using: pacf(data$R1,main="Series R1 Residuals") [[alternative HTML version deleted]]

Hello List, My question is an elementary one. I have run a fuzzy kmeans cluster using FANNY to group freshwater fish assemblages. I then went in the field to validate that classification and have retrieved new assemblage data for a new suite of streams. Therefore I would like to use Predict to determine how well the original clustering fits the new data. However I have not figured out a

Full_Name: Stephanie MAHEVAS Version: 2.4.1 OS: Windows NT Submission from: (NULL) (134.246.55.50) the two following instructions provide a synthax error : if ( 5 > 4 ) cat("ok1") else cat("ok2") and if ( 5 > 4 ){ cat("ok1")} else cat("ok2") whereas these ones don't if ( 5 > 4 ) cat("ok1") else cat("ok2") and if (

Hi all, may i know is it possible to plot a graph by first letter? for example: Name: Age: Angel 20 Amelia 20 Bernard 19 Stephanie 20 Vanessa 22 Angeline 23 Camel 21 If I want to plot the name started with letter 'A' and their Angel,

Hello, I would like to place shaded boxes on different areas of a phylogenetic tree plot. Since I can not determine how to find axes on the phylogenetic tree plot I am not able to place the box over certain areas. Below is example code for the shaded box that I have tried to use, and the first four values specify the position. rect(110, 400, 135, 450, col="grey",

Hi Everyone - Here is a draft of the first of two "pages" that I will submit to Wikipedia. The first page is a general network virtualization page. From this page, I will contribute a subpage on the Crossbow projects. These two articles are the first of a number of online articles on Crossbow that I''ll be creating on behalf of the Crossbow iteam. The article uses

Hello, I would like to remove some files which have the extension .test for example (data1.test, data2.test ....). Is there another solution to remove them instead of doing it one by one ? Thanks for your help, St?phanie Langevin -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Is there a way to set "na.rm=TRUE" as a default, so that this does not have to be re-specified for all of the functions requiring this option? Thanks in advance, Stephanie Broyles sbroyl@lsuhsc.edu [[alternative HTML version deleted]]

Hi, I have been trying to use the function cochran.test from the Outliers package to test for homogeneity of variance. This works well except when I use transformed data. Would anyone have an idea why it doesn't work and how I could do the cochran test on transformed data? Thanks, Stephanie >library(outliers) > set.seed(1234) > x=rnorm(100) >

Hello, I am using tapply to pull out data by the day of week and then perform functions (e.g. mean). I would like to have the number of values used for the calcuation for the functions, sorted by each day of week. A number of entries in any given column are NAs. I have tried the following code and simple variants with no luck. for (i in 1:length(a[1,])){ x<-tapply(a[,i],a[,1],mean,

Hi, We have 7 Dell 2850 servers with dual xeon 3 gig processors running the APF firewall version 0.9.6 http://rfxnetworks.com/apf.php They run fine for a day or two, then suddenly lock out all incoming connections, other than the backend IP, sometimes restarting the firewall resolves this, but occasionally we may have to leave it 10 mins or so before restarting where it will actually

In my own package, I want to use the default S3 method of the generic function lrtest() from the lmtest package. Since I need only one function from lmtest, I tried to use importFrom in my NAMESPACE: importFrom(lmtest, lrtest) However, this fails R CMD check in the examples: Error in UseMethod("lrtest") : no applicable method for 'lrtest' applied to an object of class

I have a loop that I would like to use to extract the "stddev" for each itteration so I can average the "stddev" for all the runs. It would be helpful to know how to extract the "stddev" for each run from the VarCorr. Thanks MCruns<-1000 sighatlvec<-rep(NA,MCruns) sighatbvec<-rep(NA,MCruns) sighatevec<-rep(NA,MCruns) for(mc in 1:MCruns) {

I use a package to contain simple functions that can be handled by unit tests for correctness and more complex functions that combine the simple functions with business logic. Where there are proposals to change either the simple functions or the business logic, a sample needs to be run before the change and then after it to understand the impact of the change. I do this currently by 1. Using

Hi all, I''m running puppet agent as a onetime process. I would expect, if it can''t create the pid file, that it would exit non-zero. However, it seems to exit zero. Is there any way to change this? I''d like to know if it couldn''t run, and it seems that this is an error. [root@db1r1 ~]# puppet agent --onetime --detailed-exitcodes [root@db1r1 ~]# puppet

Hello all, I'm trying to plot the effects of variables estimated by a regression model fit individually, and for categorical predictors, the independent variable shows up on the y-axis, with the dependent variable on the x-axis. Is there a way to prevent this reversal? Sample code with dummy data: # make dummy data set.seed(1) x1 <- runif(200) x2 <- sample(c(1,2),200, TRUE) x3 <-

Dear all, I am having trouble creating a "jagged array" in R. I have found no commentary on how to do it in the help files, but maybe I am misunderstanding the purpose of some of the array functions (e.g. tapply). I am using a longitudinal dataset where each individual has a different number of time points (observations). I need to create correlation structures for each