similar to: How to change the lengend in a DA figure

I must have missed something simple, but still, i don't know what. I obtained my basic data as follows: x <- sprintf("SELECT m_id,sale_date,YEAR(sale_date) AS sale_year,WEEK(sale_date) AS sale_week,return_type,0.0001 + DATEDIFF(return_date,sale_date) AS elapsed_time FROM `merchants2`.`risk_input` WHERE DATEDIFF(return_date,sale_date) IS NOT NULL") moreinfo <- dbGetQuery(con,

I have some data that was created for import into a Lotus 1-2-3 spreadsheet and on of the columns is time. The time is akin to Julian were the value 1 is mapped "01-Jan-00 12:00:00 AM" in Lotus 1-2-3. Is there a function in an R package that can convert this numeric vector to a POSIXct vector? With best wishes and kind regards I am Sincerely, Corey A. Moffet Instructor Department

I am trying to reproduce an analysis done in SAS from which Type III Sums of squares are used. How can you get Type III sums of squares in an aov similar to what SAS provides? With best wishes and kind regards I am Sincerely, Corey A. Moffet Instructor Department of Range, Wildlife, and Fisheries Management Mail Stop 2125 Texas Tech University Lubbock, Texas 79409-2125 (806) 742-2843 Voice

I am feeling rather dumb right now. I created what I thought was a data.frame as follows: aaa <- lapply(split(moreinfo,list(moreinfo$m_id),drop = TRUE), fun_m_id) m_id_default_res <- do.call(rbind, aaa) print("==========================================") m_id_default_res print("==========================================") ndf <- m_id_default_res[, c('mid',

Dear maintainer, I used bioconductor package, pdInfoBuilder, to make a microarray platform annotation package. I named this package, pd.pdinfo.gpl11164.ndf.txt. This "self-made package" is to be used with oligo package. Prior to using oligo package, I need to install the annotation package, pd.pdinfo.gpl11164.ndf.txt, and cannot just copy it to the library tree. I tried using the

The second beginner question. I want to create a new dataframe, where each element of the original dataframe is translated to 1 if it was "+", to 0 if it was "-" to -1 otherwise. I could do with: Lines <- "a b c d + - + + + + + - + 1 - '+ ' - + + + + N - +" DF <-

Hello! I have a function: zywnoscCalosc<- function( jedzenie, n1, n2, n3, n4, d1, d2, d3, d4 ) { ndf <- data.frame(nn1=n1,nn2=n2,nn3=n3,nn4=n4) ddf <- data.frame(dd1=d1,dd2=d2,dd3=d3,dd4=d4) for (i in 1:length(n1)){ wekt_n = ndf[i,] wekt_n_ok = wekt_n[!is.na(wekt_n)] dl_n = length(wekt_n_ok) wynik = (1*wekt_n_ok)/(1*dl_n) } } and I get an error like this: Error in 1 * wekt_n_ok :

Hello, I want to know what do you think about my function. I know that it isn't briliant :/ but what do you think? What I should do that my function will be better? (now is very slow and not ideal, sometimes I also get a mistake!) ########## My function ############################################# dzieci<-transform(dzieci, zywnosc=0) zywnoscCalosc<- function( jedzenie, sklep, n1, n2,

I have a function for plotting soil texture that I am having a problem with. The function is: plot.psd <- function(sand, clay, ...) { conv.ter <- function(x, y) { x.con <- 100 - x - y*cos(1/3*pi) y.con <- y*sin(1/3*pi) data.frame(x=x.con, y=y.con) } plot(conv.ter(sand, clay), xlim = c(0,100), ylim = c(0,100*sin(1/3*pi)), axes = FALSE, xlab="",

Hi, I have run into another problem using offsets, this time with the predict function, where there seems to be a contradiction again between the behavior and the help page. On the man page for predict.lm, it says Offsets specified by offset in the fit by lm will not be included in predictions, whereas those specified by an offset term in the formula will be. While it indicates nothings about

hi Rusers I am trying to use symbos in gtools package symbols(data1,data3,circle=data1/data3,inches=0.1,bg="lightgreen") Now I want to give a lengend to tell the reader the meaning or magnitude of these circle. How can I add these information in symbols plot just like legend in plot ? thank you . -- TANG Jie Email: totangjie@gmail.com Tel: 0086-2154896104 Shanghai Typhoon

Full_Name: Gerard Torrent Version: R version 2.8.0 (2008-10-20) OS: Linux 2.6.27.5-41.fc9.x86_64 #1 SMP Submission from: (NULL) (85.52.227.233) In some cases qt complains about NaNs and don't gives the correct result: > qt(0.1,, 0.1) [1] NaN Warning message: In qt(p, df, lower.tail, log.p) : NaNs produced But the result can be found: > pt(-1.60443e+06, 0.1) [1] 0.09999997 If I

Tim Thompson Commercial Sales Engineer http://www.amatechtel.com (806) 722-2227 Would anyone have an idea on how I would be able to take the mic in on the computer and put it as the "talking party" for a conference room. I would then be able to set up a "listen only" profile for others to get in on. Reason for doing this is for 'shut-in's' for my Church.

In one particular situation predict.coxph gives an error message. Namely: stratified data, predict='expected', new data, se=TRUE. I think I found the error but I'll leave that to you to decide. Thanks, Chris ######## CODE library(survival) set.seed(20121221) nn <- 10 # sample size in each group lambda0 <- 0.1 # event rate in group 0 lambda1 <- 0.2 # event rate in group 1

Hello, I am a new user of R and Linux (debian squeeze), and I am trying to install the ncdf library. First I don't know wich is the best to use (ncdf or ncdf4), But not matter the one I tried to install, it always failled. I try to find a solution to this common problem on the web, but I never find the good way to solve it. I hope you could help me I try many things during two days and now I

Hello compadRes, I'm developing a script that selects "cells" over a certain metabolic rate to kill them. A rate between 9 and 12 means that the cells are candidates for death. I'll show you what I mean: # a would be a vector of cell metabolic rates. a<-c(8, 7, 9, 6, 10, 11, 4, 5, 6) #now identify which cells will be candidates for death, namely those cells with metabolic

Is there a function/package that will allow you to test the hypothesis beta1 = x in a simple linear regression, where x is a constant? The AOV tests if beta1 in different from 0 but what if I want to know if beta1 one is different from 1? With best wishes and kind regards I am Sincerely, Corey A. Moffet Instructor Department of Range, Wildlife, and Fisheries Management Mail Stop 2125 Texas

Hello list, I'm new to this list, so please forgive my ignorance. I have searched R-help for some hints into what might be my problem, but I truly have no idea where to go from here. I have an object of approximately 15,000 rows and 2 columns. There are many duplicates in the first column, all with different corresponding values in the second column. For example (2 is duplicated):

How can I extract the adjusted R^2 value from an ols object (using rms package)? library(rms) x <- rnorm(10) y <- x + rnorm(10) ols1 <- ols(y ~ x) Typing "ols1" displays adjusted R^2 among other things, but how can I assign it to a variable? I tried str(ols1) but couldn't see where to go from there. Thanks, Mark Seeto

v2: use the same object for private data as with the drm debugfs functions Signed-off-by: Karol Herbst <nouveau at karolherbst.de> --- drm/nouveau/nouveau_debugfs.c | 54 +++++++++++++++++++++++++++++++++++++++++-- 1 file changed, 52 insertions(+), 2 deletions(-) diff --git a/drm/nouveau/nouveau_debugfs.c b/drm/nouveau/nouveau_debugfs.c index 5392e07..fe53743 100644 ---