Re: [R] How to obtain the estimate of baseline survival function?

See ?basehaz in 'survival' package. -Christos -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Zheng Yuan Sent: Friday, November 03, 2006 6:20 PM To: r-help@stat.math.ethz.ch Subject: [R] How to obtain the estimate of baseline survival function? Hi, If

Re: [R] multiple plots in the same graph

See ?points and ?lines for adding new groups of points or lines to an existing graph created by plot. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Li Zhang Sent: Friday, November 03, 2006 9:55 PM To: R-help@stat.math.ethz.ch Subject: [R] multiple plots

Re: [R] nlme Error: Subscript out of bounds

I can't begin to guess. If your example were self contained (and preferably simple), it would be easier to diagnose. When I have problems like this, I often try to find a simple example that produced the same error message. That process often leads me to a solution. If it doesn'

[R] multiple plots in the same graph

I'd like to plot y vs x according to the third variable "group" which has three levels. I am wondering how can I put the three plots in one graph? Thank you (http://groups.yahoo.com) ___

Re: [R] as.zoo behavior (

Mark: > temp<-ewma(actualdiff,.05), ok, because ewma() is a wrapper for filter() this returns an object of class "ts". Therefore you want to add the index again, but: > rollmeandifflogbidask<-as.zoo(temp,index(actualdiff)) is not the right thing to do. as.zoo() is a coercion generic and its "ts

[R] Ridge for logistic regression

Dear all experts, Does anyone know if there is a R function which can perform Ridge regression for logistic models? Thank in advance. Ying Yuan __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the po

Re: [R] Chronbach's Alpha

On 11/04/06 02:25, Ralf Finne wrote: > Dear all experts, > > Can anyone help me to find information on > how to calculate the Chronbach's Alpha Coefficient > (sometimes called Crombach or even just Alpha coefficient) And even more often called Cronbach. It is in the psy package. And there is s

Re: [R] (no subject)

On 03-Nov-06 Martin Maechler wrote: > [...] > But I want to make another point: > > For about a year now, for "serious" data-analysis using > data frames, I've been advocating to use the slightly more > clumsy but much more error-prone ``column indexing by name'' > instead of the quick-and-dirty "

[R] Chronbach's Alpha

Dear all experts, Can anyone help me to find information on how to calculate the Chronbach's Alpha Coefficient (sometimes called Crombach or even just Alpha coefficient) Thanks in advance Ralf Finne SYH University of Applied Science Wolffskavägen 33 P.O.Box 6 FIN-65201 VASA, Finland _

[R] plotting residuals

Does anyone know how to obtain a plot of residuals by predicted values for a main-effects aov? I want to check that the residuals are distributed equally across treatment means. Thanks, Jeff [[alternative HTML version deleted]] __

Re: [R] How to convert an existing package?

On Fri, Nov 03, 2006 at 07:59:37PM +0100, rasti matus wrote: > Dears, > > how could I convert an existing package to be able to modify it? > > Let s say, if I want to change the way how the results of the package are > display (to see the oroginal code of the code). Hi, Could you be more specif

Re: [R] question of one matrix column being treated as list rather than vector

On Fri, Nov 03, 2006 at 02:47:37PM -0800, Waverley wrote: Hi, You probably created that matrix from a list. Eg > a <- matrix(list(1,2,3,4),2,2) > a [,1] [,2] [1,] 13 [2,] 24 > a[,2] [[1]] [1] 3 [[2]] [1] 4 > b <- matrix(1:4,2,2) > b[,2] [1] 3 4 Both a and b are matrices, b

[R] How to obtain the estimate of baseline survival function?

Hi, If I fit a Cox model using "coxph", is there a R function so that I could obtain the estimate of baseline survival function? Thank you. Zheng -- Zheng Yuan Ph.D student Department of Biostatistics University of Michigan Ann Arbor, MI 48109 __ R-h

[R] as.zoo behavior (

hi all : the code pasted below runs but then, a dput on rollmeandifflogbidask gives me what is below the code. the structure of rollmeandifflogbidask is a zoo object but with a "frequency" so it's not the same structure as the original actual diff and this really causes things to blow up in later

[R] question of one matrix column being treated as list rather than vector

Hi, I have one matrix derived from my function. It is as following: Browse[1]> each.uaa.data [,1] [,2] [1,] 33 152.3530 [2,] 33 264.7059 [3,] 51 583. [4,] 33 323. [5,] 50 292.549 [6,] 33 300 [7,] 56 104.1176 [8,] 52 246.6667 [9,] 53 242.5490 [10,] 33 2

[R] finding significant change points after fitting GAM using TPRS

Does anybody know if code or functions exist for determining significant change points after fitting GAM in 'mgcv' using thin plate regression spline and 'ts' bases? I would like to identify years where the curvature or change in trajectory of a longterm annual population curve is significant. T

Re: [R] one row matrix

On 11/3/06, Waverley <[EMAIL PROTECTED]> wrote: > Hi, > For example: > > > a = matrix (12, 3, 4) > > b = a[1,] > > b > [1] 12 12 12 12 > It doesn't have anything to do with one-row matrices, at least in your example, since that isn't what you made. Instead, you created a matrix with 3 rows, 4 col

Re: [R] one row matrix

On Fri, 2006-11-03 at 14:12 -0800, Waverley wrote: > Hi, > > when I assign a one row matrix to another variable, then somehow R > automatically convert that varaible into a list. A 'vector', not a 'list' > For example: > > > a = matrix (12, 3, 4) > > b = a[1,] > > b > [1] 12 12 12 12 > > Is t

[R] one row matrix

Hi, when I assign a one row matrix to another variable, then somehow R automatically convert that varaible into a list. For example: > a = matrix (12, 3, 4) > b = a[1,] > b [1] 12 12 12 12 Is there a way to enable R automatically make b as a one row matrix , rather than explicitly assign like

Re: [R] Import problem with S-Plus 7.0 dataset

Thanks to Brian Ripley and Richard Heiberger for the solution to my problem. I will spell it out below in steps so that (hopefully) nobody else will need to ask the same question in the future (if they search the R help archive). HOW TO GET S-PLUS 7.0 DATA into R

[R] (no subject)

From: Martin Maechler <[EMAIL PROTECTED]> MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Message-ID: <[EMAIL PROTECTED]> Date: Fri, 3 Nov 2006 22:26:20 +0100 To: Barry Rowlingson <[EMAIL PROTECTED]> Cc: Sarah Goslee <[EMAIL PROTECTED]>, r-help@stat.

Re: [R] help went away from top emacs window :-(

Thanks for looking into this. On 11/3/2006 9:53 AM, Richard M. Heiberger wrote: > I located what might be the problem. There is a different misbehavior > with Gnu emacs. I don't have Xemacs. I have an idea how to fix it but > haven't worked out the details. > > > ess-inf.el in the lines > (de

Re: [R] help went away from top emacs window :-(

I concur: Thanks very much for looking into this. Spencer Graves Duncan Murdoch wrote: > Thanks for looking into this. > > On 11/3/2006 9:53 AM, Richard M. Heiberger wrote: >> I located what might be the problem. There is a different misbehavior >> with Gnu emacs. I don't have Xemacs. I have

Re: [R] Import problem with S-Plus 7.0 dataset

It looks like it works. The result you printed is the name of the file that data.store read. The name of the variable is the same as the name you called it in S-Plus. type data[1:5,] and you should see your data.frame. __ R-help@stat.math.ethz.ch m

Re: [R] Microarray analysis

Friederike - The program is doing exactly what you asked it to do, but please notice that varprobe<-c(var(schizo.matrix[x,])) creates a vector containing all the elements of the variance-covariance matrix of schizo.matrix[x,], which is a 2x150 matrix, giving a variance-covariance matrix

Re: [R] Import problem with S-Plus 7.0 dataset

As recommended, I have tried the following code to solve my problem importing data to R from S-Plus 7.0. Unfortunately, I have not had success. In S-Plus: > data.dump('data', file='C:\\temp\\ddump.sdd', oldStyle=T) This resulted in the production of a file called "ddump.sdd" that I can import i

[R] Microarray analysis

Dear List I am currently running a microarray analysis on the dataset schizo and would like to filter out all genes with a low variance. However, when running the code detailed below, I end up with 12,500 genes in my final set rather than the 10,000 I was looking for. Can anyone pinpoint where

Re: [R] identify extremes positions of values into matrix

I am glad to help ;-) Cheers, Francisco Milton Cezar Ribeiro wrote: > Francisco, > Thank you so much. It´s more than I need. > You help was of great value. > > Kind regards, > > Miltinho > > "Francisco J. Zagmutt" <[EMAIL PROTECTED]> escreveu: > Oi Milton, > > Is this wha

[R] microarray analysis

__ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] identify extremes positions of values into matrix

Francisco, Thank you so much. It´s more than I need. You help was of great value. Kind regards, Miltinho "Francisco J. Zagmutt" <[EMAIL PROTECTED]> escreveu: Oi Milton, Is this what you are after? x=read.table("miltondat.txt") >x V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 1 0 0 0 0 0 0 0

[R] How to convert an existing package?

Dears, how could I convert an existing package to be able to modify it? Let s say, if I want to change the way how the results of the package are display (to see the oroginal code of the code). Thanks [[alternative HTML version deleted]] __ R

Re: [R] identify extremes positions of values into matrix

Oi Milton, Is this what you are after? x=read.table("miltondat.txt") x V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 1 0 0 0 0 0 0 0 0 0 0 2 0 0 0 0 0 1 1 1 0 0 3 0 0 0 0 0 1 1 1 0 0 4 0 0 0 0 0 0 0 0 0 0 5 0 0 2 2 2 0 0 0 0 0 6 0 0 0 2 2 2 0 0

Re: [R] Translation of R code required

On 11/3/06, Barry Rowlingson <[EMAIL PROTECTED]> wrote: > How do you know that step works? If the dataframe 'mental' doesnt have > an 'Rx' column then mental$Rx will be 'NULL', and the factor() function > will make a factor out of that... Good point. I rarely use factors and missed that one. H

Re: [R] ssh and X11 forwarding

> Everything is enabled, but still doesn't work. > xeyes and xclock give errors (Error: Can't open display: localhost:10.0), > R produces Rplots.ps > > I also have a windows box, and go to chena using PuTTY. > I have cygwin with X server installed on it. > If I enable X11 forwarding in PuTTY, star

Re: [R] Translation of R code required

Sarah Goslee wrote: > Since this step works, > > mental$Rx <- factor(mental$Rx, levels=c("VS","IPS")) > > I think that some of the names are right, but perhaps that > one is spelled differently (Centre vs centre, maybe, since > R is case-sensitive?). How do you know that step works? If the da

Re: [R] Import problem with S-Plus 7.0 dataset

read.S is used to read S3 files (S-Plus 4 and older, not S-Plus 6 or newer) from the binary in the _Data directory. For example y <- read.S("c:/HOME/rmh/503.f99/_Data/y") data.restore is used to read S3 .sdd data.dump files. The S3 part means that S-Plus 7 users must force the oldStyle data.dump

Re: [R] Translation of R code required

Thank you for your help. Yes, it was a case-sensitive problem. I have run all the code now and have the requred output. Sarah Sarah Goslee wrote: > On 11/3/06, Sarah White <[EMAIL PROTECTED]> wrote: > >> >> I am using Windows XP, R version 2.4. >> >> The code works until I get to the fourth

Re: [R] Import problem with S-Plus 7.0 dataset

>From ?read.S 'read.S' can read the binary files produced in most recent versions of S-PLUS on either Windows (versions 3.x, 4.x, 2000) or Unix (version 3.x with 4 byte integers). It automatically detects whether the file was produced on a big- or little-endian machine

Re: [R] Translation of R code required

On 11/3/06, Sarah White <[EMAIL PROTECTED]> wrote: > > I am using Windows XP, R version 2.4. > > The code works until I get to the fourth line starting > > model.effect > > and I receive the following Error message. > > Error in eval(expr, envir, enclos) : object "Centre" not found > It sounds li

Re: [R] Translation of R code required

I am using Windows XP, R version 2.4. The code works until I get to the fourth line starting model.effect and I receive the following Error message. Error in eval(expr, envir, enclos) : object "Centre" not found When I say it works until the fourth line I mean I dont get an error message unt

Re: [R] Translation of R code required

> I am a complete novice when it comes to R. I have been sent some code > which I need to run on a dataset. I have downloaded and installed R > but I dont seem to be able to get the code to run. It would be very helpful if you would tell us your OS, version of R, and if you would copy and paste

[R] Translation of R code required

Dear All I am a complete novice when it comes to R. I have been sent some code which I need to run on a dataset. I have downloaded and installed R but I dont seem to be able to get the code to run. I was hoping someone on the list would be able to give me a 'translation' of what the code i

Re: [R] correaltion equal 1

On 11/3/06, Doran, Harold <[EMAIL PROTECTED]> wrote: > or, it could mean you need to recenter your time variable. > > > -Original Message- > From: [EMAIL PROTECTED] on behalf of Marc Bernard > Sent: Fri 11/3/2006 7:24 AM > To: r-help@stat.math.ethz.ch > Subject: [R] correaltion equal 1 > >

[R] Import problem with S-Plus 7.0 dataset

I am running R 2.3.1 on a Windows XP machine. I am trying to import some data into R that is stored as an S-Plus 7.0 .sdd file. When I run the following command, I get this error: > library(foreign) > d <- read.S(file='H:\\Research\\data.sdd') Error in read.S(file = "H:\\Research\\data.sdd") :

Re: [R] R-help in Newsreader?

On Fri, 2006-11-03 at 17:36 +0100, Matthias Voigt wrote: > Dear list, > > I am new to this all and therefore have following Newbie question: > > How can I receive and read R-help mailings in a newsreader like thunderbird? > > Thanks > Matthias If you go to the main R web site (http://www.r-proj

[R] R-help in Newsreader?

Dear list, I am new to this all and therefore have following Newbie question: How can I receive and read R-help mailings in a newsreader like thunderbird? Thanks Matthias _ __ R-hel

[R] CI -- fixed effects -- response -- mixed models

Dear Friends, I have been following the discussions of mixed-effects models with some confusion, and I realize that much of this is work in progress, and that much of the discussion is beyond my knowledge of statistics. My question is simple, though: Is there a set of commands that will pro

Re: [R] R CMD BATCH: unable to start device PNG

And on that note, here is a function that I use to get around it: -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jeffrey Horner Sent: Friday, November 03, 2006 10:01 AM To: [EMAIL PROTECTED] Cc: r-help@stat.math.ethz.ch Subject: Re: [R] R CMD BATCH: una

Re: [R] R CMD BATCH: unable to start device PNG

Whoops. I hit the wrong button! Here's the whole message: And on that note, here is a function that I use to get around it: startPNG <- function(filename, width = 680, height = 990, ...) { nativePNG <- capabilities("png") if(nativePNG) { png(file = filename, width = width, he

[R] Eilig Geschaftsvorschlag

Attn. Geschäftsvorschlag. Zuerst muß ich um Ihre Zuversicht in dieser verhandlung bitten. Das ist auf Grund seiner lage, als das sein total VERTRAULICH und Geheimnisvoll. Aber ich weiss, dass eine verhandlung dieses Ausmaßes irgendeinen Ängstlich und besorgt machen wird,aber ich versichre

Re: [R] generic arrays in R

On Fri, Nov 03, 2006 at 07:20:34AM +, Prof Brian Ripley wrote: > but > > a <- array(list(f),c(2,2)) > a[[2,1]] Thank you very much. Can you please tell me how to make apply generate arrays like this? Eg if a <- array(1:16,rep(2,4)) then a apply(a,c(1,2,4), polynomial) will just give me

Re: [R] maps. display hierararchical cluster analysis results on a map?

Besides symbols, size and color ramps come to mind Marsh Feldman -Original Message- From: David Farrar [mailto:[EMAIL PROTECTED] Sent: Thursday, November 02, 2006 9:37 PM To: R-help@stat.math.ethz.ch Subject: [R] maps. display hierararchical cluster analysis results on a map? Al

Re: [R] read file problem

Luis Ridao Cruz wrote: > R-help, > > I have the following file I want to import to R (some lines > removed) > > > Calibrated CTD data for station:00280001 > Calibrated:23/8 2001, Salinity Unsmoothed, Fluorescence Uncalibrated > Maximum observed depth:36 m

Re: [R] read file problem

On Fri, Nov 03, 2006 at 02:02:35PM +, Luis Ridao Cruz wrote: > DEPTH CTDPRS CTDTMP CTDSAL RAWFLU NUMB. QUAL > MDBAR IPTS-68 PSS-78 OBS. > *** *** *** *** > 1 1.0 2999 > 2 2.0 5.9793 35.1629.10717

Re: [R] read file problem

You should pay attention to the argument na.string. na.strings: a character vector of strings which are to be interpreted as 'NA' values. Blank fields are also considered to be missing values in logical, integer, numeric and complex fields. On 11/3/06, Luis Ridao Cru

Re: [R] read file problem

Hi If the file is tab delimited you could try something like this: a<-read.delim(file, skip = 9, header=F, na.strings=NA) Are you sure you want to skip 10 lines? (Is there a blank line somewhere?) J --- -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On

Re: [R] R CMD BATCH: unable to start device PNG

[EMAIL PROTECTED] wrote: > Dear r-helpers: > > Any ideas how to avoid problem described below? > I am having the same problem when I run R remotly (not from cgi script); > somehow png device wants to talk to X11 and X11 complains... > > Best regards, > Ryszard > > error msg from batch R process

Re: [R] help went away from top emacs window :-(

I located what might be the problem. There is a different misbehavior with Gnu emacs. I don't have Xemacs. I have an idea how to fix it but haven't worked out the details. ess-inf.el in the lines (defconst inferior-R-1-input-help (format "^ *help *(%s)" ess-help-arg-regexp)) (defconst inferior

Re: [R] help went away from top emacs window :-(

This change has created problems for me using XEmacs 21.4.19 with ESS 5.3.1 > ?help Error in print.help_files_with_topic("c:/progra~1/R/R-2.4.0/library/utils/chm/help") : CHM file could not be displayed > sessionInfo() R version 2.4.0 (2006-10-03) i386-pc-mingw32 locale: LC_COLLA

Re: [R] help went away from top emacs window :-(

I located what might be the problem. There is a different misbehavior with Gnu emacs. I don't have Xemacs. I have an idea how to fix it but haven't worked out the details. ess-inf.el in the lines (defconst inferior-R-1-input-help (format "^ *help *(%s)" ess-help-arg-regexp)) (defconst inferior

Re: [R] help went away from top emacs window :-(

On 11/3/2006 12:02 AM, Spencer Graves wrote: > This change has created problems for me using XEmacs 21.4.19 with > ESS 5.3.1 > > > ?help > > Error in > print.help_files_with_topic("c:/progra~1/R/R-2.4.0/library/utils/chm/help") > : > CHM file could not be displayed > > > sessionIn

[R] R CMD BATCH: unable to start device PNG

Dear r-helpers: Any ideas how to avoid problem described below? I am having the same problem when I run R remotly (not from cgi script); somehow png device wants to talk to X11 and X11 complains... Best regards, Ryszard error msg from batch R process (Linux/R-2.4.0) (executed on Apache server f

Re: [R] identify extremes positions of values into matrix

try this: # `mat' is your matrix apply(which(mat == 2, arr.ind = TRUE), 2, range) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/(0)16/3368

[R] identify extremes positions of values into matrix

Hi R-friends around world! I have a matrix that looks like bellow one. Each cell are filled with values 0, 1 or 2. In fact I have a lot of 512x512 matrices, with values varying from 0 to 999. I need to retrieve the left and right columns and the top and botton lines of the matrix where occours

[R] read file problem

R-help, I have the following file I want to import to R (some lines removed) Calibrated CTD data for station:00280001 Calibrated:23/8 2001, Salinity Unsmoothed, Fluorescence Uncalibrated Maximum observed depth:36 m QUAL has one digit for each of press

Re: [R] Using a deriv function in nlme

Lisa Avery sympatico.ca> writes: > > Hello, > > I have a deriv function that I am feeding to nlme. It works, and I can use > it in nls, but when I try to use it in nlme I get Error: subscript out of > bounds. I can fit the model using SSasympOrig, instead of the deriv It might by a similar pr

Re: [R] Rank transformation and the linear mixed model

Bruno L. Giordano music.mcgill.ca> writes: > > I am looking for references about mixed models built on rank transformed > data. To complement Frank H: Try lrm in package Design. Dieter __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/

[R] Error in X11(): Font Path?

Hy, I had the some problem. In my case an adjustment of the X11 font path, in the xorg.conf file did solve the problem. The fonts on my system were located in /usr/share/fonts and not in /usr/share/X11/fonts. Cheers Marc __ R-help@stat.math.ethz.c

Re: [R] Rank transformation and the linear mixed model

Bruno L. Giordano wrote: > - Original Message - From: "Frank E Harrell Jr" > <[EMAIL PROTECTED]> > To: "Bruno L. Giordano" <[EMAIL PROTECTED]> > Cc: > Sent: Friday, November 03, 2006 8:22 AM > Subject: Re: [R] Rank transformation and the linear mixed model > > >> Bruno L. Giordano wrote

Re: [R] difference in using with() and the "data" argument in glm call

It looks into the data argument first and if not found there it looks into environment(formula). Thus to pick it out of the 'with' try this: with(KK, glm(fo, data = environment()) Also you could try this: environment(fo) <- KK glm(fo) On 11/3/06, vito muggeo <[EMAIL PROTECTED]> wrote: > Dear

Re: [R] Rank transformation and the linear mixed model

- Original Message - From: "Frank E Harrell Jr" <[EMAIL PROTECTED]> To: "Bruno L. Giordano" <[EMAIL PROTECTED]> Cc: Sent: Friday, November 03, 2006 8:22 AM Subject: Re: [R] Rank transformation and the linear mixed model > Bruno L. Giordano wrote: >> Hello, >> I am looking for references

Re: [R] [R-pkgs] rcompletion: TAB completion for the R command line

Just a quick follow-up. Thanks to Simon Urbanek and hints from Alexandre Aguiar, CRAN now has a version (0.0-10) of rcompletion that should work on more systems. If it still doesn't work when you think it should, please let me know. -Deepayan On 10/30/06, Deepayan Sarkar <[EMAIL PROTECTED]> wrote

[R] difference in using with() and the "data" argument in glm call

Dear all, I am dealing with the following (apparently simple problem): For some reasons I am interested in passing variables from a dataframe to a specific environment, and in fitting a standard glm: dati<-data.frame(y=rnorm(10),x1=runif(10),x2=runif(10)) KK<-new.env() for(i in 1:ncol(dati)) assi

Re: [R] Rank transformation and the linear mixed model

Bruno L. Giordano wrote: > Hello, > I am looking for references about mixed models built on rank transformed > data. > Did anybody ever consider this topic? > > Thank you, > Bruno If you are speaking about the response variable, it's better to use a formal model such as the proportional odd

Re: [R] Compare linear regressios for significant differences of the slopes

Dieter Menne wrote: > Rainer M Krug wrote: >>> I have (8 measures * 96 groups) = 768 datasets for which I did linear >>> regressions using lm(). Sorry for my (probably incorrect) usage of "group" and "measure". Here a (probably better) description of my design: I ran 96 different simulations (sp

[R] Rank transformation and the linear mixed model

Hello, I am looking for references about mixed models built on rank transformed data. Did anybody ever consider this topic? Thank you, Bruno ~~ Bruno L. Giordano, Ph.D. CIRMMT Schulich School of Music, McGill University 555 Sherbrooke Street West Montréal, QC

Re: [R] Likelihood of multiple random processes

this way the problem is not well posed. You need to place an upper bound on the number of distributions in the mixture you want to estimate, for otherwise the likelihood is unbounded (take one normal distribution for each value in your data). - Original Message - > Do you know which obs

Re: [R] correaltion equal 1

or, it could mean you need to recenter your time variable. -Original Message- From: [EMAIL PROTECTED] on behalf of Marc Bernard Sent: Fri 11/3/2006 7:24 AM To: r-help@stat.math.ethz.ch Subject: [R] correaltion equal 1 Dear All, I wonder if this is a technical or an interpretation

[R] correaltion equal 1

Dear All, I wonder if this is a technical or an interpretation problem. I fitted an lme model including a random intercept and a random slope. The estimated correlation, from the lme, between the random intercept and the random slope is equal to 1. Does it mean that I should suppress on

[R] Forest plot

Need to modify ods ratios (like ods ratios from rmeta) so that they can line forest plot? How to do that? Branko! Branko Education inc. Brd. Icc. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do rea

[R] Using a deriv function in nlme

Hello, I have a deriv function that I am feeding to nlme. It works, and I can use it in nls, but when I try to use it in nlme I get Error: subscript out of bounds. I can fit the model using SSasympOrig, instead of the deriv function, but I am trying to reproduce an earlier analysis (done in Sp

Re: [R] Prediction intervals for predict.gls

Spencer Jones gmail.com> writes: > > I am using R 2.3.0, Is there a way to get prediction intervals using > predict.gls? > A r-project site search (alway recommended) would have given you http://finzi.psych.upenn.edu/R/Rhelp02a/archive/42932.html Dieter __

Re: [R] Installing a package - and making it work?

On Fri, 3 Nov 2006, Torsten Hothorn wrote: > > On Fri, 3 Nov 2006, Hvidberg, Martin wrote: > > > Dear list > > > > - Propoerly a beginner question, so bare with my frustration... > > > > I tried install the 'shapefiles' package into R 2.4.0 but it seems that the > > install had little effect...

Re: [R] Compare linear regressios for significant differenc es of the slopes

Rainer M Krug wrote: > > I have (8 measures * 96 groups) = 768 datasets for which I did linear > > regressions using lm(). > > > > Now I want to compare the slopes for each of the 8 measures in each of > > the 96 groups. As I understand , I can not use > >> anova(lm1, ..., lm8) > > as the lm1 ..

Re: [R] Installing a package - and making it work?

On Fri, 3 Nov 2006, Hvidberg, Martin wrote: > Dear list > > - Propoerly a beginner question, so bare with my frustration... Yes, and posting to R-help three minutes after posting to R-sig-geo, the correct list - raises the question about who's time is more important, doesn't it? What does he

Re: [R] Installing a package - and making it work?

On Fri, 3 Nov 2006, Hvidberg, Martin wrote: > Dear list > > - Propoerly a beginner question, so bare with my frustration... > > I tried install the 'shapefiles' package into R 2.4.0 but it seems that the > install had little effect... > > > > install.packages(c("shapefiles")) >

[R] Installing a package - and making it work?

Dear list - Propoerly a beginner question, so bare with my frustration... I tried install the 'shapefiles' package into R 2.4.0 but it seems that the install had little effect... > install.packages(c("shapefiles")) --- Please select a CRAN mirror for use in t

Re: [R] avoiding a loop: "cumsum-like"

Hi I have not seen any answer yet so I wil try (partly). I believe that the loop can be vectorised but I am a little bit lost in your fors and ifs. I found that first part of res is same as cumsum(tab$x.jour) until about 81st value. However I did not decipher how to compute the remaining part.