similar to: Stop Warnings for Invalid Factor Level, NAs generated?

I was wondering if someone familiar with survival analysis can help me with the following. I would like to fit a Weibull curve, that may be dependent on a covariate, my dataframe "labdata" that has the fields "cov", "time", and "censor". Do I do the following? wieb<-survreg(Surv(labdata$time, labadata$censor)~labdata$cov,

Hello, I'd like to generate a filled step-function in R as you can see in the attachement (image is generated by grapher, but I need a freely available alternative). The problem is, that I don't know, how to fill it up. My code right now looks like this: pH <- c(0,6.1,6.1,6.3,6.6,7.3,0) Tiefe <- c(0,0,10,40,80,120,10) plot(pH,Tiefe,type="s",main="Tiefenfunktion

What is this error message about and how do I troubleshoot it? > sqlTables(channel) Error in .Call("RODBCFetchRows", attr(channel, "handle_ptr"), max, buffsize, : negative length vectors are not allowed The channel was created as such channel <- odbcConnect("Labdata") I think this is telling me that my odbcConnect command worked >

Greetings! As part of my research project I am using R to study temperature data collected by a network. Each node (observation point) records temperature of its surroundings throughout the day and generates a dataset. Using the recorded datasets for the past 7 days I need to build a prediction model for each node that would enable it to check the observed data against the predicted data. How can

it''s been a nightmare trying to set it up. I keep getting a an Application Error message when I try to navigate to a url which should be taken care of by my newly created controllers. For instance, I wanted to test RoR out so I created a MyTest controller which should allow me to navigate something like: http://localhost/rubytest/MyTest but it does not. I have no problem getting

I''d like to delete some records before saving them. I ran this experiment: before_save :prune_lab_datas def prune_lab_datas lab_datas.each do |lab_data| lab_data.destroy end raise lab_datas.to_yaml end But nothing seems to be deleted in the hash: --- - &id001 !ruby/object:LabData attributes: created_at: 2010-03-21 22:16:06 unit_id: "3"

I am having trouble using get.bnddata() in the RArcInfo package. I installed the version 0.2 binary of the package for Windows and am using R 1.4.1 on Windows NT 4.0. I have had no problems importing and plotting arc data with get.arcdata() and plotarc(). plotpal() also works fine. I've been able to use get.arcdata(), get.labdata(), get.paldata(), get.tablenames(), get.tablefields(),

I find myself wanting to plot three graphs side by side 'as if' they were panels -- that is, with the same y-axis limits, no space between the graphs, and precise vertical alignment of the plot areas. However, I don't want strip titles; I want each graph to have its own x-axis label, on the bottom of the plot. The best way I have so far found to do this is to fake up a data frame that

I don't know of any R package that can match all the functionality of LEM eg fitting equality constraints to model parameters a la LISREL. WRT dumping tables, I would have thought that as.data.frame.table does pretty much what you want, [not tested] newtab <- as.data.frame(table(a,b,c)) cat("dim\n") for(i in seq(1, ncol(newtab)-1) { cat(nlevels(newtab[,1]," ") }

I have a number of plots that I'm trying to do using the lattice package, unfortunately, they involve Posix values. A small sample of the data is as follows: Time TC.1 <Snipped> 1 2003-04-24 13:29:15 26.61 2 2003-04-24 13:29:30 26.48 3 2003-04-24 13:30:00 26.53 4 2003-04-24 13:30:30 27.85 <Snipped> where str(data) `data.frame': 415 obs. of

Hello: I am using rpart and would like more background on how the splits are made and how to interpret results - also how to properly use text(.rpart). I have looked through Venables and Ripley and through the rpart help and still have some questions. If there is a source (say, Breiman et al) on decision trees that would clear this all up, please let me know. The questions below pertain to a

Hello list. I'm trying to rebuild the 2.6.9.67.0.20.EL kernel, but it fails even without modifications. How did I try it? Created a (non-root) build environment (not a mock ) Installed the kernel.scr.rpm and did a rpmbuild -ba --target=`uname -m` kernel-2.6.spec 2> prep-err.log | tee prep-out.log The build failed at the end: Processing files: kernel-xenU-devel-2.6.9-67.0.20.EL Checking