Re: [R] About RedR and RKWard

Hi Rookie, You're running RedR on Linux OR on Windows? If on Windows whether you download it on its official website? Or on R ran; install.packages(RedR) ? If on Linux whether you install RedR on repo? B.R. Stephen L - Original Message From: Rookie rookie8...@gmail.com To:

Re: [R] Learning ANOVA

Liaw, Andy andy_l...@merck.com on Mon, 16 Aug 2010 08:22:33 -0400 writes: From: Stephen Liu Hi JesperHybel, Thanks for your advice. If you're trying to follow the youtube video you have a typing mistake here: InsectSprays.aov

[R] Wilcoxon test and grouping factor with multiple levels

Dear R users, I have a dataset with two variables: $esan - a grouping factor with 8 levels and $reus. I'd like to do wilcox.test on this dataset as sugested Weiwei here: https://stat.ethz.ch/pipermail/r-help/2007-July/136627.html. I tried to adapt his recommendation but no succes. Can anyone help

Re: [R] A question of install R packages

hen chao Chang wrote: On Aug 16, 2010, at 11:05 PM, hen chao Chang wrote: My R_version is 2.11.1. I can’t install R packages. The error massages are following: package 'DBI' successfully unpacked and MD5 sums checked error in normalizePath(path) : path[1]=E:\R-2.11.1\library/DBI……

[R] Missing values

Hi folks, Could you please help me to understand; An Introduction to R http://cran.r-project.org/doc/manuals/R-intro.html#Logical-vectors The function is.na(x) gives a logical vector of the same size as x with value TRUE if and only if the corresponding element in x is NA. z -

Re: [R] Missing values

NA is a value that you can use for missing or inapplicable. It is also the value returned by R functions to indicate missing (e.g. if you use the match function to search a vector for a particular value, but the value isn't present, match will return NA). NaN means that the value wasn't missing

[R] AIC in MuMIn

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

[R] AIC in MuMIn

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

Re: [R] About RedR and RKWard

It was interesting to play with RedR, thanks for mentioning it. On Aug 16, 2010, at 10:15 PM, Stephen Liu wrote: Hi folds, RedR - visual programming for R http://www.red-r.org/ Has any folk used RedR before? What will be the difference compared with; RKWard - GUI for R

Re: [R] About RedR and RKWard

Not a problem. On Aug 16, 2010, at 11:49 PM, Stephen Liu wrote: Hi Rookie, Thanks for your info Stephen - Original Message From: Rookie rookie8...@gmail.com To: Stephen Liu sati...@yahoo.com Cc: r-help@r-project.org Sent: Tue, August 17, 2010 2:45:50 PM Subject: Re: [R] About RedR

Re: [R] About RedR and RKWard

I downloaded the Windows version. Although I have not played around with it long enough to give any meaningful feedback. At first glance, it seems that RedR would be a helpful tool to display some of the things that R is capable of doing. On Aug 16, 2010, at 11:30 PM, Stephen Liu wrote: Hi

Re: [R] A question of install R packages

Peter Dalgaard-2 wrote: hen chao Chang wrote: On Aug 16, 2010, at 11:05 PM, hen chao Chang wrote: My R_version is 2.11.1. I can’t install R packages. The error massages are following: package 'DBI' successfully unpacked and MD5 sums checked error in normalizePath(path) :

[R] ROCR predictions

Hi everybody, I am having a problem building a ROC curve with my data using the ROCR package. I have 10 lists of proteins such as attached (proteinlist.xls). each of the lists was calculated with a different p-value. The goal is to find the optimal p-value for the highest number of true

[R] Merge xts

Hi all , I have 12 xts objects of differing timeseries stamp. For example : str(s1_predict.xts) An ‘xts’ object from 1990-03-25 20:00:00 to 1990-12-15 09:00:00 containing: Data: num [1:725, 1] 11.23 10.18 9.3 9.74 10.18 ... - attr(*, dimnames)=List of 2 ..$ : NULL ..$ : NULL Indexed

Re: [R] Does anybody know how to control the appearance of the end of the line in lattice?

On Mon, Aug 9, 2010 at 1:19 AM, George Chen glc...@stanford.edu wrote: Hi All, I am plotting vertical lines using xyplot in lattice and type=h. It works well, but the problem is that the tops of the lines are convex and the bottoms are concave. Is there a way to flatten the tops and

Re: [R] creating pdf of wireframe

On Tue, Aug 10, 2010 at 1:55 AM, jculbertson culbert...@cogsci.jhu.edu wrote: Dear R list, I have written some code to produce several wireframe plots in a panel. They look good, but when I try to create a pdf, many (but not all) of the details I have specified are not reproduced. For

Re: [R] ROCR predictions

Dear Assa, I am having a problem building a ROC curve with my data using the ROCR package. I have 10 lists of proteins such as attached (proteinlist.xls). each of the your file didn't make it to the list. lists was calculated with a different p-value. The goal is to find the optimal

[R] ROCR data input

Hi there, I'm having some difficulty with the ROCR package. I've installed it fine, and the sample data works (ROCR.simple), however when I try to load my own data it complains that there is an error in prediction as the number of classes is not equal to 2. I read the data from a text file which

[R] Weird differing results when using the Wilcoxon-test

Hi, I became a little bit confused when working with the Wilcoxon test in R. As far as I understood, there are mainly two versions: 1) wilcox.test{stats}, which is the default and an approximation, especially, when ties are involved 2) wilcox_test{coin}, which does calculate the distribution

Re: [R] Lattice xyplots plots with multiple lines per cell

Here's a lattice solution using some faked data. library(lattice) testdat - data.frame(time = rep(rep(1:10, each = 4), 8), y = rnorm(320), gender = factor(rep(rep(c('F', 'M'), each = 80), 2)), grade = factor(rep(c('1-3', '4-6'),

Re: [R] Independent variables omitted in lm and glm

On Tue, 2010-08-17 at 10:29 +0800, elaine kuo wrote: Dear List, Some independent variable were missing in calculation using lm and glm (gaussian). (X= Y1+Y2+..+Y16, Independent number: 16 variable) What's this? This isn't valid R code. If it is `Elaine Kuo model notation` how on earth

[R] semiparametric fractional autoregressive model

folks, does anyone know if the SEMIFAR model has been implemented in R? i see that there's a S-FinMetrics function SEMIFAR() that does the job, but I have no access to that software. essentially, this semiparametric fractional autoregressive model introduces a deterministic trend to the

Re: [R] Wilcoxon test and grouping factor with multiple levels

On Aug 17, 2010, at 2:48 AM, Iurie Malai wrote: Dear R users, I have a dataset with two variables: $esan - a grouping factor with 8 levels and $reus. I'd like to do wilcox.test on this dataset as sugested Weiwei here: https://stat.ethz.ch/pipermail/r-help/2007-July/136627.html. I tried to

Re: [R] ROCR data input

Sorry, I'm new to R, and relatively new to statistics too so I'm still a bit unclear. The values in the post were only a sample of around 8400 rows. The label has 1 or 0 (I thought this was the two classes needed). Each label row has an equivalent probability. This is the data that I output from

[R] ROCR data format

Hi there, I'm having some difficulty with the ROCR package. I've installed it fine, and the sample data works (ROCR.simple), however when I try to load my own data it complains that the number of predictions should be the same as the number of labels for each run. I imported as a txt file with

[R] replacing values in a vector

Dear helpRs Does anyone have an elegant way of doing the following: For a given numeric vector, e.g. vec - c(3,2,6,4,7) Create a series of vectors where all but 1 of the values are replaced by 0's, e.g. vec.a - c(3,0,0,0,0) vec.b - c(0,2,0,0,0) vec.c - c(0,0,6,0,0) vec.d - c(0,0,0,4,0)

[R] Automated plot and linear regression line/data

Hi, I have an excel sheet (already imported to R) with multiple columns and I am looking for a way in R that will allow me to generate a plot for every possible pair and its linear regression line/data. Any tip/idea/script how do to so. Thanks, As hz -- View this message in context:

[R] confidence intervals for Harrell's c-index in survival setting

Dear All, Is it possible to get confidence intervals for Harrell's concordance index or, equivalently, Somer's D using the rms package or in some other way ? I have survival data it would be the c-index in the Cox model setting Many thanks Dr Bernard North Statistical Consultant Statistical

[R] Array Comparision

Hello. Can I just know how to compare two different arrays, but take only the elements from first array (both arrays are multidimensional). Its just opposite of intersect option in R. Thanks in advance -- View this message in context:

Re: [R] xyplot - I can't find colors in my picture

Thank you very much, Chris! Your advice was very helpful for me ;) -- View this message in context: http://r.789695.n4.nabble.com/xyplot-I-can-t-find-colors-in-my-picture-tp2326848p2328236.html Sent from the R help mailing list archive at Nabble.com.

[R] Trimmed Sample

I know about the trimmed mean function. Is there a function that trims values of an array, vector ect cheers -- View this message in context: http://r.789695.n4.nabble.com/Trimmed-Sample-tp2328139p2328139.html Sent from the R help mailing list archive at Nabble.com.

[R] predict.lm, matrix in formula and newdata

Dear all, I am stumped at what should be a painfully easy task: predicting from an lm object. A toy example would be this: XX - matrix(runif(8),ncol=2) yy - runif(4) model - lm(yy~XX) XX.pred - data.frame(matrix(runif(6),ncol=2)) colnames(XX.pred) - c(XX1,XX2) predict(model,newdata=XX.pred) I

Re: [R] A question of install R packages

Hi: I run into this problem occasionally (Windows). My 'solution' is to just reinstall the offending package (it's never failed for me a second time) and then reinstall/update the packages that followed it. Annoying, but no biggie. HTH, Dennis On Mon, Aug 16, 2010 at 8:05 PM, hen chao Chang

Re: [R] replacing values in a vector

Hi Karen, Try this: vec - c(3,2,6,4,7) res - diag(vec) k - length(vect) for(i in 1:k) assign(paste('vec', letters[i], sep = '.'), res[i,]) vec.a # [1] 3 0 0 0 0 vec.b # [1] 0 2 0 0 0 HTH, Jorge On Tue, Aug 17, 2010 at 6:57 AM, Karen Kotschy wrote: Dear helpRs Does anyone have an elegant

Re: [R] Wilcoxon test and grouping factor with multiple levels

Hi: As the help page for wilcox.test() states (?wilcox.test), the test is meant for one or two groups. Since you have eight groups in your data, the help page for wilcox.test() suggests using kruskal.test() instead, where the latter function applies the Kruskal-Wallis test if your intention is to

Re: [R] replacing values in a vector

Try this: mapply(function(x, y)assign(x, y, envir = globalenv()), sprintf('vec.%s', letters[1:length(vec)]), split(diag(vec), 1:length(vec))) On Tue, Aug 17, 2010 at 7:57 AM, Karen Kotschy ka...@sevenc.co.za wrote: Dear helpRs Does anyone have an elegant way of doing the following: For a

Re: [R] replacing values in a vector

Hi, You can use diag() and matrix multiplication to create a matrix with the requested rows, and assign() to store the rows as separate vectors: vec - c(3,2,6,4,7) mat - diag(vec) for (i in seq_along(vec)) assign(paste(vec, letters[i], sep=.), mat[i,]) vec.a

Re: [R] predict.lm, matrix in formula and newdata

try it better this way: XX - matrix(runif(8), ncol = 2) DF - as.data.frame(XX) DF$yy - runif(4) model - lm(yy ~ ., DF) XX.pred - as.data.frame(matrix(runif(6), ncol = 2)) predict(model, XX.pred) I hope it helps. Best, Dimitris On 8/17/2010 2:24 PM, Stephan Kolassa wrote: Dear all, I am

Re: [R] predict.lm, matrix in formula and newdata

Hi Stephan, You'll get the expected result if you pass a data.frame to the lm() function when fitting the model. yy - runif(4) XX - matrix(runif(8),ncol=2) df - data.frame(y=yy, XX) df yX1X2 1 0.52889284 0.8055476 0.6670006 2 0.09989951 0.2498907 0.8867955 3

Re: [R] replacing values in a vector

try this: vec - c(3,2,6,4,7) n - length(vec) for(i in seq_along(vec)){ r - numeric(n) r[i] - vec[i] assign(paste(vec., letters[i], sep = ), r) } I hope it helps. Best, Dimitris On 8/17/2010 12:57 PM, Karen Kotschy wrote: Dear helpRs Does anyone have an elegant way of doing

Re: [R] Automated plot and linear regression line/data

On 08/17/2010 11:47 AM, ashz wrote: Hi, I have an excel sheet (already imported to R) with multiple columns and I am looking for a way in R that will allow me to generate a plot for every possible pair and its linear regression line/data. Any tip/idea/script how do to so. Thanks, As hz

[R] dims error

Greetings, I am a very novice user with R, and in the course of running a linking procedure : P.old- function (a, c, b, xi){ #a contains a parameters #c contains c parameters #b contains b parameters #xi is a one column vector

[R] remove rows ??

Hello I have a file with this format 2005-01-03 09:05 0.00 2005-01-03 09:10 0.01 2005-01-03 09:15 0.02 2005-01-03 09:20 0.03 2005-01-03 09:25 0.04 2005-01-03 09:30 0.05 2005-01-03 09:35 0.06 2005-01-03 09:40 0.07 2005-01-03 09:45 0.08 2005-01-03

Re: [R] dims error

Hi, Maybe knowing what you provided to the function call as a, b, c, and xi would help! Since you're at it, think about using the dput() function to provide copy/pastable objects ;) HTH, Ivan Le 8/17/2010 15:00, Yaacov Petscher a écrit : Greetings, I am a very novice user with R, and in

Re: [R] remove rows ??

Try this: #DF - read.table( ) DF[!as.POSIXlt(DF$V1)[['min']] %% 15,] On Tue, Aug 17, 2010 at 10:02 AM, skan juanp...@gmail.com wrote: Hello I have a file with this format 2005-01-03 09:05 0.00 2005-01-03 09:10 0.01 2005-01-03 09:15 0.02 2005-01-03 09:20 0.03

Re: [R] ROCR data input

Anneley, Sorry, I'm new to R, and relatively new to statistics too so I'm still a bit unclear. That's OK - everyone started some time and was new. However, it is really important to post a reproducible example here. If you are so new that you don't know how to do that exactly, you should

Re: [R] confidence intervals for Harrell's c-index in survival setting

On Aug 17, 2010, at 6:10 AM, North, Bernard V wrote: Dear All, Is it possible to get confidence intervals for Harrell's concordance index or, equivalently, Somer's D using the rms package or in some other way ? I have survival data it would be the c-index in the Cox model setting I

Re: [R] Missing values

- Original Message From: Michael Bedward michael.bedw...@gmail.com To: Stephen Liu sati...@yahoo.com Cc: r-help@r-project.org Sent: Tue, August 17, 2010 3:57:05 PM Subject: Re: [R] Missing values NA is a value that you can use for missing or inapplicable. It is also the value returned

Re: [R] AIC in MuMIn

On Tue, 2010-08-17 at 16:05 +0800, elaine kuo wrote: Hello, Why did you decide to post the exact same message from two different email addresses?? 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

Re: [R] dims error

Yaacov, I cannot really help you, I don't even understand what you're trying to do! What's the link between your function and all the information you've just sent (by the way, you forgot to send it to the list)? By the way, attachments are not really advised, you'd better use dput() I've

Re: [R] Specify decimal places for parameters in BUGS output

On 16.08.2010 23:01, Anamika Chaudhuri wrote: Hi All: I had a basic question to ask. I am running R2WinBUGS so that I could automate the running of my model using 1000 simulated datasets. Below is the code I am using. The only problem I am having is the bugs output that comes out shows my

Re: [R] Odd crash with tcl/tk

Hi @ Andrew: Did you get this problem solved? I am having similar problems and have tried to work around it using options(gsubfn.engine = R) as suggested by Gabor. However, this don't solve on my machine. R just freezes up trying to execute read.csv.sql() On OSX 10.6.4, R 2.11.1 (2010-05-31)

[R] using color palette in persp

Hello! I am using persp to get a plot of a 32x48 matrix. data - as.matrix(read.table(file=mySample.dat)) persp(data) And I'm having trouble with trying to give a z-dependent palette to it like rainbow (or any built in one would do). How to combine that with the col argument? nice day Szilvi

Re: [R] Odd crash with tcl/tk

Sorry to chime in late (for some reason, I was not at the Mac when this came in originally). I can't reproduce the situation. One guess is that you haven't installed the tcltk bits, as indicated on the http://cran.at.r-project.org/bin/macosx/ page. -pd On Aug 17, 2010, at 4:12 PM, Lars

[R] Missing values in function calls

Hello, I would like to ask you how to determine if function value is missing for particular function call. I mean if I define function f __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] \ell symbol (log-likelihood)

Dear R experts---is it possible to plot the \ell symbol in R under the pdf device?  the following did not work: pdf(file=ell.pdf); plot( c(0,1), c(0,1) ); text( 0.5, 0.5, \u2113 ) dev.off() my guess is that this cannot be done, but I thought I would ask. Ivo Welch (ivo.we...@brown.edu,

Re: [R] Automated plot and linear regression line/data

Hi, I presume the easy way is using plot and lm. Thx, As hz -- View this message in context: http://r.789695.n4.nabble.com/Automated-plot-and-linear-regression-line-data-tp2328027p2328372.html Sent from the R help mailing list archive at Nabble.com.

[R] [R-pkgs] Major update of Deducer, including GSOC funded work

A new version of Deducer (0.4) has been sent to CRAN, and should be propagating to all the mirrors in due course. Also, a new plug-in package (DeducerExtras) has also been released to CRAN, containing additional dialogs and functionality. I would like to take this opportunity to thank

Re: [R] remove rows ??

I get this message character string is not in a standard unambiguous format and if I try to use a zoo object Error en `$.zoo`(final, V2) : not possible for univariate zoo series cheers -- View this message in context: http://r.789695.n4.nabble.com/remove-rows-tp2328312p2328487.html Sent from

Re: [R] remove rows ??

Works for me: Lines - 2005-01-03 09:05;0.00 2005-01-03 09:10;0.01 2005-01-03 09:15;0.02 2005-01-03 09:20;0.03 2005-01-03 09:25;0.04 2005-01-03 09:30;0.05 2005-01-03 09:35;0.06 2005-01-03 09:40;0.07 2005-01-03 09:45;0.08 2005-01-03 09:50;0.09 2005-01-03 09:55;0.10 2005-01-03 10:00;0.00 2005-01-03

Re: [R] remove rows ??

Now it works I was reading from a file DF - read.table(prueba.txt) DF[!as.POSIXlt(DF$V1)[['min']] %% 10,] didn't do anything, then I tried zz[!as.POSIXlt(zz$V2)[['min']] %% 10,] that caused an error Error en as.POSIXlt.character(as.character(x)) : character string is not in a standard

Re: [R] Delete rpart/mvpart cross-validation output

On Fri, 2010-08-13 at 14:52 -0300, Marie-Hélène Ouellette wrote: Dear all, I've been away for several weeks but I don't see an answer on the list so... I was wondering if there is a simple way to avoid printing the multiple cross-validation automatic output to the console of recursive

Re: [R] question about bayesian model selection for quantile regression

Hi All: the package MuMIn can be used to select the model based on AIC or AICc. The code is as follows: data(Cement) lm1 - lm(y ~ ., data = Cement) dd - dredge(lm1,rank=AIC) print(dd) If I want to select the model by BIC, what code do I need to use? And when to select the best model based on

[R] Plotproblem, histogram, y-axis label?

How can I control label on y-axis? Specifically how high relative to the y-axis it is plotted? At the moment I am using mtext appearing in the function quotet below. The function is applied to a dataframe and plots histograms for all numerical variables of the frame generating several

Re: [R] Odd crash with tcl/tk

Hi Peter Spot on! You where absolutely right. After installing the tcltk bits it worked just fine. Thank you very much! Lars On Aug 17, 4:36 pm, peter dalgaard pda...@gmail.com wrote: Sorry to chime in late (for some reason, I was not at the Mac when this came in originally). I can't

[R] how to paste a vector to expression?

hello, i need to annotate a plot with an expression and a variable value - like: plot(1:4) dat-1:2 text(2:3,2:3,expression(paste(bar(x)==dat))) ... but with the dat values plotted, 1 at x=2,y=2, and 2 at x=3,y=3. maybe someone can help? yours, kay - Kay

Re: [R] Weird differing results when using the Wilcoxon-test

After fixing the parentheses in your code so it does run, it seems that the difference is that wilcox.test defaults to using a continuity correction and your manual calculation does not. Use wilcox.test(big1, big2, correct=FALSE). -thomas On Tue, 17 Aug 2010, Cedric Laczny wrote: Hi,

Re: [R] Odd crash with tcl/tk

That should not be necessary if you used options(gsubfn.engine = R) since it should not be using tcltk in the first place. Perhaps you are using an old version of gsubfn -- the latest on CRAN is 0.5-3. On Tue, Aug 17, 2010 at 11:49 AM, Lars Dalby lars.da...@gmail.com wrote: Hi Peter

Re: [R] using color palette in persp

On Aug 17, 2010, at 10:35 AM, szisziszilvi wrote: I am using persp to get a plot of a 32x48 matrix. data - as.matrix(read.table(file=mySample.dat)) persp(data) And I'm having trouble with trying to give a z-dependent palette to it like rainbow (or any built in one would do). How to combine

Re: [R] Odd crash with tcl/tk

On Tue, Aug 17, 2010 at 11:49 AM, Lars Dalby lars.da...@gmail.com wrote: Hi Peter Spot on! You where absolutely right. After installing the tcltk bits it worked just fine. Thank you very much! Lars On Aug 17, 4:36 pm, peter dalgaard pda...@gmail.com wrote: Sorry to chime in late (for

Re: [R] Missing values in function calls

On Aug 17, 2010, at 8:22 AM, vo...@cbox.cz wrote: Hello, I would like to ask you how to determine if function value is missing for particular function call. I mean if I define function ?is.na f David Winsemius, MD West Hartford, CT __

[R] how to selection model by BIC

Hi All: the package MuMIn can be used to select the model based on AIC or AICc. The code is as follows: data(Cement) lm1 - lm(y ~ ., data = Cement) dd - dredge(lm1,rank=AIC) print(dd) If I want to select the model by BIC, what code do I need to use? And when to select the best model based on

Re: [R] ROCR predictions

Dear Claudia, thank you for your fast answer. I add again the table of the data as an example. Protein ID Pfam Domain p-value Expected Is Expected True Postive False Negative False Positive True Negative NP_11.2 APH 1.15E-05 APH TRUE 1 0 0 0 NP_11.2 MutS_V 0.0173 APH FALSE 0 0 1 0

Re: [R] predict.lm, matrix in formula and newdata

-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Stephan Kolassa Sent: Tuesday, August 17, 2010 5:25 AM To: R-help@r-project.org Subject: [R] predict.lm, matrix in formula and newdata Dear all, I am stumped at what should

[R] frailtypack

Hello everyone, I have a question about using frailtypack to run recurrent events data. I have a file with 1,580 events. The file structure fits the andersen gill model. When I run the file with all records, the R stops and shuts down. Once I run only part of the file, about 400 events, R runs

Re: [R] Wilcoxon test and grouping factor with multiple levels

Thank you, David Winsemius and Dennis Murphy! I want to compare pairs one by one. 2010/8/17 Iurie Malai iurie.ma...@gmail.com: Dear R users, I have a dataset with two variables: $esan - a grouping factor with 8 levels and $reus. I'd like to do wilcox.test on this dataset as sugested Weiwei

Re: [R] grid.table and expression in table body?

Hi - I can't get this figured out ... Thanks for any hint. Joh load(/tmp/AbsoluteTable.Rdata) absolutetable library(gridExtra) grid.table(absoluteTable)#Works grid.table(absoluteTable,parse=TRUE) Error in parse(text = d[ii]) : unexpected symbol in Survey Scans sessionInfo() R version

[R] navel-gazing

month - 2010-August list - r-help ##list - r-sig-ecology ##list - r-sig-mixed-models ## month - 2010q3 n - 50 baseurl - https://stat.ethz.ch/pipermail/; library(RCurl) z - getURL(paste(baseurl,list,/,month,/author.html,sep=)) zz - strsplit(z,LI)[[1]] namefun - function(x) {

Re: [R] \ell symbol (log-likelihood)

On Aug 17, 2010, at 8:41 AM, ivo welch wrote: Dear R experts---is it possible to plot the \ell symbol in R under the pdf device? the following did not work: pdf(file=ell.pdf); plot( c(0,1), c(0,1) ); text( 0.5, 0.5, \u2113 ) dev.off() my guess is that this cannot be done, but I thought

Re: [R] how to selection model by BIC

Xin__ Li msf09xl at mail.wbs.ac.uk writes: Hi All: the package MuMIn can be used to select the model based on AIC or AICc. The code is as follows: data(Cement) lm1 - lm(y ~ ., data = Cement) dd - dredge(lm1,rank=AIC) print(dd) If I want to select the model by BIC, what code do I

[R] R Send an Email

Just out of curiosity, has anyone ever written a function that sends an email to you when an R process has finished? For instance, I often work with very large data sets and certain tasks (e.g., merging records, lmer runs) can take a long time and I find myself constantly looking over to see if

Re: [R] navel-gazing

Ben, I change the line: z - getURL(paste(baseurl,list,/, month,/author.html,sep=)) to z - getURL(paste(baseurl,list,/, month,/author.html,sep=), ssl.verifypeer = FALSE) because don't work for me. Nice! On Tue, Aug 17, 2010 at 1:47 PM, Ben Bolker bbol...@gmail.com wrote: month -

Re: [R] R Send an Email

See sendmailR package. On Tue, Aug 17, 2010 at 11:54 AM, Doran, Harold hdo...@air.org wrote: Just out of curiosity, has anyone ever written a function that sends an email to you when an R process has finished? For instance, I often work with very large data sets and certain tasks (e.g.,

Re: [R] R Send an Email

On Tue, Aug 17, 2010 at 3:54 PM, Doran, Harold hdo...@air.org wrote: Just out of curiosity, has anyone ever written a function that sends an email to you when an R process has finished? For instance, I often work with very large data sets and certain tasks (e.g., merging records, lmer runs)

Re: [R] R Send an Email

Hi Harold, I have not used it recently, but it seems promising: http://cran.r-project.org/web/packages/sendmailR/index.html HTH, Jorge On Tue, Aug 17, 2010 at 10:54 AM, Doran, Harold wrote: Just out of curiosity, has anyone ever written a function that sends an email to you when an R

Re: [R] 2 not defined because of singularities appearing after introducing Fixed Effects

Sorry, please nevermind. It seems to have been an econometrical problem after all. (A variable consisting of 2 or more of the 150 Country-Dummies for the fixed effects causes perfect multicollinearity as well. So does a variable that differs over cross-sections but not over periods, seemingly. I

Re: [R] how to selection model by BIC

On Tue, 2010-08-17 at 17:12 +0100, Xin__ Li wrote: Hi All: the package MuMIn can be used to select the model based on AIC or AICc. The code is as follows: data(Cement) lm1 - lm(y ~ ., data = Cement) dd - dredge(lm1,rank=AIC) print(dd) If I want to select the model by BIC, what code

Re: [R] Weird differing results when using the Wilcoxon-test

Thanks for the hint. I tested this generic example and had the same behavior also in a special example, that can be found below. This example does not involve continuity correction but exibits the same unexpected behavior: GDS_example = function() { print(--- GDS EXAMPLE ---, quote=F)

Re: [R] how to paste a vector to expression?

On Aug 17, 2010, at 11:51 AM, Kay Cichini wrote: hello, i need to annotate a plot with an expression and a variable value - like: plot(1:4) dat-1:2 text(2:3,2:3,expression(paste(bar(x)==dat))) ... but with the dat values plotted, 1 at x=2,y=2, and 2 at x=3,y=3. One way ... actually

[R] Problems building own package (Error: package has been build before R-2.10.0)

Dear List, I’m doing my first baby steps towards developing own R Packages and ran into the following problem: R CMD check mypackage works fine (no errors, no warnings) R CMD build mypackage works fine (no errors, no warnings) R CMD INSTALL –library=”C:\R\R-2.11.1\library”

Re: [R] using color palette in persp

yes, I had, and failed somehow, because I couldn't work out some points. I don't have my work with me at the moment, but will give a detailed description about my porblems. I don't remember for sure what I tried and what not, but I failed at this line: zfacet - z[-1, -1] + z[-1, -ncz] + z[-nrz,

[R] Tables in Packages

I have a 100x2 matrix to go into a package. I could type it into the R code but, more generally, what is the *best* way to read a table from disk into, say, package 'foo' when it is loaded into memory via library(foo)? I thought of having .onLoad create an environment and assign a matrix read from

Re: [R] Missing values

- Original Message - From: Stephen Liu sati...@yahoo.com To: r-help@r-project.org Sent: Tuesday, August 17, 2010 8:25 AM Subject: Re: [R] Missing values - Original Message From: Michael Bedward michael.bedw...@gmail.com To: Stephen Liu sati...@yahoo.com Cc:

Re: [R] compare gam fits

On Thu, 2010-08-05 at 13:49 -0300, Mike Lawrence wrote: Hi folks, I've been pondering this for a little while, drumming up courage to reply, whilst I've been away on fieldwork with patchy email access... so here goes with the neck sticking out thing... ;-) I originally tried R-SIG-Mixed-Models

Re: [R] grid.table and expression in table body?

On 17 August 2010 18:38, Johannes Graumann johannes_graum...@web.de wrote: Hi - I can't get this figured out ... Thanks for any hint. The parse=TRUE argument means that all the table content will be parsed and interpreted as expression. In your content you have invalid expressions (e.g. text

Re: [R] navel-gazing

Since Peter Dalgaard is splitting his considerable contributions between Peter Dalgaard and peter dalgaard, I made the following changes (which shouldn't be a problem unless e e cummings becomes a regular poster): # from base::chartr documentation capwords - function(s, strict = FALSE) {

[R] Combining values from a data.frame row with matrix rownames.

Apologies, I neglected to check that the email was plain text. Here is a the text of the previous email. Hello, I am trying to iteratively build a build a panel of variables to discriminate between two groups. My starting position is a matrix of experimental data and I have a function that

Re: [R] TM Package - Corpus function - Memory Allocation Problems

I'm using R 2.11.1 on Win XP (32-bit) with 3 GB of RAM. My data has (only) 16.0 MB. I want to create a VCorpus object using the Corpus function in the tm package but I'm running into Memory allocation issues: Error: cannot allocate vector of size 372 Kb. My data is stored in a csv file which

Re: [R] TM Package - Corpus function - Memory Allocation Problems

On Aug 17, 2010, at 3:45 PM, Guelman, Leo wrote: I'm using R 2.11.1 on Win XP (32-bit) with 3 GB of RAM. My data has (only) 16.0 MB. Probably more than that. Each numeric is 8 bytes even before overhead, so a csv file that was all single digit integers and commas would more that double

Re: [R] navel-gazing

I think that gsub example on help page is more clear: library(XML) # could be used the XML package to get the names cnames - gsub('\n', '', head(tail(sapply(getNodeSet(htmlParse(z, asText = TRUE), //i), xmlValue), -3), -3)) gsub((\\w)(\\w*), \\U\\1\\L\\2, cnames, perl=TRUE) On Tue, Aug 17,

  1   2   >