similar to: P-values from survreg (survival package) using a clusterterm

Hi all. After upgrading to R2.0.1 i get > "??" [1] "??" > "??" [1] "\370" > "??" [1] "??" Whereas under R1.9.1 i get > "??" [1] "??" > "??" [1] "??" > "??" [1] "??" Any hints apreciated. Steen Steen Ladelund, statistician +4543233275

Hi all. I think it would be nice to be able to combine levels of a factor on creation a la x <- rep(0:5,5) y <- factor(x,levels=0:5,labels=c('1','1',2:5)) ## (1) y [1] 1 1 2 3 4 5 1 1 2 3 4 5 1 1 2 3 4 5 1 1 2 3 4 5 1 1 2 3 4 5 Levels: 1 1 2 3 4 5 I thougt this would (should?) create a factor with 5 levels, ie combining 0 and 1 in x into one level in y. I find it

Hi helpeRs. When I do two-ways tables dimnames in columns are almost left-adjusted: <- factor(rbinom(30,1,.3),labels=c("first","second")) > b <- a > table(a,b) b a first second first 21 0 second 0 9 If I do a oneway table however, dimnames are right-adjusted: table(a) a first second 21 9 Actually I must admit I

Hi all. To automate the production of a document I want to have something like: test.tex : test.Rnw "Sweave("test.Rnw") "| R --no-restore --no-save in a makefile. This however gives an error as R states: Sweave("test.Rnw") : not found. I use R-1.9.0pat (040424 I think) on Win XP. Any suggestions, hints or plain help appreciated. Steen Ladelund

Dear R users, I have a response variable in a csv file called "y" and a matrix of predictor variables in a csv file called "mat". I have used the function "nls" I have specified the nonlinear relation between these variable.The code I have witten is called Rprog which begins with the phrase: L.minor.m1<-nls(Y~a ....etc.. The program when I execute the program, I

Hi everyone, I have dataset which I make a sample of it couple of times and each time I get the mean and standard deviation of each row for each sample. I have a function for that, which takes the name of the file and number of times to sample and then returns the mean and standard deviation for each row in each sample. Sample=function(name, n){

Hello, I'm new here, but will try to be as specific and complete as possible. I'm trying to use “lm“ to first estimate parameter values from a set of calibration measurements, and then later to use those estimates to calculate another set of values with “predict.lm”. First I have a calibration dataset of absorbance values measured from standard solutions with known concentration of

Hi there, I have a the following function which takes a name(fname), total number of rows (Total), number of times to make a sample(nSample), number of records/rows to be included in each sample(nPatient). then it generates required samples and calculates the mean and standard deviation of the columns for each sample and finally in calculates the grand mean and standard deviation. The problem

Hello, I am using package MuMIn to calculate AIC for a full model with 10 explanatory variables. Thanks in advance in sharing your experience. Q1 In the AIC list of all models, each model is differentiated by model number. Please kindly advise if it is possible to find the corresponding explanatory variable(s) for the model number. Q2 error message I tried to display sub-model with only

Hello, I have a simple problem (checked the archives and the appropriate help pages, to no avail). I want to creat a vector that is length(2000). The vector is to consist of two strings( "std" and "dev") with 1760 "std"s and 240 "dev"s. Furthermore, for each element of the vector, i want the selection of one of the strings to be approx. random. The

Full_Name: Grant Izmirlian Version: 2.2.0 OS: SuSe Linux version 9.2 Submission from: (NULL) (156.40.34.177) Sorry about the last submission, my bug-fix had an error in it because ifelse doesn't vectorize. I'll repost with the correct bug-fix. ------------------------------------------------------------------------------- The option exists to include all parameters, including the

Hi I am trying to print out means, STDs and histograms under two sets of factors. I can manage it for one set - see below but not for two sets. That is, I want ot print out the mean STD and Histogram for each ITEM code within each DELIVERABLE code. In addition I can only get to view the histogram for the last item. How do you get R to stop overriding the histogram for eg level 1 for factor 1

Hi, I need to produce a standard report for several variables in Sweave and thus would need the possibility to define a TeX-command which includes R-code like \newcommand{\meansd}[1]{The mean is \Sexpr{mean(#1)} and the standard deviation is \Sexpr{sd(#1)} . } and then just write \meansd{age} in the latex code to get the whole sentence. The above does not work, since

under Ubuntu 9.04 R seems to be very slow at plotting. the example below illustrates with a plot of error bars of sample means where i watch as each error bar is plotted one at a time. very annoying and pain in the neck when running Sweave repeatedly. running R 2.9.1 under Windoze on the same machine the error bars plotted in the code below appear instantaneously. has anyone else noticed this

Hi All- I have SuSE 8.2 installed on two different Intel machines, with (as far as I can determine) all required packages for doing what I'm trying to do here. Both machines suffer the same problems. I didn't find any SuSE 8.2 rpms, so I'm trying to build Samba 3.0 (stable) (24th Sep, 2003) on either of these machines but am having problems with it. As root, I did a

Hi I noticed a small inconsistency when using sum() vs cumsum() I have a char-based series > tryjpy$long [1] "0.0022" "-0.0002" "-0.0149" "-0.0023" "-0.0342" "-0.0245" "-0.0022" [8] "0.0003" "-0.0001" "-0.0004" "-0.0036" "-0.001" "-0.0011"

Consider this situation: > x1 <- runif(100); x2 <- runif(100); y <- 2 + 3*x1 - 4*x2 + rnorm(100) > m1 <- summary(lm(y ~ x1)) > m2 <- summary(lm(y ~ x2)) > m3 <- summary(lm(y ~ x1 + x2)) Now you have estimated 3 different "competing" models, and suppose you want to present the set of models in one table. xtable(m1) is cool, but doing that thrice would give

>>>>> Tomas Kalibera >>>>> on Tue, 25 Aug 2020 09:29:05 +0200 writes: > On 8/23/20 5:02 PM, Rory Winston wrote: >> Hi >> >> I noticed a small inconsistency when using sum() vs cumsum() >> >> I have a char-based series >> >> > tryjpy$long >> >> [1]

Hi all, I'm going nuts. Ken Weiss and Steve Holstead have been very nice to me and supplied me with info on how to set my configure files and such. It helped me a lot but I still get stuck a bit later in the process. Namely during pass_check.c. I'm on a AIX 4.1.5 box. Using the export version of AFS 3.4. Help is appreciated. I wouldn't ask you all this if I hadn't already spent

My data frame looks like: SightingID PA1 PA2 PlotID InOverlap Area1 2001 1 -99 392 Y 0.22 2002 1 -99 388 Y 0.253 2008 1 NA 104 N 0.344 2010 1 NA 71 N 0.185 2012 1 NA 61 N 0.166 2013 1 NA 61 N 0.227 2014 1 NA 62