Re: [R] Trouble with 'smooth' using xyplot in lattice

Peter, The error message comes from panel.loess, which is the panel function that draws the Loess smoothing line. It means one of your groups has too few data points or is too clustered in its domain for the Loess smoother to converge. You could - give up on the pretty-looking smoother -

Re: [R] Convert ASCII string to Hex in R (vice versa)

Dear all, For example I have the following ASCII strings: ascii_str - ORRRIROOQRQQOLORRQRROKKIKFAA ascii_str2 - FFDDC Each character in the above string represent a hexadecimal value. I want to translate those string into an array (of size == string length) which

Re: [R] Convert ASCII string to Hex in R (vice versa)

Gundala Viswanath gundalav at gmail.com writes: For example I have the following ASCII strings: ascii_str - ORRRIROOQRQQOLORRQRROKKIKFAA ascii_str2 - FFDDC Each character in the above string represent a hexadecimal value. I want to translate those string into

[R] Convert ASCII string to Decimal in R (vice versa) was: Hex

Hi Dieter, Sorry my mistake. I wanted to convert them into Decimal (not Hexadecimal). Given this string, the desired answer follows: ascii_str - ORQIK 79 82 81 62 73 75 ascii_str2 - FDC 70 68 67 - Gundala Viswanath Jakarta - Indonesia On Mon, Dec 22, 2008 at 5:49 PM, Dieter Menne

[R] R 2.8.1 is released

I've rolled up R-2.8.1.tar.gz a short while ago. This is a maintenance release and fixes a number of mostly minor bugs and platform issues. See the full list of changes below. You can get it (in a short while) from http://cran.r-project.org/src/base/R-2/R-2.8.1.tar.gz or wait for it to be

Re: [R] Convert ASCII string to Decimal in R (vice versa) was: Hex

?raw On Mon, Dec 22, 2008 at 4:22 AM, Gundala Viswanath gunda...@gmail.com wrote: Hi Dieter, Sorry my mistake. I wanted to convert them into Decimal (not Hexadecimal). Given this string, the desired answer follows: ascii_str - ORQIK 79 82 81 62 73 75 ascii_str2 - FDC 70 68 67 -

Re: [R] Convert ASCII string to Decimal in R (vice versa) was: Hex

Gundala Viswanath gundalav at gmail.com writes: Sorry my mistake. I wanted to convert them into Decimal (not Hexadecimal). Given this string, the desired answer follows: ascii_str - ORQIK 79 82 81 62 73 75 ascii_str2 - FDC 70 68 67 You might have a look at Henrik Bengtsson's

Re: [R] Convert ASCII string to Decimal in R (vice versa) was: Hex

Gundala Viswanath gundalav at gmail.com writes: Sorry my mistake. I wanted to convert them into Decimal (not Hexadecimal). Given this string, the desired answer follows: ascii_str - ORQIK 79 82 81 62 73 75 ascii_str2 - FDC 70 68 67 lapply(ABCD,charToRaw)[[1]] Dieter

[R] row sum question

Dear helpers, I'm using R version 2.8.0. Suppose that I have a small data set like below. [,1] [,2] [,3] [,4] a1100 b0110 c1110 d0111 First, I'd like to find row sum of values uniquely present in each row, but only

Re: [R] How can I get the interpolated data?

Thank you very much! It's very helpful to me! David Winsemius wrote: If you look at the CR.rsm object with str() you will see that it inherits from the lm class of objects. Therefore the predict.lm method will be available and if you further look at: ?predict.lm You see that all

Re: [R] error bars

On Fri, 2008-12-19 at 13:06 +, Kelly-Gerreyn B.A. wrote: Dear Help I'm new to R (from matlab)...using windows XP. I've tried to work out, to no avail, 4 things: 1) rotating the numbers on axes...something to do with par(str) ?? ?par and argument 'las' for basic control. There is a

Re: [R] sorting variable names containing digits

Note that mysort2 is slightly more general as it handles the case that the strings begin with numerics: u - c(51a2, 2a4) mysort(u) [1] 51a2 2a4 mysort2(u) [1] 2a4 51a2 On Mon, Dec 22, 2008 at 12:32 AM, John Fox j...@mcmaster.ca wrote: Dear Gabor, Thank you (again) for this second

Re: [R] row sum question

Will this do it for you: nrows - 10 ncols - 10 mat - matrix(sample(0:1, nrows * ncols, TRUE), nrow=nrows, ncol=ncols) mat [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,]101110000 0 [2,]000000011 0

[R] row sum question

Dear helpers, I'm using R version 2.8.0. Suppose that I have a small data set like below. [,1] [,2] [,3] [,4] a1100 b0110 c1110 d0111 First, I'd like to find sum of each row uniquely present in each row, but only sequentially

[R] Hmisc Dotplot with confidence intervals and panel.points problem

Hello useRs, I have a question regarding the function Dotplot from the Hmisc package: I want two things: 1) confidence intervals around the dots 2) some additional annotation points plotted in the graphic I can easily achieve (1) by constructing an appropriate object with Cbind. But for (2)

[R] AR(2) coefficient interpretation

I am a beginner in using R and I need help in the interpretation of AR result by R. I used 12 observations for my AR(2) model and it turned out the intercept showed 5.23 while first and second AR coefficients showed 0.40 and 0.46. It is because my raw data are in million so it seems the

Re: [R] row sum question

Assuming DF is a data frame like this: DF - data.frame(V1 = c(1, 0, 1, 0), V2 = c(1, 1, 1, 1), V3 = c(0, 1, 1, 1), V4 = c(0, 0, 0, 1)) # try this: head(rowSums((rbind(0, cummax(DF)) rbind(cummax(DF), 0))), -1) On Sun, Dec 21, 2008 at 8:36 PM, Keun-Hyung Choi khc...@sfsu.edu wrote: Dear

[R] imputing the numerical columns of a dataframe, returning the rest unchanged

Hi R-experts, how can I apply a function to each numeric column of a data frame and return the whole data frame with changes in numeric columns only? In my case I want to do a median imputation of the numeric columns and retain the other columns. My dataframe (DF) contains factors, characters and

[R] Matching

I understand this is an easy question, and have been playing around with grep and the match function, but was hoping for a little incite: I have one .csv with the following data: names values A1 B2 C3 D4 The second .csv is: names A C I am hoping to match all

Re: [R] Matching

Try this: Lines1 - names,values + A,1 + B,2 + C,3 + D,4 Lines2 - names + A + C DF1 - read.csv(textConnection(Lines1)) DF2 - read.csv(textConnection(Lines2)) merge(DF1, DF2) names values 1 A 1 2 C 3 On Mon, Dec 22, 2008 at 11:03 AM, vpas vic.pas...@gmail.com wrote: I

[R] post hoc comparisons on interaction means following lme

Dear Colleagues, I have scoured the help files and been unable to find an answer to my question. Please forgive me if I have missed something obvious. I have run the following two models, where category has 3 levels and comp has 8 levels: mod1 - lmer(x~category+comp+(1|id),data=impchiefsrm)

Re: [R] Matching

You don't need grep for this. Use the merge() function and make sure the arguments all.x and all.y are considered depending on whether this is a left or right merge. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of vpas Sent:

[R] Treatment of Date ODBC objects in R (RODBC)

Dear all, Retrieving an Oracle Date data type by means of RODBC (version 1.2-4) I get different classes in R depending on which operating system I am in: On MacOSX I get Date class On Windows I get POSIXt POSIXct class The problem is material, as converting the POSIXt POSIXct object

Re: [R] Treatment of Date ODBC objects in R (RODBC)

I tend to avoid the issue by asking Oracle for a character string representation of the date. I use sql like this: to_char( thedatefield, 'mmdd' ) as thedate Then in R: d - as.Date( as.character( thedate, '%Y%m%d') ) Hope this helps, Avram --Original Message-- From: Ivan

[R] package sn

This message is of interest only to users of package sn. In early 2007, I have posted an announcement on the package web site http://azzalini.stat.unipd.it/SN/announce2.html about a forthcoming version 1. This will be a deeply revised version, with substantial chances both in the internal

Re: [R] queue simulation

On Mon, 22 Dec 2008, Gerard M. Keogh wrote: Hi all, I have a multiple queing situation I'd like to simulate to get some idea of the distributions - waiting times and allocations etc. Does R has a package available for this - many years ago there used to be a language called simscript for

Re: [R] Globbing Files in R

Gundala Viswanath wrote: Typically Perl's idiom would be: __BEGIN__ @files = glob(/mydir/*.txt); foreach my $file (@files) { # process the file } __END__ Something like this has been suggested in R-help before: files - dir() results - lapply(files, yourprocessing()) The dir function

Re: [R] Globbing Files in R

On 12/22/2008 1:14 PM, Earl F Glynn wrote: Gundala Viswanath wrote: Typically Perl's idiom would be: __BEGIN__ @files = glob(/mydir/*.txt); foreach my $file (@files) { # process the file } __END__ Something like this has been suggested in R-help before: files - dir() results -

[R] questions about read datafile into R

Dear all: I have been thinking to import below one data file (.txt)into R by read.table(..,skip=1, header=T). But How can I deal with the repeated rows of TABLE NO.1 and names of data variables in the middle of this data file. The similar block will be repeated 100 times, here only show 4 of

Re: [R] questions about read datafile into R

This kind of problem - reading a messy data file - was discussed at great length just last week. The ideas from that discussion may help you with your version: http://tolstoy.newcastle.edu.au/R/e5/help/08/12/10404.html Sarah On Mon, Dec 22, 2008 at 2:38 PM, Lu, Zheng zheng...@mpi.com wrote:

Re: [R] How can I avoid nested 'for' loops or quicken the process?

On Mon, 22 Dec 2008, Brigid Mooney wrote: Hi All, I'm still pretty new to using R - and I was hoping I might be able to get some advice as to how to use 'apply' or a similar function instead of using nested for loops. Unfortunately, you have given nothing that is reproducible. The details

[R] sem package fails when no of factors increase from 3 to 4

I checked through every 3 factor * 3 loading case. While, 4 factor * 3 loading failed. the data is 6 factor * 3 loading require(sem); cor18-read.moments(); 1 .68 1 .60 .58 1 .01 .10 .07 1 .12 .04 .06 .29 1 .06 .06 .01 .35 .24 1 .09 .13 .10 .05 .03 .07 1 .04 .08 .16 .10 .12 .06 .25

Re: [R] queue simulation

Date: 22-Dec-2008 10:11:28 GMT From: Gerard M. Keogh gmke...@justice.ie Subject: [R] queue simulation To: r-help@r-project.org Hi all, I have a multiple queing situation I'd like to simulate to get some idea of the distributions - waiting times and allocations etc. Does R has a

[R] How can I avoid nested 'for' loops or quicken the process?

Hi All, I'm still pretty new to using R - and I was hoping I might be able to get some advice as to how to use 'apply' or a similar function instead of using nested for loops. Right now I have a script which uses nested for loops similar to this: i - 1 for(a in Alpha) { for (b in Beta) { for (c

[R] offlist Re: How can I avoid nested 'for' loops or quicken the process?

I do agree with Dr Berry that your question failed on several grounds in adherence to the Posting Guide, so this is off list. Maybe this will give you guidance that you can apply to your next question to the list: alist - list(a,b,c) blist - list(ab,ac,ad) expand.grid(alist, blist)

Re: [R] Error compiling R.2.8.1 with gcc 4.4 on Mac OS 10.5.6

Mike Lawrence wrote: Hi all, I've encountered a build error with the latest R source (2.8.1). This is a relatively fresh install of OS Leopard (10.5.6), latest developer tools installed, gcc/g++/gfortran version 4.4 installed (via http://hpc.sourceforge.net/, after which I updated the gcc g++

[R] Summary information by groups programming assitance

All - I have data that looks like psd Species Lake Length WeightSt.weightWr Wr.1 vol 432 substock SMB Clear150 41.00 0.01 95.12438 95.10118 0.0105 433 substock SMB Clear152 39.00 0.01 86.72916 86.70692 0.0105 434 substock

Re: [R] Summary information by groups programming assitance

On Mon, Dec 22, 2008 at 3:51 PM, Ranney, Steven steven.ran...@montana.edu wrote: All - I have data that looks like psd Species Lake Length WeightSt.weightWr Wr.1 vol 432 substock SMB Clear150 41.00 0.01 95.12438 95.10118 0.0105 433 substock

Re: [R] Summary information by groups programming assitance

Maybe summaryBy (or lapplyBy/splitBy) in the doBy package might help you. Regards Søren Fra: r-help-boun...@r-project.org på vegne af Ranney, Steven Sendt: ma 22-12-2008 22:51 Til: r-help@r-project.org Emne: [R] Summary information by groups programming assitance

Re: [R] Error compiling R.2.8.1 with gcc 4.4 on Mac OS 10.5.6

On Mon, 22 Dec 2008, Peter Dalgaard wrote: Mike Lawrence wrote: Hi all, I've encountered a build error with the latest R source (2.8.1). This is a relatively fresh install of OS Leopard (10.5.6), latest developer tools installed, gcc/g++/gfortran version 4.4 installed (via

[R] newbie question on tcltk

Hi List, Can anyone tell me how could i put the BACK button in the following code, just under the AAA menu? I want this button to go back to the previous page, and since it has nothing to do with the 1 and 2 buttons, i want it somehow separated from these two buttons, but i don't know how. I

Re: [R] newbie question on tcltk

*** WARNING *** Just a warning to anyone thinking of copying the code below into their workspace -- the second line in the code below will erase your entire workspace. On Mon, Dec 22, 2008 at 5:57 PM, eugen pircalabelu eugen_pircalab...@yahoo.com wrote: Hi List, Can anyone tell me how could i

[R] Error: cannot allocate vector of size 1.8 Gb

dim(data) [1] 2228319 dm=dist(data, method = euclidean, diag = FALSE, upper = FALSE, p = 2) Error: cannot allocate vector of size 1.8 Gb Hi Guys, thank you in advance for helping. :-D Recently I ran into the cannot allocate vector of size 1.8GB error. I am pretty sure this is not a

Re: [R] Integrate function

If these messages you're hearing are warnings, then the answer might be: ?warnings -- David Winsemius On Dec 22, 2008, at 6:07 PM, glenn roberts wrote: Quick One if any one can help please. On use of integration function ‘integrate’; how do I get the function to return just the value

Re: [R] Summary information by groups programming assitance

Yet another suggestion is describe.by in the psych package. At 11:25 PM +0100 12/22/08, Søren Højsgaard wrote: Maybe summaryBy (or lapplyBy/splitBy) in the doBy package might help you. Regards Søren Fra: r-help-boun...@r-project.org på vegne af Ranney, Steven

Re: [R] Error: cannot allocate vector of size 1.8 Gb

I hope the following info will help, thanks again! sessionInfo() R version 2.7.1 (2008-06-23) x86_64-redhat-linux-gnu locale:

Re: [R] Integrate function

.. but it turned out he wanted; integrate(integrand)$value -- David Winsemius On Dec 22, 2008, at 6:26 PM, David Winsemius wrote: If these messages you're hearing are warnings, then the answer might be: ?warnings -- David Winsemius On Dec 22, 2008, at 6:07 PM, glenn roberts wrote:

Re: [R] Summary information by groups programming assitance

Here are two solutions assuming DF is your data frame: # 1. aggregate is in the base of R aggregate(DF[c(Length, vol)], DF[c(Lake, psd)], max) or the following which is the same except it labels psd as Category: aggregate(DF[c(Length, vol)], with(DF, list(Lake = Lake, Category = psd)), max)

Re: [R] sem package fails when no of factors increase from 3 to 4

Dear Xiaoxu LI, sem.mod(mod4, cor18, 500, debug=TRUE) will show you what went wrong with the optimization. Since the three-factor solutions look reasonable, I tried using them to get better start values for the parameters in the four-factor model, producing the solution shown below. As well, I

Re: [R] svyglm and sandwich estimator of variance

On Fri, 19 Dec 2008, Roberta Pereira Niquini wrote: Hi, I would like to estimate coefficients using poisson regression and then get standard errors that are adjusted for heteroskedasticity, using a complex sample survey data. Then I will calculate prevalence ratio and confidence intervals. Can

Re: [R] question about read datafile

Read in the data using readLines to read the complete line. Use grep/regexpr to scan for valid lines and then convert them to numeric by using strsplit/as.numeric. On Mon, Dec 22, 2008 at 2:13 PM, Lu, Zheng zheng...@mpi.com wrote: Dear all: I have been thinking to import below one data file

Re: [R] Summary information by groups programming assitance

Thank you all for your help. I appreciate the assistance. I'm thinking I should have been more specific in my original question. Unless I'm mistaken, all of the suggestions so far have been for maximum vol and maximum Length by Lake and psd. I'm trying to extract the max vol by Lake and

[R] Problem in passing on an argument via ... how do I access it?

Hi r-experts, I want to check if a certain argument has been passed on in a function call via ... ftest - function(x1, ...) { if(hasArg(y2)==TRUE) print(y2) } Now I call the function passing y2 via ... but I cannot access or use the object. ftest(y2= 2, x= 1) error in print(y2) : object

[R] Simulate dataset using Parallel Latent CTT model in R

All, I want to simulate dataset using Parallel Latent CTT model in R however dont know how to start. Is there anyone who have done the same? Any help on this will be greatly appreciated. Regards -NK -- View this message in context:

[R] nlsrob fails with puzzling error message on input accepted by nls

I have a nonlinear model estimation problem with ~50,000 data records and a simple 3 parameter model (logistic type - please don't tell me that there are linear methods for such a problem). I run nls with constraints once to get a good initial parameter guess, then try to run nlrob to get

Re: [R] Summary information by groups programming assitance

On Mon, Dec 22, 2008 at 7:15 PM, Ranney, Steven steven.ran...@montana.edu wrote: Thank you all for your help. I appreciate the assistance. I'm thinking I should have been more specific in my original question. Unless I'm mistaken, all of the suggestions so far have been for maximum vol and

Re: [R] Summary information by groups programming assitance

Just sort the data first and then apply any of the solutions but with tail(x, 1) instead of max, e.g. DFo - DF[order(DF$Lake, DF$Length, DF$vol), ] aggregate(DFo[c(Length, vol)], DFo[c(Lake, psd)], tail, 1) On Mon, Dec 22, 2008 at 8:15 PM, Ranney, Steven steven.ran...@montana.edu wrote: Thank

Re: [R] Summary information by groups programming assitance

The sorting should have been by Lake, psd and vol (not what I had) so it should be revised to: DFo - DF[order(DF$Lake, DF$psd, DF$vol), ] aggregate(DFo[c(Length, vol)], DFo[c(Lake, psd)], tail, 1) This is the same as before except DF$psd is used in place of DF$Length in the first line. On Mon,

Re: [R] Problem in passing on an argument via ... how do I access it?

Try: ftest - function(x1, ...) { yargs =list(...) ; if (hasArg(y2) == TRUE) print(YES); return(yargs) } ftest(2, y2 = 3) [1] YES $y2 [1] 3 yt - ftest(2, y2=3) [1] YES yt $y2 [1] 3 On Dec 22, 2008, at 6:16 PM, Mark Heckmann wrote: Hi r-experts, I want to check if a certain

[R] sorting regression coefficients by p-value

Hi, Is there a way to get/extract a matrix of regression variable name, coefficient, and p values? (for lm and glm; which can be sort by p value?) thanks Dhruv [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] Tabular output: from R to Excel or HTML

What is the equivalent for formatted tabular output of the various very sophisticated plotting tools in R (plot, lattice, ggplot2)? In particular, I'd like to be able to produce formatted Excel spreadsheets (using color, fonts, borders, etc. -- probably via Excel XML) and formatted HTML tables

Re: [R] sorting regression coefficients by p-value

Assuming that you are using the example in the lm help page: ctl - c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) trt - c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69) group - gl(2,10,20, labels=c(Ctl,Trt)) weight - c(ctl, trt) lm.D9 - lm(weight ~ group) # The coefficients are just :

Re: [R] Tabular output: from R to Excel or HTML

You should be looking at odfWeave. It has support for the OpenOffice table formatting and once those are created, the conversion to Excel should proceed smoothly. -- David Winsemius On Dec 22, 2008, at 11:58 PM, Stavros Macrakis wrote: What is the equivalent for formatted tabular output

Re: [R] AR(2) coefficient interpretation

As I need your urgent help so let me modify my question. I imported the following data set to R and run the statements i mentioned in my previous reply Year Month Period ab c 1 2008 Jan 2008-Jan 105,536,785 9,322,074 9,212,111 2 2008 Feb 2008-Feb 137,239,037 10,986,047

[R] newbie problem using Design.rcs

Hi, I read data from a file. I'm trying to understand how to use Design.rcs by using simple test data first. I use 1000 integer values (1,...,1000) for x (the predictor) with some noise (x+.02*x) and I set the response variable y=x. Then, I try rcs and ols as follows: m = ( sqrt(y1) ~ (

Re: [R] QCA adn Fuzzy

Dear Gott and Prof Adrian DUSA , I am learning fuzzy set QCA and recently, I just write a function to construct a truthTable, which can be passed to QCA:::eqmcc to do the Boolean minimization. The function is here: http://code.google.com/p/asrr/source/browse/trunk/R/fs_truthTable.R and the help