Re: [R] R code to generate Sweave input ? Sweave squared ?

Hi, Your idea reminds me of an example in the documentation of the brew package, featuring the generation of a template. You might want to check it out. baptiste On Wed, Nov 12, 2008 at 2:01 PM, Kem Phillips <[EMAIL PROTECTED]>wrote: > Dirk, > > > > I came upon your message below in searching f

[R] R not responding

Dear R help team, I am sorry to put urgent in the subject but to be honest I am in big troubles if I cannot make R to work in my computer. I am a PhD student in statistics and most of my homeworks and projects are in R. I have been using R in my mac book pro, mac os X 10.5.5 without problems unti

Re: [R] R not responding

Daniel Fernandez wrote: > Dear R help team, > I am sorry to put urgent in the subject but to be honest I am in big > troubles if I cannot make R to work in my computer. > > I am a PhD student in statistics and most of my homeworks and projects are > in R. > > I have been using R in my mac book pr

[R] Urgent: R not responding

Dear R help team, I am sorry to put urgent in the subject but to be honest I am in big troubles if I cannot make R to work in my computer. I am a PhD student in statistics and most of my homeworks and projects are in R. I have been using R in my mac book pro, mac os X 10.5.5 without problems unti

[R] Negative prediction by gam

Gam in mgcv package is predicting negative values which should not be the case despite all the predictors and response variables are positive. Tried to use log link function but it did not help. Please help S

Re: [R] Citing R in journal articles (or the failure to)

> > Could I suggest that citation() be modified to include > > the URL automatically? > > > I second this suggestion. I experienced similar case once as well. Thanks for pointing that function out to me. When I run it, the URL is included: citation() To cite R in publications use: R Devel

[R] REXcel problem

> Hi, > > I am trying to run a code in Excel using VBA from R, using > Rexcel.but everytime I am getting this error. Error - 2147220502 in module Recel.Rserver Error running expression Eval(parse(text="setwd(\"H:\\ > I am using R-2.8.0 and R.2.7 > Please help. > > > > << OLE Object

[R] (no subject)

Hi, Browse[1]> d4 EVDO_Rev Session_Setup FCA bin counts 50 NA 0 5 1 1 51 NA 0 5 2 1 52 NA 0 5 3 1 53 NA 0 5 4 1 54 NA 0 5 5 1 55 NA

[R] CROSSTABULATION

I want to form a 3x3 crosstabulation for the signs of two vectors (i.e. Negative, Zero, Positive). The problem is that I am simulating the data so for some iterations one of the categories is absent. Thus the resulting table shrinks to 3x2. I want it to be 3x3 with zero column corresponding to the

Re: [R] R-square in robust regression

yes thank you! it is perfect. I was using lmrob in package robustbase and it did not have that option in the summary. Laura 2008/11/13 Mark Difford <[EMAIL PROTECTED]> > > Hi Laura, > > >> I was searching for a way to compute robust R-square in R in order to > get > >> an > >> information simila

[R] how to plot a variable's histogram in the levels of a second variable

Goodevening to everyone, I am trying to create a clustered barplot for the following 2 variables who have more cases but as an example i am giving you 10 items: sex  socio-economic status 1     1 2     2 2     2 2     3 2     2 1     4 1     1 1     1 2    

Re: [R] Citing R in journal articles (or the failure to)

Achim Zeileis wrote: On Wed, 12 Nov 2008, Michael Friendly wrote: In two recent in press articles where I cited R packages, I was asked by the copy-editor to add a location or url to the references. Could I suggest that citation() be modified to include the URL automatically? That's a good

Re: [R] different results with plot.lm vs. plot.lm(which=c(2))

This was AFAICS a bug introduced in 2.7.1 by o plot(, which=5) uses more correct Cook's distance contours; (fix to fix to PR#9316). Thanks to Greg for the diagnosis. Will be fixed in R-patched later today. On Wed, 12 Nov 2008, Greg Snow wrote: Just a clarification on one of my

Re: [R] R-square in robust regression

> "LP" == Laura Poggio <[EMAIL PROTECTED]> > on Thu, 13 Nov 2008 10:43:14 + writes: LP> yes thank you! it is perfect. LP> I was using lmrob in package robustbase and it did not have that option in LP> the summary. Yes lmRob() from "robust" is from a company which

Re: [R] R-square in robust regression

I was searching for a way to compute robust R-square in R in order to get an information similar to the "Proportion of variation in response(s) explained by model(s)" computed by S-Plus. This post is dealing with that. Would be possible to have some hints on how to calculate this parameter within R

Re: [R] Urgent: R not responding

THere is a separate list for R and the mac: https://stat.ethz.ch/mailman/listinfo/r-sig-mac On Thu, Nov 13, 2008 at 3:38 AM, Daniel Fernandez <[EMAIL PROTECTED]> wrote: > Dear R help team, > I am sorry to put urgent in the subject but to be honest I am in big > troubles if I cannot make R to work

Re: [R] R-square in robust regression

Hi Laura, >> I was searching for a way to compute robust R-square in R in order to get >> an >> information similar to the "Proportion of variation in response(s) >> explained >> by model(s)" computed by S-Plus. There are several options. I have had good results using wle.lm() in package wle and

Re: [R] (no subject)

"Kurapati, Ravichandra \(Ravichandra\)" <[EMAIL PROTECTED]> napsal dne 13.11.2008 11:37:58: > But I observed that it was throwing an error. Really. If what you send is all your data then the problem is that EVDO_Rev is all NA. The problem is in your data, change your data and you will get resu

[R] Odp: (no subject)

Hi [EMAIL PROTECTED] napsal dne 13.11.2008 10:51:39: > Hi, > > > > Browse[1]> d4 > >EVDO_Rev Session_Setup FCA bin counts > > 50 NA 0 5 1 1 > > 51 NA 0 5 2 1 > > 52 NA 0 5 3 1 > > 53 NA

Re: [R] separate a variable in several variables

Thank you DAVID, GABOR, PHIL, WACEK, PETER and GREG I tried: Data2 <- data.frame(x=DATA$x, y = gsub('[^[:digit:]]', '', DATA$x), z= gsub('[^[:alpha:]]', '', DATA$x)) and everything is ok Wacek Kusnierczyk wrote: > > CE.KA wrote: >> >> x >> 1 12F >> 2 13 AD >> 3 356PO >> 4 1D >> 5

Re: [R] CROSSTABULATION

Hello... Which code are you using to perform the breakup into the three classes? Can you be more specific on that? Best regards, Carlos J. Gil Bellosta http://www.datanalytics.com On Thu, 2008-11-13 at 09:57 +, Sohail wrote: > I want to form a 3x3 crosstabulation for the signs of two vecto

Re: [R] how to plot a variable's histogram in the levels of a second variable

Argyro Antaraki yahoo.gr> writes: > I am trying to create a clustered barplot for the following 2 variables > who have more cases but as an example i am giving you 10 items: > > sex  socio-economic status > 1     1 > 2     2 > 2     2 > 2     3 .. > where sex: 1:man, 2:woman >

Re: [R] R-square in robust regression

I am aware of the limits of the parameter R^2 in this case. However often it is required for many different reasons. And it is helpful to have a function that does it. The most important is to know the drawback of the"number", I think. Laura 2008/11/13 Martin Maechler <[EMAIL PROTECTED]> >

[R] Surface plotting

Hi, I would like to plot a matrix which comprises 3 columns. So first column should be the x-axis. Second the y-axis and the third one should be the z-axis. I know that I can plot such data using the scatterplot3d, but I would like to have a surface plot like you can do with wireframe and persp.

[R] corfdrci and kendallfdrci from GeneNT

Dear all, I believe there is a problem with the corfdrci and the kendallfdrci function from the GeneNT package. The 2nd screening, that is taking what is not in the intersect of the CI and the user-defined MAS criterium does not work properly on my dataset. It seems to work properly on the da

[R] Calculate Specificity and Sensitivity for a given threshold value

Hi list, I'm new to R and I'm currently using ROCR package. Data in input look like this: DIAGNOSIS SCORE 1 0.387945 1 0.50405 1 0.435667 1 0.358057 1 0.583512 1 0.387945 1 0.531795 1 0.527148 0 0.526397 0 0.372935 1 0.8610

Re: [R] Surface plotting

On 11/13/2008 8:51 AM, mentor_ wrote: Hi, I would like to plot a matrix which comprises 3 columns. So first column should be the x-axis. Second the y-axis and the third one should be the z-axis. I know that I can plot such data using the scatterplot3d, but I would like to have a surface plot li

[R] Unflatten a table in R

Hi All, I'm pretty new to R, so would really appreciate it if someone could point me in the right direction on this problem. I am trying to "unflatten" a table in R, and can't seem to find a function or method to complete this task, (hopefully efficiently). My data table is full of historical st

Re: [R] Calculate Specificity and Sensitivity for a given threshold value

Kaliss wrote: Hi list, I'm new to R and I'm currently using ROCR package. Data in input look like this: DIAGNOSIS SCORE 1 0.387945 1 0.50405 1 0.435667 1 0.358057 1 0.583512 1 0.387945 1 0.531795 1 0.527148 0 0.526397 0 0.372935

Re: [R] Retrieving x argument name from a trellis object in R 2.8.0

"Deepayan Sarkar" <[EMAIL PROTECTED]> wrote on 13/11/2008 02:23:25: > I will try to track down the reason for the odd results that we > currently get if I have time, but it's not going to be a priority. > > -Deepayan > In any case, thanks. Vincenzo

[R] How to get txt file out of .rda file and vice versa ?

Ive got problems with preparing my own experiment data. I need to see what is the format, so there is a need of extracting rda file to tab delimited txt. After that I will probably have to get txt again to rda. How to do that ? I can use load() option but no clue how to save it (either txt or rda)

[R] readPDF() -- unsure how to install xpdf to make this work?

Dear R-Help, I need to convert a set of '.pdf' files into an equivalent set of '.txt' files. This is so that i can do some text mining on the content. In the latest R-News letter (http://cran.r-project.org/doc/Rnews/ Rnews_2008-2.pdf), the package 'tm' for text mining is mentioned. In that lovely

[R] Help with matplot

Hi, two days ago I have asked the list if there is a better way to plot a matrix using the following plot command: plot(matrix[1,], type="l") for (i in 2:dim(matrix)[1]) { lines(matrix[i,], type="l") } I have been told to use the matplot function, but unfortunately I matplot does not plo

[R] Boxplots with different variables and different ranges in R

Dear all, I am trying to make boxplots in one diagram but it is not working. Actually I have a large dataset, around 35 water quality variables with different units and of course in different ranges (some are orders of magnitude bigger than others). I want to produce one diagram with all the b

Re: [R] Boxplots with different variables and different ranges in R

We need to have some idea of what the data looks like and what code you have been trying to make any sensible suggestions. Can you supply a small example ? --- On Thu, 11/13/08, Lathouri, Maria <[EMAIL PROTECTED]> wrote: > From: Lathouri, Maria <[EMAIL PROTECTED]> > Subject: [R] Boxplots w

[R] uninitialized pdMat object

Dear all, I am using the nlme package and trying to modelate the variance- covariance matrix of random effects. However, I do not understand why do I get this error message: > M4CompSym <- update(M4,random=A~pdCompSymm(~CED)) Error in pdMatrix.pdCompSymm(x) : Cannot extract the matrix from

Re: [R] Unflatten a table in R

on 11/13/2008 09:49 AM Brigid Mooney wrote: > Hi All, > > I'm pretty new to R, so would really appreciate it if someone could point me > in the right direction on this problem. > > I am trying to "unflatten" a table in R, and can't seem to find a function > or method to complete this task, (hopef

Re: [R] How to get txt file out of .rda file and vice versa ?

?write.table ?read.table --- On Thu, 11/13/08, memento1984 <[EMAIL PROTECTED]> wrote: > From: memento1984 <[EMAIL PROTECTED]> > Subject: [R] How to get txt file out of .rda file and vice versa ? > To: r-help@r-project.org > Received: Thursday, November 13, 2008, 10:32 AM > Ive got problems wi

Re: [R] Help with matplot

on 11/13/2008 09:39 AM mentor_ wrote: > Hi, > > two days ago I have asked the list if there is a better way to plot a matrix > using > the following plot command: > > plot(matrix[1,], type="l") > for (i in 2:dim(matrix)[1]) { >lines(matrix[i,], type="l") > } > > I have been told to use t

[R] Date conversion

Hi List, If I have a date format as: d <- "2001/1/1" I can easily convert it to number by using as.Date(d). But if I have d<-"1/1/2001", it does not work. Does anyone know how I can convert it using pre-written function in R? Regards, Alireza [[alternative HTML version deleted]] ___

Re: [R] Calculate Specificity and Sensitivity for a given threshold value

Hi Frank, Thank you for your answer. In fact, I don't use this for clinical research practice. I am currently testing several scoring methods and I'd like to know which one is the most effective and which threshold value I should apply to discriminate positives and negatives. So, any idea for my

[R] replacing Na's in a vector with the preceding integer element

Hi, I have a vector where an integer is followed by a randomly varying number of NA elements. I want to replace the NA's with the preceding integer. Let me explain with a simple example : avec<- c(1,NA,NA,5,NA,3,4,NA,NA,NA,NA,9,NA,0,NA,NA) I want to get the following vector from a : bvec<-c(1,1,1

Re: [R] Date conversion

Dr. Alireza Zolfaghari wrote: > Hi List, > If I have a date format as: > d <- "2001/1/1" > I can easily convert it to number by using as.Date(d). Yes. If it means "January the 1st" and not "1st of January", that is... > But if I have d<-"1/1/2001", it does not work. Does anyone know how I can >

Re: [R] replacing Na's in a vector with the preceding integer element

See na.locf in the zoo package. It does exactly that. On Thu, Nov 13, 2008 at 12:07 PM, ravi <[EMAIL PROTECTED]> wrote: > Hi, > I have a vector where an integer is followed by a randomly varying number of > NA elements. I want to replace the NA's with the preceding integer. Let me > explain wit

Re: [R] R-square in robust regression

Hi Laura, The fact that you had copied Prof. Ripley's response suggested to me that this was so. Nevertheless, I think Martin Maechler was wise to emphasize the problem, just in case. Bye, Mark. Laura POggio wrote: > > I am aware of the limits of the parameter R^2 in this case. However often

Re: [R] Calculate Specificity and Sensitivity for a given threshold value

Hi Pierre-Jean, Sensitivity (Se) and specificity (Sp) are calculated for cutoffs stored in the "performance" x.values of your prediction for Se and Sp: For example, let's generate the performance for Se and Sp: sens <- performance(pred,"sens") spec <- performance(pred,"spec") Now, you can have a

[R] logical operators in 'subset'

Dear All, Am I correct in believing that logical operators such as && and || cannot be used as part of a 'subset' logical expression? Example: >foo = as.data.frame(matrix(c(1:9), nrow=3, ncol=3, dimnames=list(c("r1","r2","r3"),c("c1","c2","c3" >foo c1 c2 c3 r1 1 4 7 r2 2 5 8 r3

Re: [R] logical operators in 'subset'

2008/11/13 Reitsma, Rene - COB <[EMAIL PROTECTED]>: >>subset(foo, c1 > 1 && c2 > 5) > [1] c1 c2 c3 > <0 rows> (or 0-length row.names) Too many ampersands: > subset(foo, c1 > 1 & c2 > 5) c1 c2 c3 r3 3 6 9 See help("&") Barry __ R-help@r-projec

[R] problem installing and updating packages

I am having a problem installing and updating packages. It doesn't happen with every package I try to install or update, but most of them. I get this warning message: Warning: unable to move temporary installation 'C:\Program Files\R\R-2.8.0\library\file5f906952\rpart' to '[EMAIL PROTECTED] $2'

Re: [R] function that uses a variable name as the parameter

On 12/11/2008 7:59 AM, David Croll wrote: > > Hello dear R people! > > > Several times it occurred to me that a function that uses a variable name as a parameter would be helpful, but I did not find out how to write such a function. I have experience in several programming language, but I did no

Re: [R] Calculate Specificity and Sensitivity for a given threshold value

[EMAIL PROTECTED] wrote: Hi Frank, Thank you for your answer. In fact, I don't use this for clinical research practice. I am currently testing several scoring methods and I'd like to know which one is the most effective and which threshold value I should apply to discriminate positives and neg

[R] Plotting from a list

Hello, I am working with version 2.7.2 on a PC and have 2.8.0 available, but have not upgraded completely yet. I have 6 species for which I've run 6 unique environmental management scenarios against for comparison purposes. Each scenario is run for 36 years (1965 - 2000). Some of the species h

Re: [R] Calculate Specificity and Sensitivity for a given threshold value

N. Lapidus wrote: Hi Pierre-Jean, Sensitivity (Se) and specificity (Sp) are calculated for cutoffs stored in the "performance" x.values of your prediction for Se and Sp: For example, let's generate the performance for Se and Sp: sens <- performance(pred,"sens") spec <- performance(pred,"spec")

Re: [R] CROSSTABULATION

On Thu, 13 Nov 2008, Sohail wrote: I want to form a 3x3 crosstabulation for the signs of two vectors (i.e. Negative, Zero, Positive). The problem is that I am simulating the data so for some iterations one of the categories is absent. Thus the resulting table shrinks to 3x2. I want it to be 3x3

[R] Error in inherits(x, "ordered") : subscript out of bounds (MLSV, according to Wang et al., 2007)

Dear, Im trying to run the R program, provided by Wang (2007), for the analysis of nested frailty models. While running this program I get the following error: Error in inherits(x, "ordered") : subscript out of bounds Could someone tell me what the problem is in this case? Thanks in advance.

[R] Error in Quantile function

If anyone can assist with this problem you have my great thanks: I am trying to establish and plot confidence intervals on a bootstrapped function. I have a more complicated function that has no problems with determining the confidence intervals using the quantile command. This is outside the b

Re: [R] Plotting from a list

Yes you can automate it. The exact way would depend on how your data is structured and what you want the graphs to look like. It is hard to say without an example of the data and at least the command that your are using to create a single plot. On Thu, Nov 13, 2008 at 2:16 PM, <[EMAIL PROTECTED

Re: [R] Plotting from a list

Hi again, Yes you are right I should have included an example of the file It looks like this after reading in the data and using tapply I arrive at the following > Alt5rimc_mean A B C D E F 1965 8.423645 0.208 1.6296296 0.000 0.000

[R] 2^2 problem revisited

Dear R gurus: Here is the following from Montgomery's Design and Analysis of Experiments, 5th edition. > str(rout1.df) 'data.frame': 16 obs. of 3 variables: $ resp: num 18.2 18.9 12.9 14.4 27.2 24 22.4 22.5 15.9 14.5 ... $ A : Factor w/ 2 levels "-1","1": 1 1 1 1 2 2 2 2 1 1 ... $ B :

[R] growing a list sequentially -- memory management

hello, i have something like: out<-list() for(i in 1:n){ data<-gen(...) #fixed size data out[[i]]<- fun(data) } > object.size(out[[1]]) 6824 In principle 1 GB should allow n = 1024^3/6824 = 157347? i have about 2GB are not taken by other processes. however, I can see the memory shrinking q

[R] can you help how to solve the equation with intregration and complex ,

library(elliptic) f <- function (x) {myintegrate(function(t){Im(as.complex(1+(x*t)i))}, 1,2)-1} uniroot(f,lower = -2, upper = 2) > f <- function (x) {myintegrate(function(t){Im(as.complex(1+(x*t)i))}, > 1,2)-1} Error: unexpected symbol in "f <- function (x) {myintegrate(function(t){Im(as.comp

Re: [R] 2^2 problem revisited

Edna Bell wrote: Dear R gurus: Here is the following from Montgomery's Design and Analysis of Experiments, 5th edition. str(rout1.df) 'data.frame': 16 obs. of 3 variables: $ resp: num 18.2 18.9 12.9 14.4 27.2 24 22.4 22.5 15.9 14.5 ... $ A : Factor w/ 2 levels "-1","1": 1 1 1 1 2 2 2

Re: [R] Error in Quantile function

Please show us what you get when you execute: str(spead) typeof(spread) class(spread) It is doubtful that either of your theories is correct. It's more likely that spread is a more complex object than you realize. We need to know what spread looks like to the R interpreter in the environment

Re: [R] 2^2 problem revisited

Please read and understand the response that I already sent you on this issue. Rolf Turner On 14/11/2008, at 10:44 AM, Edna Bell wrote: Dear R gurus: Here is the following from Montgomery's Design and Analysis of Experiments, 5th edition. str(rout1.df) 'data.frame': 16 obs. of

Re: [R] Plotting from a list

on 11/13/2008 03:30 PM [EMAIL PROTECTED] wrote: > Hi again, > > Yes you are right I should have included an example of the file > > It looks like this after reading in the data and using tapply I arrive at > the following > >> Alt5rimc_mean > A B C D

Re: [R] Plotting from a list

Convert it to long form and plot like this (or look at ggplot2 for another approach): DF <- structure(list(A = c(8.423645, 30.881773, 29.783251, 17.827586, 27.507389, 17.46798, 8.167488, 32.502463, 44.610837), B = c(0.208, 27.344, 15.824, 31.776, 30.52, 27.232, 0, 19.736, 11.208), C = c(1.6296296,

Re: [R] read SAS file

You need to load the foreign package first. library(foreign) - Original Message - From: "b g" <[EMAIL PROTECTED]> To: Sent: Monday, November 10, 2008 12:28 PM Subject: [R] read SAS file Looked at a lot of documentation and listserv postings and still can't solve this problem. I ne

[R] lines.formula with subset

Dear list, When I try to use lines.formula with subset and another argument I get an error. e.g. x<-1:5 y<-c(1,3,NA,2,5) plot(y~x, type="n") # set up frame lines(y~x, subset=!is.na(y)) # works OK lines(y~x, type="o", col="blue") # works OK # but lines(y~x, subset=!is.na(y), col="red") # gi

[R] Splitting a lattice stripplot across several pages

Hello, I have a stripplot with 200+ labels i.e i have network connections. The y-axis are the server port numbers and am graphing the number of packets in a connection. Roughly stripplot(totpacks~portnum,data=network) Due to the large number of server ports in my dataset

[R] R run under Emacs fatal errors halt execution

I am running R.exe/Rterm.exe 2.7.1 under Emacs 22.1 under the DOS shell under Windows XP. (My Cygwin installation is broken!) I am not currently using ESS: > R.exe --vanilla R version 2.7.1 ... ... > 1+1 2 This works fine until I hit an R error -- pretty much any error -- when th

[R] radial.plot(plotrix) - plotting multiple polygons?

Is it possible to plot multiple polygons on radial.plot(plotrix)? The new=FALSE argument (as used to add additional points or lines with the plot function to an existing plot) doesn't appear to work. I would like to overlay the outlines of multiple polygons of different colors on the same radial pl

Re: [R] Splitting a lattice stripplot across several pages

On Thu, Nov 13, 2008 at 4:57 PM, Saptarshi Guha <[EMAIL PROTECTED]> wrote: > Hello, >I have a stripplot with 200+ labels i.e i have network connections. > The y-axis are the server port numbers and am graphing the number of > packets in a connection. Roughly >stripplot(totpacks~port

[R] aov help

Please pardon an extremely naive question. I see related earlier posts, but no responses which answer my particular question. In general, I'm very confused about how to do variance decomposition with random and mixed effects. Pointers to good tutorials or texts would be greatly appreciated. To

Re: [R] TimeZone Help - Finding TimeZone codes

On Monday 10 November 2008 01:23:30 pm Prof Brian Ripley wrote: > See ?Sys.timezone, which help.search("timezone") points you to. > > On Mon, 10 Nov 2008, stephen sefick wrote: > > I have looked at > > ?as.POSIXct > > That also links to Sys.timezone from the 'tz' argument. > > > ?POSIXct > > and ma

Re: [R] aov help

Rather than estimating the variance components via the method-of-moments estimators, have a look at the 'nlme' and 'lme4' packages, which provide likelihood-based tools for estimating random and mixed models (and in the case of nlme, gls models too). Advantages of the likelihood-based approaches in

Re: [R] sample size estimation for count (poisson?) data?

The notion that you can just add or subtract 0.03 from estimate is obviously incorrect. Presuming you meant to call you lower bound q05 and the upper bound q95, the numbers I get are in your 10,000 iteration loop are 4.97 and 5.18 (around a mean of 5.08). So roughly a .1 swing on each side

Re: [R] lines.formula with subset

> -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] > On Behalf Of John Field > Sent: Thursday, November 13, 2008 4:41 PM > To: r-help@r-project.org > Subject: [R] lines.formula with subset > > Dear list, > > When I try to use lines.formula with subset and anot

Re: [R] TimeZone Help - Finding TimeZone codes

On Thu, Nov 13, 2008 at 9:10 PM, J Dougherty <[EMAIL PROTECTED]> wrote: > On Monday 10 November 2008 01:23:30 pm Prof Brian Ripley wrote: >> See ?Sys.timezone, which help.search("timezone") points you to. >> >> On Mon, 10 Nov 2008, stephen sefick wrote: >> > I have looked at >> > ?as.POSIXct >> >>

[R] how to plot a variable's histogram in the levels of a second variable(1)

Goodmorning to everyone, I am trying to create a clustered barplot for the following 2 variables who have more cases but as an example i am giving you 10 items: sex  socio-economic status 1     1 2     2 2     2 2     3 2     2 1     4 1     1 1     1 2     3

[R] ellipse in pairs plot

Hi dear list, I'm scrambling with this: i have a p-dimensional data frame for which i want to have a pairs() plot, with ellipse added in each lower.panel. I want the ellipse's parameter (mu, sigma) to be extracted from the appropriate entries in a pre-computed matrix ("dab" in the code below)

[R] Symbols to use in text

Hello all, I'm trying to use text symbols in R, such are greek letters and mathematical operators. Do you know where I can find a list telling me how to use all those symbols? I'd like to have them displayed in my graphs. Best, Guillaume __ R-help

Re: [R] sample size estimation for count (poisson?) data?

Thanks. I did the search before I posted and found those threads. However, it does not seem to do what I want. All I want to do is estimate the sample size for a point estimate, not do a GLM. I just want the mean within a margin of error, and to a given CI. I've tried writing some code to do

Re: [R] Error in Quantile function

This is strange the first function that I am doing this for is saying in the typeof() is a double for some reason this is coming up as a list. > str(spread) 'data.frame': 18 obs. of 1 variable: $ bootdataframe.age: num 1.46e-09 1.67e-09 1.46e-09 1.45e-09 1.37e-09 ... > typeof(spread) [1] "

[R] negative prediction by gam (mgcv package)

Hi Gam in mgcv package is predicting negative values which should not be the case despite all the predictors and response variables are positive. Tried to use log link function but it did not help. Please help sunil -- View this message in context: http://www.nabble.com/negative-prediction

[R] question on clustered barplot

Thank you all for your help about the encoding of my yahoo.gr. I am sending you now the original data. Hope tha you can read them. I am trying to create a clustered barplot for the following 2 variables: > h$rocbiopsy1   [1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0

Re: [R] negative prediction by gam (mgcv package)

On Thu, 2008-11-13 at 21:32 -0800, sunil_POUDEL wrote: > Hi > > Gam in mgcv package is predicting negative values which should not be > the case despite all the predictors and response variables are positive. > Tried to use log link function but it did not help. Please help > > > sunil Pleas