similar to: logistic regression for a data set with perfect separation

Hi I'm sorry for my so basic question, but it's so basic that I can't find anwser anywhere... I have data like this: [1] 16.8000000 6.5333333 5.0666667 3.9333333 2.2000000 1.6666667 [7] 1.2000000 1.0666667 0.7333333 0.6666667 [1] 35.6113946 6.9576953 4.5271667 2.3744674 1.4735768 1.1751393 [7] 1.2649111 0.7988086 0.7037316 0.7237469 [1] 32.9333333 15.6666667

Hi, I have this data.frame with two variables in it, > z V1 V2 1 10 8 2 NA 18 3 9 7 4 3 NA 5 NA 10 6 11 12 7 13 9 8 12 11 and a vector of means, > means <- apply(z, 2, function (col) mean(na.omit(col))) > means V1 V2 9.666667 10.714286 My intention was substracting means from z, so instictively I tried > z-means V1 V2 1 0.3333333

Hello all! I have a problem with ggplot2 library. I want to do an heat map and the y variables are the year months. If I use the following code, he y values are in alphabetical order, but I want it in month order. The code is: library(reshape) library(ggplot2) library(scales) p <- ggplot(data.m, aes(variable, Month)) + geom_tile(aes(fill = value),

Hi all, I have an issue with the lm() function regarding the listing of the coefficients. My data are below, showing a list of hours (HR) relating to the time spent resting (R) by an individual animal. Simply i want to run a lm() to run in an anova() to see if there is a significant difference in resting between hours. HR R 1 2 0.6666667 2 2 0.4666667 3 2 0.8000000 4 2

Hello, I am having some trouble with the 'censoring-adjusted C-index' by Uno et al, in the package survAUC. The relevant function is UnoC. The question has to do with what happens when I specify a time point t for the upper limit of the time range under consideration (we want to avoid using the right-end tail of the KM curve). Copying from the example in the help file: TR <-

Dear R-users: I have been working in R for few years and never I'd been the problems like this. I'm running the last R version and when I try to save a graph in PS format I have the next error Error: Invalid font type In addition: Warning messages: 1: font family not found in PostScript font database 2: font family not found in PostScript font database I thinks it's curious

Hi This is my first question.. (sorry for my poor english.. ) Is this right way to ask something about installing xen, build kernel ...etc using mailing list?? I want know "How or Where can I ask something & How can I get your smart answers" I have followed this tutorial " http://wiki.xen.org/xenwiki/Fedora13Xen4Tutorial" exactly on my Core 2 Dou. after config the

Hi This is my first question.. (sorry for my poor english.. ) Is this right way to ask something about installing xen, build kernel ...etc using mailing list?? I want know "How or Where can I ask something & How can I get your smart answers" I have followed this tutorial " http://wiki.xen.org/xenwiki/Fedora13Xen4Tutorial" exactly on my Core 2 Dou. after config the

Hello, I have a matrix (class matrix) composed of GridCell (row and column). The matrix value is the beta diversity index value between two grids. Now I would like to get the average value of each GridCell. Please kindly advise how to make the calculation. Thank you. Elaine The matrix looks like (cited from Michael Friendly) I would like to get the average value of each color. Obs

Hi As you may''ve seen from my last posts I''m struggling to get the setup on my subject line to work without flaws. I''ve tried compiling the Kernel from Source, I''ve also tried Fedora 14, which has a flaky pv_ops Dom0 kernel which DOESN''T work at all, but none of these seem to work. I''m a mid-time XEN user, I just love it and want to stay

Hi As you may''ve seen from my last posts I''m struggling to get the setup on my subject line to work without flaws. I''ve tried compiling the Kernel from Source, I''ve also tried Fedora 14, which has a flaky pv_ops Dom0 kernel which DOESN''T work at all, but none of these seem to work. I''m a mid-time XEN user, I just love it and want to stay

Hi I need to define the following c("one group" = class.weight[2], "other group" = class.weight[1]) #class.weight = c(1,2) but I don't like the hard-coded way and would like to use my.group <- array(c("one group", "other group")) but now c(my.group[1] = class.weight[2], my.group[2] = class.weight[1]) gives an error how can I solve this

again, under another subject: sorry, maybe an all too trivial question. But we have power data from J frequency spectra and to have the same range for the data of all our subjects, we just transformed them into % values, pseudo-code: power[i,j]=power[i,j]/sum(power[i,1:J]) of course, now we have a perfect linear relationship in our x design-matrix, since all power-values for each subject sum up

Dear R-users How can I get the eigenvalues out of an lda analysis? thanks a lot christoph -- Christoph Lehmann <christoph.lehmann at gmx.ch>

Hi I need to overlay two pixmaps (library (pixmap)). One, a pixmapGrey, is the basis, and on this I need to overlay a pixmapIndexed, BUT: the pixmapIndexed has set only some of its "pixels" to an indexed color, many of its pixels should not cover the basis pixmapGrey pixel, means, for this "in pixmapIndexed not defined pixels" it should be transparent. What would you

Hi there, I've created a litle video explaining on how to setup some group policy to redirect folders using Active Directory Tools and Group Policy Management tools, it works and I've tested it so i decided to include it on my small youtube Linux collection. Hope you enjoy it and like; might be it can make it to the wiki video library. Thanks for Samba 4, it really _IS_ the thing.

Hello list, i've got some problems getting a E1 line with a E100P up and running (germany). # cat /proc/zaptel/1 Span 1: WCT1/0 "Digium Wildcard E100P E1/PRA Card 0" HDB3/CCS/CRC4 YELLOW RED "YELLOW RED" sounds not so good. When launching asterisk, enabling pri debug on that span, i see outgoing attempts: ;---------- snip -------- > [00 01 7f ] > Unnumbered

am I right in the assumption, that for calculation of the condition numbers I have to use the correlation matrix of X, and not t(x) %*% x? > e <- eigen(t(x) %*% x) better (x must not have a first column of ones): > e <- eigen(cor(x)) > e$val [1] 6.6653e+07 2.0907e+05 1.0536e+05 1.8040e+04 2.4557e+01 2.0151e+00 > sqrt(e$val[1]/e$val) [1] 1.000 17.855 25.153 60.785 1647.478

Hi Since the 'equal slope' assumption doesn't hold in my data I cannot use a proportional odds model ('Design' library, together with 'Hmisc'). I would like to try therefore a partial proportional odds model Please, could anybody tell me, where to find the code and how to specify such a model ..or any potential alternatives many thanks for your kind help christoph

Hi It seems, I just miss something. I defined treshold <- function(pred) { if (pred < 0.5) pred <- 0 else pred <- 1 return(pred) } and want to use apply it on a vector sapply(mylist[,,3],threshold) but I get: Error in match.fun(FUN) : Object "threshold" not found thanks for help cheers chris -- Christoph Lehmann <christoph.lehmann at gmx.ch>