Re: [R] R CMD check Error

Prof Brian Ripley wrote: > One thing I should add for the archives -- R 2.7.0 and later have > > oNew function showNonASCII() in package tools to aid detection > of non-ASCII characters in .R and .Rd files. Is there any chance of having the output of this show up in the test log wh

Re: [R] eval and parent.frame [was: Error in Design package: dataset not found for options(datadist)]

Charles C. Berry wrote: > On Tue, 22 Apr 2008, Gad Abraham wrote: > >> Charles C. Berry wrote: >>> On Sat, 19 Apr 2008, Gad Abraham wrote: >>> Charles C. Berry wrote: > On Fri, 18 Apr 2008, Gad Abraham wrote: > >> Frank E Harrell Jr wrote: >>> Gad Abraham wrote: H

Re: [R] matrix problem

Thanks very much Petr and Rold for your helpful replies. Cheers Bill __ R-help@r-project.org 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, mi

Re: [R] another matrix question

Thanks Julien & Rolf for your help. The double commas were just typos. Bill __ R-help@r-project.org 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 commen

Re: [R] Named arrays

I am using Bioconductor to pull chromosomal locations for affymetrix probe ids. library("annotate") z <- buildChromLocation("hgu95av2") vec <- chromLocs(z)[["Y"]] vec["266_s_at"] 266_s_at -19611913 Array vec is index by the ids and contain the respective locations. I would like to write to a fil

Re: [R] eval and parent.frame [was: Error in Design package: dataset not found for options(datadist)]

On Tue, 22 Apr 2008, Gad Abraham wrote: > Charles C. Berry wrote: >> On Sat, 19 Apr 2008, Gad Abraham wrote: >> >>> Charles C. Berry wrote: On Fri, 18 Apr 2008, Gad Abraham wrote: > Frank E Harrell Jr wrote: >> Gad Abraham wrote: >>> Hi, >> Design isn't strictly

Re: [R] Fwd: how to convert non numeric data into numeric?

"kakul modani" <[EMAIL PROTECTED]> wrote in news:[EMAIL PROTECTED]: > *@ Mr. Alspach and @ Mr. Winsemius : Thanks for replying. > > reqdIRR contains returns in lognormal form.I am trying to estimate > parameter for type-1 Gumbel Distribution and fit it into reqdIRR. > > If i use sapply its retu

Re: [R] eval and parent.frame [was: Error in Design package: dataset not found for options(datadist)]

Charles C. Berry wrote: > On Sat, 19 Apr 2008, Gad Abraham wrote: > >> Charles C. Berry wrote: >>> On Fri, 18 Apr 2008, Gad Abraham wrote: >>> >>> > Frank E Harrell Jr wrote: >>> > > Gad Abraham wrote: >>> > > > Hi, >>> > > > > > > Design isn't strictly an R base package, but maybe >>> someo

Re: [R] Multidimensional contingency tables

Now that is simple and elegant. Thanks! PS. Is there a course available for learning how to read R help information? :) At 10:52 PM 4/21/2008, Gabor Grothendieck wrote: > xtabs(count ~., prob1) > >On Mon, Apr 21, 2008 at 10:46 PM, Robert A. LaBudde <[EMAIL PROTECTED]> wrote: > > How does one ide

Re: [R] Multidimensional contingency tables

xtabs(count ~., prob1) On Mon, Apr 21, 2008 at 10:46 PM, Robert A. LaBudde <[EMAIL PROTECTED]> wrote: > How does one ideally handle and display multidimenstional contingency > tables in R v. 2.6.2? > > E.g.: > > > prob1<- data.frame(victim=c(rep('white',4),rep('black',4)), > + perp=c(rep('whit

[R] Fwd: how to convert non numeric data into numeric?

*@ Mr. Alspach and @ Mr. Winsemius : Thanks for replying. reqdIRR contains returns in lognormal form.I am trying to estimate parameter for type-1 Gumbel Distribution and fit it into reqdIRR. If i use sapply its returning me NA values. * I am having the following error in my function function(the

[R] Multidimensional contingency tables

How does one ideally handle and display multidimenstional contingency tables in R v. 2.6.2? E.g.: > prob1<- data.frame(victim=c(rep('white',4),rep('black',4)), + perp=c(rep('white',2),rep('black',2),rep('white',2),rep('black',2)), + death=rep(c('yes','no'),4), count=c(19,132,11,52,0,9,6,97)

Re: [R] Named arrays

Erhan <[EMAIL PROTECTED]> wrote in news:[EMAIL PROTECTED]: > Hi! I need to access all the strings used as indexes in an array. > Can someone tell me how to do that? > I'm not sure this is the right place to ask, but I couldn't find > another group. In R dataframe indices can be either row or col

Re: [R] Named arrays

Could you send an example please? Thanks, Erin On Mon, Apr 21, 2008 at 8:47 PM, Erhan <[EMAIL PROTECTED]> wrote: > Hi! I need to access all the strings used as indexes in an array. > Can someone tell me how to do that? > I'm not sure this is the right place to ask, but I couldn't find > another

Re: [R] how to convert non numeric data into numeric?

"kakul modani" <[EMAIL PROTECTED]> wrote in news:[EMAIL PROTECTED]: > I am having the following error in my function > > function(theta,reqdIRR) > { > theta1<-theta[1] > theta2<-theta[2] > n<-length(reqdIRR) > constant<- n*(theta1+theta2) > sum1<-lapply(reqdIRR*exp(theta1),FUN = sum) > sum2<-lap

Re: [R] how to convert non numeric data into numeric?

Kakul lapply() returns a list - as stated in the help file (in general, the help files and documentation provided with R are very good once one gets used to them). Lists are, potentially, complex things and it makes no sense to add them. Presumably in your case it does make sense to add them - i

[R] read variable global in R tcltk

I have any problem with my code. I build a small GUI in R with tcltk package. the proolem is I don't understand how to make a variable value can be read by R from a function. the variable value can only read if it's called via tcltk widgets and pass it to another function. i think the point is

[R] Named arrays

Hi! I need to access all the strings used as indexes in an array. Can someone tell me how to do that? I'm not sure this is the right place to ask, but I couldn't find another group. Thanks, Erhan __ R-help@r-project.org mailing list https://stat.ethz.ch

[R] how to convert non numeric data into numeric?

I am having the following error in my function function(theta,reqdIRR) { theta1<-theta[1] theta2<-theta[2] n<-length(reqdIRR) constant<- n*(theta1+theta2) sum1<-lapply(reqdIRR*exp(theta1),FUN = sum) sum2<-lapply(exp(theta2 - reqdIRR*exp(theta1)),FUN = sum) sum = sum1 + sum2 log.fcn = constant - as

Re: [R] How to turn a string into a variable name ?

Hello, following on - on this subject: i have tried the following: natural_nums <- 1:10 even_nums <- seq(2,10, by = 2) types <- c("log(natural_nums)~time", "even_nums") types <- lapply(types, as.name) ## list of variable names types [[1]] `log(natural_nums)~time` [[2]] even_nums as you

Re: [R] Documentation General Comments

Like all learning- you have to learn how to think like the people that create: language, poem, song, etc. In my own learning curve I have found that I am too often popping emails off to the list when in an hour or so I can figure out what was troubling me (using the provided examples). This may b

Re: [R] Documentation General Comments

Thank you Rob, I just downloaded it and it looks very useful. In the meantime I think I solved my immediate problem (and while pluggin' away also deepened my understanding - or so I will at least claim :-) Esmail > I don't know if you will find this helpful, but one of the better online > chap

Re: [R] Coding methods for factors

On 22/04/2008, at 11:59 AM, Paul Lynch wrote: > Is there a way in R to control how factors are coded? It looks like > reference cell coding is used by default. Is there a way to get > effects coding? > > (I could just set up the factor dummy variables as I wish and use lm, > but I'm wondering i

Re: [R] Documentation General Comments

The root of the problem is that R is a voluntary/cooperative project and those who develop and maintain R are (generously) contributing their time and probably have little-to-no time left over to devote to the improvement of the documentation. The generic answer to such a dilemma is that thos

Re: [R] Documentation General Comments

You can find all methods for an S3 class by doing this: methods(class = "recordedplot") For S4 see ?showMethods On Mon, Apr 21, 2008 at 5:56 PM, Beck, Kenneth (STP) <[EMAIL PROTECTED]> wrote: > I realize the R developers are probably overwhelmed and have little time > for this, but the documen

[R] Coding methods for factors

Is there a way in R to control how factors are coded? It looks like reference cell coding is used by default. Is there a way to get effects coding? (I could just set up the factor dummy variables as I wish and use lm, but I'm wondering if there is a better way.) Thanks, --Paul

[R] subset and plot

create a subset for "brook_dis". When I plot (week, R) I get a nice boxplot, but along the x axis, there are weeks a, b, c along with h and nh. Thank you ahead of time. keith rm(list=ls()) cond.exp1<-read.csv("condition/test.csv",header=TRUE) sub<-subset(cond.exp1, Species=="brook_d

Re: [R] Documentation General Comments

>> I realize the R developers are probably overwhelmed and have little time >> for this, but the documentation really needs some serious reorganizaton. >> A good through description of basic variable types would help a lot, >> e.g. the difference between lists, arrays, matrices and frames. > > Agre

Re: [R] rbugs on linux and wine

Hi Alexander! You are mixing WinBUGS and OpenBUGS. R package Rbugs works with OpenBUGS, but the later does not work with Rbugs under Linux! Gregor __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the post

Re: [R] Documentation General Comments

> I realize the R developers are probably overwhelmed and have little time > for this, but the documentation really needs some serious reorganizaton. > A good through description of basic variable types would help a lot, > e.g. the difference between lists, arrays, matrices and frames. Agree

Re: [R] getting off this mailing list

Go to https://stat.ethz.ch/mailman/listinfo/r-help There you will find how to do it by yourself. I hope this helps, Jorge On Mon, Apr 21, 2008 at 6:46 PM, Melissa Maxa <[EMAIL PROTECTED]> wrote: > Could I please be removed from this mailing list? > > Thanks, > > Melissa > > > ---

[R] getting off this mailing list

Could I please be removed from this mailing list? Thanks, Melissa - [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Documentation General Comments

On 21/04/2008 6:36 PM, Dr. Jeff Miller wrote: > I agree completely. Maybe it's time for an exhaustive manual (with weekly > downloadable updates, of course). > > It would also be nice if it were cross-referenced. For example, to get what > I wanted last weekend from a simple 2x2 contingency analy

Re: [R] Documentation General Comments

I agree completely. Maybe it's time for an exhaustive manual (with weekly downloadable updates, of course). It would also be nice if it were cross-referenced. For example, to get what I wanted last weekend from a simple 2x2 contingency analysis, I had to bounce between 4 different libraries. (All

Re: [R] Attempting to sync x-axis labels with grid lines in lattice bwplot

On 4/20/08, Mark Coletti <[EMAIL PROTECTED]> wrote: > I have created a lattice bwplot that suffers from axis label > overplotting. The box plot contains grid lines, and I thought that > only plotting the labels for corresponding vertical grid lines would > not only cure the overplotting but als

Re: [R] "spreadsheet" in plot

Try this: library(gplots) #after install.packages("gplots") m <- matrix(sample(9), 3) matplot(m, type='l') par(fig = c(0.5, 1, 0, 0.5), new = T) textplot(m, show.rownames = F, show.colnames = F, cex = 0.9) On 4/21/08, Stephan Ripke <[EMAIL PROTECTED]> wrote: > Hello, > how can I print a table (

Re: [R] change column names of several data frames

On Mon, 21 Apr 2008, Henrik Parn wrote: > Dear all, > > I have several data frames for which I want to change the column names. > > Example data: > data.1 <- data.frame(x1 = rnorm(5)) > data.2 <- data.frame(x1 = rnorm(5)) > . > . > > > What I want to achieve: > names(data.1) <- "y1" > names(data.1

[R] Use of recordPlot

I am trying to save plots that I create using calls to windows() folowed by multiple calls to plot to place several plots on a page. I want to store the plots as variables to be able to recall them later. The documentation states: "Recorded plot histories are of class "SavedPlots". They have a pri

[R] Documentation General Comments

I realize the R developers are probably overwhelmed and have little time for this, but the documentation really needs some serious reorganizaton. A good through description of basic variable types would help a lot, e.g. the difference between lists, arrays, matrices and frames. And, it appears ther

[R] "spreadsheet" in plot

Hello, how can I print a table (lets say my table xtable contains 3 rows and 3 Columns, the fileds consists of text) directly in a plot? when I use mtext(xtable), he tries to write everything on the same place. Thank you in advance, Stephan -- [[alternative HTML version deleted]] ___

[R] Summary of importing vector graphic formats written by R into Open Office

PROBLEM: OpenOffice (specifically Impress) imports emf files as the ones generated by the R GUI in Windows with poor quality. Is there another vector format written by an R graphic GUI that could be imported into OpenOffice (specifically Impress)? SOLUTIONS: Currently None. The suggested vector f

Re: [R] How to do survival analysis with time-related IVs?

Thanks. The first two links are helpful. For the log-reg, I was initially thinking that we need to re-arrange the dataset to make ID-Year observations (I think some books are suggesting this), and then run regression on "hazard rate". I did not know whether we should write code to calculate the h

Re: [R] another matrix question

On 21/04/2008, at 10:05 PM, William Simpson wrote: > Hi everyone, > > I would like to do the following. > > Given matrix m and matrix n, I would like to compute mn[i,,j]= m[i,,j] > + n[i,,j] if either of these elements is 0. (In other words, whichever > number is nonzero.) > Else I want mn[i,,j]=

Re: [R] matrix problem

On 21/04/2008, at 9:54 PM, William Simpson wrote: > Hi Everyone, > > I am running into a problem with matrices. I use R version 2.4.1 and > an older version. > > The problem is this: > m<-matrix(ncol=3,nrow=4) > m[,1:3]<-runif(n=4) > > That does what I expect; it fills up the rows of the matrix w

Re: [R] Change the core code

I often just download the source, find the appropriate function, create an file with an alternate version of it (i.e. plotMeans becomes plot.Means) and modify it to suit. I load all custom functions like that from my .rprofile. I guess you _could_ recompile the library, but, that might break oth

Re: [R] graphics history

Duncan Murdoch <[EMAIL PROTECTED]> wrote: > >> One thing I'd like to do, but didn't have time to implement before > >> 2.7.0, is to have history set to some finite size, e.g. a default might > >> be the last 3 or 10 plots. The problem with record=TRUE is that it > >> keeps a record of all the plo

Re: [R] finding an unknown distribution

andrea previtali wrote: > Hi, > I need to analyze the influences of several factors on a variable that is a > measure of fecundity, consisting of 73 observations ranging from 0 to 5. The > variable is continuous and highly positive skewed, none of the typical > transformations was able to normaliz

Re: [R] Regression inclusion of variable, effect on coefficients

Hello! I was thinking again about the possible interaction between x1 and x4. Theoretically it makes sense, that the influence of x4 on y is the stronger, the less informative is x1. It can be argued that the higher x1, the less informative it is x1. How could I incorporate this relationship in

Re: [R] Regression inclusion of variable, effect on coefficients

Hello :) I am happy to hear that I am not necessarily asking stupid questions. The thing is, that I have data on x1 and x4 for the whole sample. However, theoretically, it is clear that the informational content of x1 is not as high as of x4. x4 provides more accurate information to the subjects

Re: [R] Analysis of Epidemiological Data Using R

José Ignacio Bustos Melo wrote: > Hi everyone, > > I'm studying the manual name: Analysis of Epidemiological Data Using > R and Epicalc, maked by: Virasakdi Chongsuvivatwong and Edward McNeil. > > And I can't find the data base that they use in some examples, this are > the names: > > Chapter7.Rdat

[R] finding an unknown distribution

Hi, I need to analyze the influences of several factors on a variable that is a measure of fecundity, consisting of 73 observations ranging from 0 to 5. The variable is continuous and highly positive skewed, none of the typical transformations was able to normalize the data. Thus, I was thinking

[R] Change the core code

HI, pretty basic question: is that possible to change the code of the function within library? If so what should I do? I work on R linux (ubuntu), thanks a lot -- View this message in context: http://www.nabble.com/Change-the-core-code-tp16808285p16808285.html Sent from the R help mailing list

Re: [R] Adding number of non-NAs to boxplot

Have a look at the addtable2plot function in the plotrix package. It should do what you want. --- stephen sefick <[EMAIL PROTECTED]> wrote: > boxplot(x[,c(2,15,28,41,54,67,80,93,106)], > ylab="mg/s", names=c("RM215", > "RM202", "RM198", "RM190", "RM185", "RM179", > "RM148", "RM119", "RM61")) >

[R] jpeg legend space issues

Dear R community, I am printing a jpeg file (using plot) and my y-axis label becomes partly cut (at the left) by a very close margin of document. See example: http://www.igm.jhmi.edu/~gehret/progr_collect_data/beta.jpg Can you please help me fix this? I tried "din" and "fig" in the parameter an

Re: [R] Re ad From EXCEL

Hi there, Try this: # Function to read data in R from Excel FromExcel=function(yourfile,spreadsheet){ require(RODBC) channel=odbcConnectExcel(yourfile) sqlTables(channel) mydata=sqlFetch(channel, spreadsheet) attach(mydata) mydata } mydata=FromExcel("C:/mydata/2008/yourfile.xls","yourspreadsheet

[R] Analysis of Epidemiological Data Using R

Hi everyone, I'm studying the manual name: Analysis of Epidemiological Data Using R and Epicalc, maked by: Virasakdi Chongsuvivatwong and Edward McNeil. And I can't find the data base that they use in some examples, this are the names: Chapter7.Rdata,Chapter8.Rdata,Chapter9.Rdata Somebody can t

Re: [R] Avoiding a loop

Thank you so much to Jim and Mark for their advices. Now I solved the problem I had using a new approach. Best, Jorge On Mon, Apr 21, 2008 at 12:53 PM, jim holtman <[EMAIL PROTECTED]> wrote: > Will this do it for you: > > # Seed and data frames X, Y and Z > set.seed(123) > X=matrix(rnorm(300),

Re: [R] Data labels in barchart (lattice)

Hi again, Deepayan Sarkar wrote (21.4.2008): > Write your own panel function (which may or may not be a simple > exercise depending on your level of expertise in R). You could use > panel.barchart as a starting point. Basically, you need to insert > some calls to panel.text() (or something equival

Re: [R] Re ad From EXCEL

fp0 = "C:\\Documents and Settings\\myname\\My Documents\\Research\\" #upper file path fp1 = "PIname\\PIproject\\Working00\\" #middle file path fp2 = "DateData\\Dates.xls" #lower file path dbase = f

[R] Using the 'by' function withing a 'for' loop

Dear R experts, I am trying to optimize my script, because right now it requires a lot of memory. The goal is to generate four plots in one page. Every plot corresponds to the means and sem's calculated for a given variable at different days. In order to obtain the means and sem's I apply the

[R] Adding number of non-NAs to boxplot

boxplot(x[,c(2,15,28,41,54,67,80,93,106)], ylab="mg/s", names=c("RM215", "RM202", "RM198", "RM190", "RM185", "RM179", "RM148", "RM119", "RM61")) this is the code I am using to make a standard box plot. Is there a way to get the number of NA observations plotted onto the graph easily. I can alwa

Re: [R] UTF-8 or Unicode on Windows PC

On Mon, 21 Apr 2008, Hans-Joerg Bibiko wrote: On 21 Apr 2008, at 12:33, Prof Brian Ripley wrote: Is it possible to download a compiled snapshot of 2.7.0 for Windows XP? Yes, http://cran.r-project.org/bin/windows/base/rtest.html And it is due for release tomorrow. Many thanks! I can see th

Re: [R] Avoiding a loop

Will this do it for you: # Seed and data frames X, Y and Z set.seed(123) X=matrix(rnorm(300),ncol=5) Y=matrix(rpois(300,10),ncol=5) Z=matrix(rexp(300,1),ncol=5) index <- seq(1, by=3, length=nrow(X)) FINAL <- matrix(ncol=5, nrow=3*nrow(X)) FINAL[index,] <- X FINAL[index + 1,] <- Y FINAL[index + 2,

Re: [R] Regression inclusion of variable, effect on coefficients

This is not a dump question. This is a serious problem and it depends on what you know or assume about the relastionship between x1 and x4. If you assume linear interaction, you might want to introduce some interaction term to the model for example. Uwe Ligges Thiemo Fetzer wrote: > Hello dea

Re: [R] Data labels in barchart (lattice)

On 4/21/08, K. Elo <[EMAIL PROTECTED]> wrote: > Dear all, > > I use the barchart-function (lattice) for plotting stacked barcharts. > The data is a summary table (data frame) of likert-scale-evaluations > (strongly agree, agree...strongly disagree) to different issues > constructed as follows (

Re: [R] Creation of dialog box

There is a package in the R souces called "windlgs" in .../src/gnuwin32/windlgs which has examples, but you may want to go some way that is cross platform and avoid Windows specific programming at all. Best wishes, Uwe Ligges Ingrida B wrote: > Dear, List members, > > My student are creating

Re: [R] UTF-8 or Unicode on Windows PC

On 21 Apr 2008, at 12:33, Prof Brian Ripley wrote: >> Is it possible to download a compiled snapshot of 2.7.0 for Windows >> XP? > > Yes, http://cran.r-project.org/bin/windows/base/rtest.html > And it is due for release tomorrow. Many thanks! I can see the progress :) But please forgive my in

Re: [R] Design and analysis of mixture experiments

A summary, for those interested and posterity... Thanks to Christos Hatzis who is correct, the package 'AlgDesign' (which I'd overlooked) has gen.mixture which "Creates a candidate list of mixture variables". > gen.mixture(4,c("egg", "flour", "butter")) Thanks also to a private e-mailer who sugge

Re: [R] ANCOVA error again

On Mon, 2008-04-21 at 17:21 +0200, Birgit Lemcke wrote: > Hello Gavin, > > thanks for you answer. > If I use it without " with" I get back the same error. > The "with" thing was only to try out for functions that do not > contain a data-argument. I still try to learn and therefor I > sometimes

[R] Problem with graphic exhibition in Java using JRI

Dear all I am working on a Java class which calls R functions. To do that, I use JRI. Everything works fine up to now. The only problem I have is the following: when I call the command plot(….) from java, a blank "R Graphics" window appears and if I click on it the program does not answer any more.

[R] logit GLM without intercept

Dear Statisticians, I would like to analyse my data with a GLM with binomial error distribution and logit link function. The point is that I want a model fitted without intercept, i.e. the fitted curve should start at y=0.5 for x=0. I tried it with the following code: glm(value~0+ppm, bi

[R] Avoiding a loop

Dear R-users, I've been working with three different data sets (X, Y and Z) with the same dimension (i.e, n \times k). What I needed to do was to conform a 4th data set, i.e. FINAL, which first row was the X's first row, its second row was the Y's first row, and its third row was the Z's first row

[R] Bar Charts

Hello, Im trying to create a bar chart. I have a file with two different columns. I have created bar charts before where I am only reading from one column but now i wish to read from two columns. Here is the code i use when creating bar charts when reading from one column: satisfaction <-read.ta

Re: [R] ANCOVA error again

Hello Gavin, thanks for you answer. If I use it without " with" I get back the same error. The "with" thing was only to try out for functions that do not contain a data-argument. I still try to learn and therefor I sometimes just try. It is understood that I am on the way to simplify the mode

Re: [R] Re ad From EXCEL - Question

Erich: Past posts on this list have pointed out various pros and cons of different methods of data transer to R from Excel, in particular, loss of precision, formatting problems, etc. Do you have any comments about to what degree any of these alternatives may be susceptible or immune from these di

Re: [R] ANCOVA error again

On Mon, 2008-04-21 at 15:43 +0200, Birgit Lemcke wrote: > Hello R users! > > I got again an error message. Something here is causing compiled code to segfault ("crash"). I don't know what the problem is here exactly --- I'll let those much more acquainted with R look into that --- but you seem to

Re: [R] Labelling a secondary axis in R

If you just want the title, look at ?mtext. Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3265 http://www.StatisticalEngineering.com -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Monday, April 21, 2

[R] Array within an array

Hello, I need help to build an array within an array, i. e., I have this: tt[,,c(1,2)] , , metadados.class_7.R CS WRC LRA Inicial 1.000 1.000 1.000 Final 0.5974482 0.6095162 0.5866560 Indep 0.4335460 0.4799575 0.4169591 Inicial 0.9925572 0.9925572 0.99

Re: [R] Labelling a secondary axis in R

Apologies for the private mail, Nabble has not yet updated the thread so I can write another post in it. I think I have confused things. I don't mean the labels are incorrect. They are fine. What I am referring to is a title for the secondary axis, which is currently entitled as "c(-100,200)".

[R] pairs diagram of qq plots?

Hello everyone, for some exploratory analysis I would like to compare the distribution of an observable WERT pairwise between several samples identified by STICHPROBE (which differ in size). > str(stichproben_o1o4_20080327ff[c("STICHPROBE", "WERT")]) 'data.frame': 6087 obs. of 2 variables:

Re: [R] Labelling a secondary axis in R

On 4/21/2008 9:02 AM, Nakamura wrote: > Hello, > > How can I label a secondary axis in R? At the moment it's labelled as > c(-100,200). Obviously I would like it to be more sensible. > > Here is the code I am using > > newx = -100+37.5*((1:9)-1) > axis(4,at=newx,labels=(newx+100)/3750) I don't

[R] External regressors in GARCH

Hello to all R users, up to my knowledge, neither garch(tseries) nor garchFit(fGarch) support including external regressors in the regression, which, for example, arima(stats) can do by setting xreg. Is there a package or any other way that can do this? To be precise, I want to estimate a varianc

Re: [R] How to do survival analysis with time-related IVs?

"Tianxu" <[EMAIL PROTECTED]> wrote in news:[EMAIL PROTECTED]: > I am wondering how to do survival analysis with time-related IVs in > R. For example, See section 4 of Fox's contribution:

Re: [R] ANCOVA error again

Hello R users! I got again an error message. I used this code: with (FemMal85_Sex, { ModelFemMal85<- glm (Sex~outLatTep_like_other*outLatTep_like_conduplicate*outLatTep_keeled_w inged*spathellae_co

[R] Labelling a secondary axis in R

Hello, How can I label a secondary axis in R? At the moment it's labelled as c(-100,200). Obviously I would like it to be more sensible. Here is the code I am using newx = -100+37.5*((1:9)-1) axis(4,at=newx,labels=(newx+100)/3750) Thanks, Rob -- View this message in context: http://www.nabb

Re: [R] Equivalent of intervals() in lmer

Sorry, I meant to say: "For the moment I wonder if the solution is not to use CIs based on the two low SEs produced by the ~ time model, and to treat them as least-significant difference intervals." _ Professor Michael Kubovy University of Virginia Department of Ps

Re: [R] Equivalent of intervals() in lmer

Thanks Doug, You write: "If you want to examine the three means then you should fit the model as lmer(rcl ~ time - 1 + (1 | subj), fr)" I do just that (which is what Dieter just sent). But the CIs are much too big compared to the CIs for differences between means (which should be bigger tha

[R] estimate of overdispersion with glm.nb

Dear R users, I am trying to fully understand the difference between estimating overdispersion with glm.nb() from MASS compared to glm(..., family = quasipoisson). It seems that (i) the coefficient estimates are different and also (ii) the summary() method for glm.nb suggests that overdispersion is

Re: [R] Equivalent of intervals() in lmer

Douglas Bates stat.wisc.edu> writes: > If you want to examine the three means then you should fit the model as > lmer(rcl ~ time - 1 + (1 | subj), fr) > True, but for the notorious "error bars" in plots that reviewers always request the 0.35 is probable more relevant than the 1.87. Which I thin

Re: [R] Deleting rows with missing data

?complete.cases On Mon, Apr 21, 2008 at 8:43 AM, Charles Vetterli <[EMAIL PROTECTED]> wrote: > Hi folks: > > I have a data set v1, v2, ... v10. > Can anyone tell me how to create a new data set where the entire row is > deleted if, say, v5, is missing, is NA on that row? > > Thanks, > Charles > >

Re: [R] Deleting rows with missing data

?na.omit Weidong Gu, Department of Medicine University of Alabama, Birmingham 1900 University Blvd., Birmingham, Alabama 35294 Email: [EMAIL PROTECTED] PH: (205)-975-9053 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Charles Vetterli Sent: Monday, Ap

Re: [R] How to insert a vector or matrix into an existing matrix

On Mon, Apr 21, 2008 at 12:50:08PM +, David Winsemius wrote: [...] > > Am I correct in assuming that after the creation of m by way of a > temporary matrix that the temporary matrix would then be available for > garbage collection, whereas if both m and m2 were created, there would > be mor

Re: [R] ANCOVA

Hello John, I am really sorry about that. I wanted to include the code but I forgot and you are completely right, I forgot the family-argument. Thanks for the help. B. Am 21.04.2008 um 14:50 schrieb John Fox: > Dear Brigit, > > My guess is that you forgot to specify the argument family=binom

Re: [R] graphics history

On 4/21/2008 8:16 AM, Prof Brian Ripley wrote: > On Mon, 21 Apr 2008, Duncan Murdoch wrote: > >> On 21/04/2008 4:59 AM, Norbert NEUWIRTH wrote: >>> dear useRs and developeRs, >>> >>> I am afraid it is a very basic question, but I did not find anything alike >>> in the literature. >>> >>> The R st

[R] Regression inclusion of variable, effect on coefficients

Hello dear R users! I know this question is not strictly R-help, yet, maybe some of the guru's in statistics can help me out. I have a sample of data all from the same "population". Say my regression equation is now this: m1 <- lm(y ~ x1 + x2 + x3) I also regress on m2 <- lm(y ~ x

Re: [R] How to insert a vector or matrix into an existing matrix

Gabor Csardi <[EMAIL PROTECTED]> wrote in news:[EMAIL PROTECTED]: > On Sun, Apr 20, 2008 at 08:16:11PM +, David Winsemius wrote: >> Gabor Csardi <[EMAIL PROTECTED]> wrote in >> news:[EMAIL PROTECTED]: >> >> > Hmm, my understanding is different, >> > >> > m <- matrix(sample(10*10), ncol=10)

Re: [R] ANCOVA

Dear Brigit, My guess is that you forgot to specify the argument family=binomial in the call to glm(). Had you included the commands that you used as well as the error that was produced, it wouldn't be necessary to guess. I hope this helps, John On Mon, 21 Apr 2008 14:23:13 +0200 Birgit Lemck

[R] Deleting rows with missing data

Hi folks: I have a data set v1, v2, ... v10. Can anyone tell me how to create a new data set where the entire row is deleted if, say, v5, is missing, is NA on that row? Thanks, Charles __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

[R] Trend test for survival data

Hello, is there a R package that provides a log rank trend test for survival data in >=3 treatment groups? Or are there any comparable trend tests for survival data in R? Thanks a lot Markus -- Dipl. Inf. Markus Kreuz Universitaet Leipzig Institut fuer medizinische Informatik, Statistik und Epid

Re: [R] another matrix question

Hi, > Given matrix m and matrix n, I would like to compute mn[i,,j]= m[i,,j] > + n[i,,j] if either of these elements is 0. (In other words, whichever > number is nonzero.) > Else I want mn[i,,j]=(m[i,,j] + n[i,,j])/2 > I need a fast method. m <- matrix(c(0,1,2,3,4,0,5,6,0),nrow=3,ncol=3) n <- mat

  1   2   >