similar to: scan problems -- what can "what" be?

Hi all, Using: R version 2.8.1 Patched (2009-03-07 r48068) on OSX (10.5.6) with survival version: Version: 2.35-3 Date: 2009-02-10 I get the following using the first example in ?summary.survfit: > summary( survfit( Surv(futime, fustat)~1, data=ovarian)) Call: survfit(formula = Surv(futime, fustat) ~ 1, data = ovarian) time n.risk n.event survival

I dont get it: > is.vector(c(mydf[1])) [1] TRUE > unique(c(mydf[1])) Error in unique(c(mydf[1])) : unique() applies only to vectors > Is it a vector or not? This stuff is driving me nuts. I'm simply trying to convince R that my grouping vector is actually a vector so that unique will work. Its just a vector of numbers, so why shouldnt it work? --

I am trying to get summary statistics from WinBUGS/JAGS output in the form of mcmc.list objects, using the summary() function. However, I get odd warning messages: Warning messages: 1: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : algorithm did not converge 2: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : algorithm did

Can one perform a QDA in R? I do not see it anywhere within the mda package. Any pointers here would be appreciated. Thanks, cjf -- ~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~* Christopher J. Fonnesbeck Ph.D. Student Georgia Cooperative Wildlife Unit University of Georgia Athens, GA 30602 Email: cjf at fonnesbeck.net Yahoo: fonnesbeck_chris ~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*

Dear list, I have the following matrix. How can I make the following plot? 1. The x-axis has index 1:7, and the first column is plotted against index 1, second against 2, and so on. 2. I want the points from the left upper conner including the antidiagonal to be plotted with col=2, and the rest with col=3 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.589 0.857 0.923 0.944 0.954 0.963

I'm hoping someone with some experience generating tables in Sweave will be able to solve this problem for me. I'm experiencing some inconsistency with the way floating point numbers are displayed -- for some tables, they will be in standard notation, e.g. 11±26 while in others, it uses exponential notation, e.g. 1.1e+01±2.6e+01 I am using code such as the following; note I have

Dear R-helpers, I'm stuck with a little problem that surely has an easy solution but I can't think of a way to solve it. I'd really appreciate any help you can offer me! I'll provide a small example. Given a dataframe data.txt that looks like this: ID freq Var Var_mean Ratio_mean Var_median Ratio_median Var_sum Ratio_min Var_max Ratio_max Var_min

I was wondering if there is a way of predicting factor scores of new data for factor analysis in R (similar to "predict" in S-plus). So far I have not been able to find it, nor found reference to it. Thanks, Chris Fonnesbeck -- ~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~* Christopher J. Fonnesbeck Ph.D. Student Georgia Cooperative Wildlife Unit University of Georgia Athens, GA 30602

I am trying to build 1.6.2 from source RPM's on RedHat 9, but encounter a failure when running the base tests. An inspection of base-Ex.Rout.fail reveals the following: > ##___ Examples ___: > > var(1:10)# 9.166667 [1] 9.166667 > > var(1:5,1:5)# 2.5 [1] 2.5 > > ## Two simple vectors > cor(1:10,2:11)# == 1 [1] 1 > > > stopifnot( is.na(var(1)), +

Hi, I am trying to make a multi-class classification tree by using rpart. I used MASS package'd data: fgl to test and it works well. However, when I used my small-sampled data as below, the program seems to take forever. I am not sure if it is due to slowness or there is something wrong with my codes or data manipulation. Please be advised ! The data is described as the output from str()

Hi, When I run the survfit function, I want to get the restricted mean value and the standard error also. I found out using the "print" function to do so, as shown below, print(km.fit,print.rmean=TRUE) Call: survfit(formula = Surv(diff, status) ~ 1, type = "kaplan-meier") records n.max n.start events *rmean *se(rmean) median 200.000

OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) install.packages("systemfit")

... and the mailing list is picky about attachments... whatever you attached did not conform to the stringent requirements mentioned in the Posting Guide. Pasting the code right into the email is usually safest, though you DO have to post using plain text (as the Posting Guide indicates) or your code may get mangled by the automatic html format removal. On May 15, 2018 7:04:31 AM PDT, Bert Gunter

Sadly you failed to set your email program to send plain text and the data is corrupted at my end. I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis. I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need

Unless there is good reason not to, always cc the list -- there are lots of smarter folks than I on it who can help. I may or may not have time to look at this. Hopefully someone else will. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip

I noticed in the collection of contributed packages for the OSX release of R, there is no RMySQL package. Does anyone know if one exists, or otherwise, how to connect to relational databases from R on Mac? Any help mode appreciated. cjf

Dear list, i'm a new R user, so I apologize if the topic is already being addressed by some other user. I'm trying to determine if the reproductive success of a species of bird is related to a list of covariates. These are the covariates: ? elev: elevation of nest (meters) ? seadist: distance from the sea (meters) ? meanterranova: records of temperature ? minpengS1: records

Hi I have the following array: test <- array(c(1:16), dim = c(3,4,3)) test ## I call some enries using an index array test.ind <- array(rbind(c(1,2,1), c(3,3,2)), dim = c(2,3)) test[test.ind] ## suppose I want all values in the 2nd row and 4th col over ## all three 3rd dimensions test[2,4,] how to specify a test.ind array with the last index left with ',' i.e test.ind should be

The previously posted benchnarks for gcc vs llvm-gfortran had one mistake. I was actually had the gfortran for 4.3.1 installed instead of that from gcc 4.2.4. Below are the polyhedron benchmark results for all three compilers... gfortran 4.2.4 Benchmark Compile Executable Ave Run Number Estim Name (secs) (bytes) (secs) Repeats Err % --------- -------

Hello Where could I find examples on how to work with the time index in a timeseries or zoo series? Let say I've got this series DATA 1990-01-01 10:00:00 0.900 1990-01-01 10:01:00 0.910 1990-01-01 10:03:00 0.905 1990-01-01 10:04:00 0.905 1990-01-01 10:05:00 0.890 ....................... 2000-12-31 20:00:00 0.992 How do I make simple calculations such as ... ? Calculate the