[R] readOGR workaround for Japanese UTF-8 geojson

I am working on a project https://github.com/AlanInTsukuba/jpucd that involves extracting shapefiles and property data from Japanese geojson files. When reading with readOGR(ibarakipath1 , encoding="UTF-8", use_iconv=TRUE), I find that the subsets of cannot be written with writeOGR without losing t

Re: [R] [FORGED] Newbie Question on R versus Matlab/Octave versus C

On 1/29/2019 11:50 PM, Jeff Newmiller wrote: On Tue, 29 Jan 2019, Alan Feuerbacher wrote: After my failed attempt at using Octave, I realized that most likely the main contributing factor was that I was not able to figure out an efficient data structure to model one person. But C lent itself

Re: [R] Newbie Question on R versus Matlab/Octave versus C

On 1/30/2019 7:44 AM, Martin Møller Skarbiniks Pedersen wrote: On Mon, 28 Jan 2019 at 22:58, Alan Feuerbacher <mailto:alan...@comcast.net>> wrote: [...] > These all have advantages and > disadvantages for certain tasks, but as I'm new to R I hardly know how >

Re: [R] [FORGED] Newbie Question on R versus Matlab/Octave versus C

On 1/29/2019 11:50 PM, Jeff Newmiller wrote: Thanks very much for providing these coding examples! I think this is a good way to learn some R. Alan On Tue, 29 Jan 2019, Alan Feuerbacher wrote: On 1/28/2019 7:51 PM, Jeff Newmiller wrote: If you forge on with your preconceptions of how such

Re: [R] how to ref data directly from bls.gov using their api

On 1/29/2019 11:24 AM, Bert Gunter wrote: Please search on "Bureau of Labor Statistics" at rseek.org. You will find several packages and other resources there for doing what you want. Wow! Thanks! Again most likely severe overkill for my simple need. Alan -- Ber Bert Gunter &q

Re: [R] Newbie Question on R versus Matlab/Octave versus C

and is particularly suitable in terms of performance for iterative computations where one iteration depends on the prior one. Given my response to Jeff Newmiller, do your comments still apply? Alan On Mon, Jan 28, 2019 at 6:32 PM Gabor Grothendieck wrote: R has many similarities to Octave

Re: [R] [FORGED] Newbie Question on R versus Matlab/Octave versus C

ix C++ with R quite easily and then you get the best of both worlds. (C and Fortran are supported, but they are a bit more finicky to setup than C++). I don't know the answer to that, but perhaps you can help decide. Alan On January 28, 2019 4:00:07 PM PST, Alan Feuerbacher wrote: On

Re: [R] [FORGED] Newbie Question on R versus Matlab/Octave versus C

ation fun. At this point I think you're right on all counts. Alan Bill Dunlap TIBCO Software wdunlap tibco.com <http://tibco.com> On Mon, Jan 28, 2019 at 4:00 PM Alan Feuerbacher <mailto:alan...@comcast.net>> wrote: On 1/28/2019 4:20 PM, Rolf Turner wrote: >

Re: [R] [FORGED] Newbie Question on R versus Matlab/Octave versus C

ore complicated than my simple problem, and has a large learning curve. Alan Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )

Re: [R] [FORGED] Newbie Question on R versus Matlab/Octave versus C

On 1/28/2019 4:20 PM, Rolf Turner wrote: On 1/29/19 10:05 AM, Alan Feuerbacher wrote: Hi, I recently learned of the existence of R through a physicist friend who uses it in his research. I've used Octave for a decade, and C for 35 years, but would like to learn R. These all have advan

[R] Newbie Question on R versus Matlab/Octave versus C

Hi, I recently learned of the existence of R through a physicist friend who uses it in his research. I've used Octave for a decade, and C for 35 years, but would like to learn R. These all have advantages and disadvantages for certain tasks, but as I'm new to R I hardly know how to evaluate t

Re: [R] Passing a Data Frame Name as a Variable in a Function

Much thanks, Hadley! Cheers, Alan On Jan 29, 2015, at 12:36 PM, Hadley Wickham wrote: > On Thu, Jan 29, 2015 at 11:43 AM, Alan Yong wrote: >> Much thanks to everyone for their recommendations! I agree that fishing in >> the global environment isn't ideal & only shows

Re: [R] Passing a Data Frame Name as a Variable in a Function

en, (2) Bind the rows of all the files. Any more feedback will be much appreciated! Cheers, Alan On Jan 29, 2015, at 3:11 AM, peter dalgaard wrote: > > On 29 Jan 2015, at 07:34 , Jeff Newmiller wrote: > >> This approach is fraught with dangers. >> >> I recommend that yo

[R] Passing a Data Frame Name as a Variable in a Function

1]) however, when I run: > length(DFName[, 1]) Error in DFName[, 1] : incorrect number of dimensions and length(df.1001[, 1]) [1] 104 do not provide the same expected answer. How can I successfully pass the data frame name of df.1001 as a variable named DFName in a

Re: [R] Symbolic equations to R code?

Have you thought about programming a function builder into the Shiny applet? It might simplify the process for your students a bit. A highly simplified version of the page David suggested could be built that generates the R code. Best of luck, Alan -Original Message- From: Scott

[R] kripp.alpha - Getting Data into It

I have used read.csv put three columns (of raters) who rated 188 things (rows) into the object PublicKnow After loading irr and calling kripp.alpha(PublicKnow), I get what I see on searching is the common error: Error in sort.list(y) : 'x' must be atomic for 'sort.list' Have you called 'sort' o

[R] Unexpected input error

Hi, I am newbie to R so pardon me if this is not the forum for asking this question. I am getting the error "unexpected input in "evalq({‚" in the Boost Tutorial example section 13.3 of Graham William's, "Data Mining with Rattle and R. Below is the line as I have entered it and version i

Re: [R] Advice on use of R for Generalised Linear Modelling

Interesting...interesting... Thank you very much for this. I'll play with this tomorrow & see what happens. Regards, Alan. On 12 August 2013 17:15, Yvonnick Noel wrote: > Alan, > > You might want to have a look to the R2STATS package on CRAN. It is a GUI > for GLM an

[R] Advice on use of R for Generalised Linear Modelling

Hi, Not an expert R user, something of a novice - please be gentle with me! I have a particular interest in generalised linear models (GLMs) and I'm experienced in fitting them using other bits of software. R can fit GLMs of course, using the glm() command. I have some large multivariate data s

Re: [R] lmer( ) parameter estimates changing depending on dummy-coded reference level

Thank you, I'll post there. Regards, Alan On Wed, Aug 7, 2013 at 1:14 PM, Bert Gunter wrote: > I think it's fair to say that this should really be posted on the > mixed-model specific (especially using lme4) list r-sig-mixed-models > and not here. > > Cheers, > Ber

[R] lmer( ) parameter estimates changing depending on dummy-coded reference level

respond correctly. However, I didn't get any warnings about model convergence, or any warnings at all.) Any clues as to why I'm getting this results would be very much appreciated. Thanks in advance, Alan Mishler Research Assistant University of Maryland -- ## Model 1 output: Blo

Re: [R] Looking for a test of standard normality

Shapiro-Wilk normality test data: qq W = 0.9831, p-value = 0.2286 Alan Herschtal Senior Biostatistician Peter MacCallum Cancer Centre Phone +61 3 9656 3639 Fax +61 3 9656 1420 Email alan.hersch...@petermac.org From: Mark Lamias

Re: [R] Looking for a test of standard normality

the percentiles, checking for equal numbers in each decile bin. This would certainly work, and the only disadvantage that I can see is that the selection of the bin boundaries is somewhat arbitrary. Alan Herschtal Senior Biostatistician Peter MacCallum Cancer Centre Phone +61 3 9656 3639 Fax +61 3

[R] Looking for a test of standard normality

that the data are not standard normal. My back up option is to use a Kolmogorov-Smirnov test, but my impression is that that is not a very powerful test. Any suggestions? Thanks, Alan Herschtal Senior Biostatistician Peter MacCallum Cancer Centre Phone +61 3 9656 3639 Fax +61 3 9656 1420 Email

[R] AUTO: Alan Chalk has left RSA (returning 30/11/2012)

I am out of the office until 30/11/2012. Please send work related emails to laura.jor...@uk.rsagroup.com or personal emails to alanch...@gmail.com. Note: This is an automated response to your message "R-help Digest, Vol 117, Issue 3" sent on 03/11/2012 11:00:07. This is the only notification y

[R] ggplot2 with separate average lines

(Type, ID, Size), aes(x=ID, y=Size, shape=Type, color=Type)) + geom_point(shape=1) + scale_fill_hue(l=40) + geom_line(stat="hline", yintercept=mean, linetype="dashed") > Regards, Alan [[alternative HTML version deleted]]

Re: [R] Force evaluation of a symbol when a function is created

Thank you both, this was very helpful. I need to study environments more. Do either of you know a good source? -Original Message- From: Bert Gunter [mailto:gunter.ber...@gene.com] Sent: Monday, August 06, 2012 6:03 PM To: William Dunlap Cc: Schoenfeld, David Alan,Ph.D.,Biostatistics

[R] Force evaluation of a symbol when a function is created

I am porting a program in matlab to R, The problem is that Matlab has a feature where symbols that aren't arguments are evaluated immediately. That is: Y=3 F=@(x) x*Y Will yield a function such that F(2)=6. If later say. Y=4 then F(2) will still equal 6. R on the other hand has lazy evaluation.

[R] Memory limit for Windows 64bit build of R

Alan Alan Simpson Technical Lead, Retail Model Development Retail Models Project National Australia Bank Level 15, 500 Bourke St, Melbourne VIC Tel: +61 (0) 3 8697 7135 | Mob: +61 (0) 412 975 955 Email: alan.x.simp...@nab.com.au The information contained in this email and its attachments may be

Re: [R] How to test for the difference of means in population, please help

try anova if they have same observe points. conditions<-rep(c(1,2),each=1600) response<-sample (1:20,1600*2, replace= TRUE) points<-rep(rep(1:400, each=4),2) replicates<-rep(1:4,800) obs<-data.frame(response,conditions,points,replicates) fit <- lm(response~conditions+points,data=obs) anova(fit)

[R] plot points using circles filled half in red and half in blue.

I want to plot many points and want to use circles. The filling color depends on variable a. if a=1, then not fill if a=2 then fill with red, if a=3 then fill with blue, if a=4, fill half with red and half with blue. Can anyone tell me how to plot the case "a=4"? Thanks a lot _

Re: [R] Actual vs. predicted survival times

Have you tried plotting obvserved survival against X*Beta? I believe the usual predictions from a cox model are just monotonic transformations of this. -Alan -Original Message- From: Bonnett, Laura [mailto:l.j.bonn...@liverpool.ac.uk] Sent: Wed 2/8/2012 1:52 AM To: 'David Winsemiu

[R] Looping and actual variable names

eader=T) for (j in 2:5) print(summary(lm(formula = dd[,1] ~ dd[,j]))) However, although it produces all the analyses that I want, the output is annotated with dd[,1] and dd[,j] rather than the actual variable names. Is there any way around this? Best wishes, Alan. __

[R] scatterplotting stock returns using quantmod and pairs()

I want to get data for a set of ticker symbols and compute the daily return of the adjusted close using quantmod, and then scatterplot returns using pairs(). The following gets data for the list of tickers: tickers <- c("SHY","TLT","SPY","IWM","GLD","IEV","ILF","EWJ","EPP","SAF","ASA") AdjClosePr

[R] take me off the list

thank you. Alan Gao University of Michigan '12 Statistics, B.S. 239.682.3509 __ 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 pr

Re: [R] How to use IML with R and SAS

I will respond off list since the issue is mostly related to SAS. Alan Mitchell, MSc Biostatistician al...@crab.org -Original Message- From: Cheryl Johnson [mailto:johnson.cheryl...@gmail.com] Sent: Sunday, October 30, 2011 5:02 PM To: r-help@r-project.org Subject: [R] How to use IML

Re: [R] How to test if two C statistics are significantly different?

?Hmisc::rcorrp.cens -Alan -Original Message- From: Eik Vettorazzi [mailto:e.vettora...@uke.de] Sent: Tue 10/11/2011 2:25 AM To: Yujie Wang Cc: r-help@r-project.org Subject: Re: [R] How to test if two C statistics are significantly different? Hi Yujie, there is still a lot of work in

[R] How to get the last modified time of a file from R

Hi All, I want to add the last modified time of the dataset I used in the analysis. Any idea about how to get it from R in windows XP system? Thanks Alan __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] issue with image plot in 2.13.1

sized. This occurs using a range of graphic devices (windows, png, bmp, jpeg, tiff, pdf), on two different machines (windows 7-64 bit and windows XP-32 bit) but not when using 2.13.0 on either machine. Do others have this issue? Any possible solutions? Regards, Alan x <- y <- seq(-4*pi, 4

[R] dlmSum(...) and non-constant state space models

ng to find an answer. Please find my code below, perhaps there is a mistake which can easily be rectified. Kind regards, Alan > library(dlm) > mydata <- read.csv("Data.csv", header=TRUE)> pc <- mydata[,c(2,3,5,22)]> > > pc[,1] <- log(pc[,1]/pc[,3]) ; pc[,2] <

Re: [R] Plotting problems directional or rose plots

e(x=cos(t),y=sin(t)) par(mar=c(0,0,0,0)) plot(circle$x,circle$y,type='l') for(idx in 1:length(dist)){ t0 = seq(angle[idx],angle[idx]+10,length.out=50)*pi/180 petal0 = data.frame(x = dist[idx]*cos(t0),y = dist[idx]*sin(t0)) polygon( x=c(0,petal0$x,0),y=c(0,petal0$y,0),col=color.list

Re: [R] Plotting problems directional or rose plots

ist[1]*cos(t0) petal0$y = dist[1]*sin(t0) polygon( x=c(0,petal0$x,0),y=c(0,petal0$y,0),col=color.list[1]) I hope this helps. Alan Mitchell -Original Message- From: kitty [mailto:kitty.a1...@gmail.com] Sent: Tue 7/26/2011 2:20 PM To: r-help@r-project.org Subject: [R] Plotting probl

Re: [R] cumulative incidence plot vs survival plot

ent ways to do this, but the crr function in cmprsk will perform a competing risks regression. Alan Mitchell, MSc Ph: (206) 839-1708 al...@crab.org __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] cumulative incidence plot vs survival plot

cmprsk package is a great place to start. Gooley has a great paper on the difference between CI and 1-KM (See Statistics in Medicine, 18, 695-706 (1999)). HTH, Alan Mitchell, MSc Biostatistician al...@crab.org -Original Message- From: array chip [mailto:arrayprof...@yahoo.com] Sent

[R] AUTO: Chalk, Alan has a new email address

I am out of the office until 01/05/2012. Please note that I have changed email address. You can contact me at alan.ch...@uk.rsagroup.com Note: This is an automated response to your message "R-help Digest, Vol 100, Issue 26 [Virus checked]" sent on 25/06/2011 11:00:11. This is the only notif

[R] Conditional distribution plot using Model-based Recursive Partitioning

a way which I would be able to plot the tree and have either a conditional spinogram in the leaf or even just text with a number showing a marginal effect or coefficient please? Thanks, Alan -- Alan Fernihough IRCHSS Scholar UCD School of Economics [[alternati

Re: [R] Lists of tables and conditional statements

Thanks Henrique - that worked like a charm - I had tried lots of other combinations before seeing your reply - wished I had asked sooner! Alan -Original Message- From: David Winsemius [mailto:dwinsem...@comcast.net] Sent: Wednesday, March 30, 2011 10:09 PM To: Herbert, Alan G Cc: r

[R] Lists of tables and conditional statements

Hi R-users, I have a list containing numeric tables of differing row length. I want to make a new list that contains only rows from tables with a "Sum" greater than 3, plus the names of each table. I was wondering whether there is an elegant way to do this using apply of related functions as th

[R] No Support for plot=F in plclust()

Hello all, I am working on converting a set of S+ functions to R. Can anyone tell me how to extract data from an hclust or dendrogram object in R that is similar to that generated by P1 = plclust(tree, plot=F) Any suggestions would be appreciated. Alan Mitchell, MSc

Re: [R] MCMCpack combining chains

Thank you Uwe for the clarification. Alan Kelly On 24.02.2011 17:48, Alan Kelly wrote: > Deal all, as MCMClogit does not allow for the specification of several > chains, I have run my model 3 times with different random number seeds and > differently dispersed multivariate normal prio

[R] MCMCpack combining chains

chains So, how may I combine the separate mcmc objects so that I have multiple chains? Many thanks, Alan Kelly [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] Waaaayy off topic...Statistical methods, pub bias, scientific validity

Bert, consider the short rebuttal offered by George Musser in Scientific American: http://www.scientificamerican.com/blog/post.cfm?id=in-praise-of-scientific-error-2010-12-20 Perhaps a more realistic assessment of the (acknowledged) problem. Regards, Alan Kelly Trinity College Dublin On 7 Jan

[R] : unusual combinations of categorical data

Regarding unusual combinations of factors in categorical data. Are there any R packages that can be used to identify the outliers i.e. unusual combinations in categorical datasets ? Thanks. Notice of Confidentia

[R] re. Mathematica and R

David - information on calling R from within Mathematica can be found at the following link: http://www.mofeel.net/1164-comp-soft-sys-math-mathematica/13022.aspx HTH, Alan Kelly __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo

Re: [R] Avoiding Loops When Iterating Over Statement That Updates Its Input

Is there a performance advantage to doing this, as opposed to growing the vector within the loop? I suppose R could have to dynamically reallocate memory at some point? Alan 2010/5/30 Uwe Ligges : > > > On 26.05.2010 08:52, Alan Lue wrote: >> >> Come to think of it, we can

Re: [R] Data Frame as Hash Table

Thanks, guys! Alan On Sun, May 30, 2010 at 5:35 AM, Marshall Feldman wrote: > Besides data.table, there's the hash package. It does not use data.frame > type structures but is a bit more flexible. > > Marsh Feldman > > On 5/30/10 [May 30, 10] 6:00 AM, r-help-requ

[R] Data Frame as Hash Table

0.002720434 Then I'd like to be able to quickly retrieve the "value" of "key" 1.5 to get -0.53. How would one go about doing this? Yours, Alan Lue __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help P

Re: [R] Counting indexes

x <- rle(id) cumsum(x$lengths) - (x$lengths - 1) cumsum(x$lengths) Alan On Tue, May 25, 2010 at 10:00 PM, Gabor Grothendieck wrote: > This assumes that for a particular id they all occur together in a run: > >> cbind(start = which(!duplicated(id)), end = which(!duplicat

Re: [R] Avoiding Loops When Iterating Over Statement That Updates Its Input

Come to think of it, we can't save the output of each invocation and concatenate it later, since we need the output as input for the next iteration. Alan On Tue, May 25, 2010 at 11:43 PM, Alan Lue wrote: > Since `for' loops are slow in R, and since `apply' functions a

[R] Avoiding Loops When Iterating Over Statement That Updates Its Input

' and concatenate it later, but is there a better way to write this (i.e., to execute more quickly) while updating `r.seq' in each iteration? Alan __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do rea

Re: [R] Determining Index of Last Element in Vector

Those are great solutions. Thanks so much for your help. Yours, Alan On Sun, Apr 25, 2010 at 11:43 AM, David Winsemius wrote: > > On Apr 25, 2010, at 2:22 PM, Chuck Cleland wrote: > >> On 4/25/2010 2:10 PM, Alan Lue wrote: >>> >>> Hi, >>> >>&

Re: [R] Determining Index of Last Element in Vector

Sorry -- I meant `v(end)' and `v[length(v)]' in the first examples of my message. Alan On Sun, Apr 25, 2010 at 11:10 AM, Alan Lue wrote: > Hi, > Is there a way to specify the last element of a vector, similar to "end" in > MATLAB? >   v[end] > would be

[R] Determining Index of Last Element in Vector

r the first and last? (I know you can achieve this using two lines, but I'm writing because I want to do it using one.) Alan -- Alan Lue Master of Financial Engineering UCLA Anderson School of Management [[alternative HTML version deleted]] _

[R] nls.lm & AIC

lihood. Does someone have any idea of how I could proceed to get such informations about my models? Thanks for your help. Best regards, Alan Baudron The University of Aberdeen is a charity registered in Scotland, No SC013683. [[alternative

[R] Strange results from Windows 7

ot;s" } However, if I run it in another PC that uses Windows 7, the following error occurs: > source("c:\\Rwork\\test.txt") Read 290 items Error in diff[i] = "b" : object of type 'closure' is not subsettable I am using R 2.9.2

Re: [R] question about availability of splines library for 2.10.0

Phil and Peter, thanks - figured it was something obvious. I had uninstalled splines to check the behavior of some other code. I just reinstalled R 2.10.0 and it shows up. Thanks! Alan Peter Ehlers wrote: Have you searched your computer? Rhome/library/splines It's one of the pac

[R] question about availability of splines library for 2.10.0

.10.0 according to the description) on another machine. Where can I find this library? Thanks in advance, Alan > sessionInfo() R version 2.10.0 (2009-10-26) i386-pc-mingw32 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinf

[R] non-intuitive behaviour after type conversion

ext books on R. Can anyone comment on why this should occur? Many thanks, Alan Kelly Dr. Alan Kelly Department of Public Health & Primary Care Trinity College Dublin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-he

[R] Rounding and printing

="") [1] 56.5 1.6 0.2 3.9 0.1 2.2 0.2 2.6 1.5 4.1 1.1 6.1 21.0 [1] "56.5 (100%)" "1.6 (2.8%)" "0.2 (0.4%)" "3.9 (6.9%)" "0.1 (0.2%)" "2.2 (3.9%)" [7] "0.2 (0.4%)" "2.6 (4.6%)" "1.5 (2.

[R] Drawing lines in margins

Hi all, Quick question: What function can I use to draw a line in the margin of a plot? segments() and lines() both stop at the margin. In case the answer depends on exactly what I'm trying to do, see below. I'm using R v. 2.8.1 on Windows XP. Cheers, Alan I'm trying to ma

[R] Long to wide format without time variable

ssor-efficient way to replicate "reshape" without a time variable. Thanks in advance for any help you can provide. A worked example and some code I've tried are below. I'm working with R v2.8.1 on Windows XP Professional. Cheers, Alan Cohen Here's what my data look like

[R] nested cross-sectional design using lmer or nlme

random effects for unit and timecat X unit. I would much appreciate any suggestions on how to code the above in lmer or nlme. Alan Kelly Trinity College Dublin __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-hel

[R] help installing Rmpi

have included the processes I attempted below. Thank you, Alan Smith # what I have attempted so far #I have tried two different ways to install openMPI in fedora in the linux shell #in the linux shell cd Download wget http://www.open-mpi.org/software/ompi/v1.3

[R] Axis Limits in Scatterplot3d

Hi, How do you obtain the limits of the plotting region in a scatterplot3d plot? `par('usr')' does not seem to give sensible values, and that vector only has 4 elements (not the expected 6). Alan __ R-help@r-project.org mai

[R] Mixed Models

Hello, I've run a Proc Mixed function on a set of data in SAS. The data was a result of an experiment that measured % viability over time and I wanted to compare a Large sample lets say 50L to a small sample say 5L. And compare the % viability between the 2 sizes to see if I get the same answ

[R] Two-way Anova

Hello, I'm trying to do a comparsion on a large scale say 10L bottle of liquid and a small scale bottle of liquid 0.5L, I have 5 different samples from each and they are measured over the space of 8 days as % viability and the % viability decreases over time. However not all 10 samples got me

[R] Problem with "apply"

h mds[200:300,], mds[250:350,] and mds[1000:1100,], only ages with 3 digits are correctly assigned - all ages <100 are returned as NA. I'm using R v 2.8.1 on Windows XP. Cheers, Alan Cohen Centre for Global Health Research, Toronto,ON > ageassign <- function(x){ + y <- NA +

[R] Weighted principal components analysis?

v2.8.1 on Windows XP. I've searched Help and the R site and had no luck. Thanks for any help you can provide. Cheers, Alan Cohen Centre for Global Health Research Toronto, Ontario __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

[R] Fwd: 'for Loop'

Hello, A nice guy call Jun Shen was helping me out with this, but I require a bit more help. Below is my data set or list called 'test'. I'm trying to calculate the %RSD for each pair of index and keep it in cronological order if you can imagine a 3rd column with 'date' beside index. Res

[R] Using apply to get group means

in function? Any suggestions? Also, any way to avoid the hassles with having to convert to a data frame and then again to numeric when one variable is character? Cheers, Alan Cohen __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listin

[R] 'for Loop'

Hello, I'm trying to create a for loop for a data set, I have a list of results in this data set and I want to take the 1st two add them together and divide by the mean of the 1st to, then do the same for the 3rd and 4th values in the list and so on and each time return a value for the calcul

[R] Burt table from word frequency list

Maybe not terribly hard, depending on exactly what you need. Suppose you turn your text into a character vector 'mytext' of words. Then for a table of words appearing delta words apart (ordered), you can table mytext against itself with a lag: nwords=length(mytext) burttab=table(mytext[-(1:d

Re: [R] Is there any difference between <- and =

thanks for your supportive comments! by that time r programs will be scanned directly from your head, i suppose, and the intelligent scanner will as gladly take <- as it will =, so the problem will rather vanish. Yes, and maybe the scanner will be more intelligent than the programmer so when th

Re: [R] Is there any difference between <- and =

d keyboard that key code will have been assigned to the locally popular sheep icon.) Alan Zaslavsky __ 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/p

[R] basketball (was: R spam)

rhouses in the U.S. who probably know much more about this than I do. Alan Zaslavsky __ 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.ht

[R] understanding how R determines numbers and characters when creating a data frame

t;-data.frame(cbind(species=as.character(paste(s)),obsnum)) # number converted to character results<-rbind(out1,results) } results #fix(results) # can now convert obsnum to numeric using fix ## Thank you, Alan Smith [[alternative HTML version deleted]]

[R] AT&T Researchers and the New York Times

I've submitted a posting to the NY Times blog mentioned by David Smith, specifically mentioning John, Rick and Allan. (Now awaiting moderation.) Of course there are many others deserving credit but far to numerous to list so I didn't get started. On Thu, Jan 8, 2009 at 5:58 PM, Robert Wilkins

[R] R in the NY Times

This article is accompanied by nice pictures of Robert and Ross. Data Analysts Captivated by Power of R http://www.nytimes.com/2009/01/07/technology/business-computing/07program.html January 7, 2009 Data Analysts Captivated by R’s Power By ASHLEE VANCE To some people R is just the 18th letter

[R] Failure to subset in R v 2.8.0

cter(c$month)) does not help. Nor does restarting R, restarting the computer, or trying the operation on smaller subsets of c. I'd appreciate any help you an provide. Sincerely, Alan Cohen __ R-help@r-project.org mailing list https://stat.ethz.c

[R] Math Expression in 3D Plots

quot;black") title3d(xlab="dist", pos=c(NA, 0, -3), color="black") title3d(ylab="delta_yrsed", pos=c(12, NA, -3), color="black") title3d(zlab="delta_dist", pos=c(-3, 0, NA), color="black") rgl.bg(color=rep("white",

[R] ubuntu system refuses to update packages

't show what the issue is. -- --- | Alan K. Jackson| To see a World in a Grain of Sand | | [EMAIL PROTECTED] | And a Heaven in a Wild Flower, | | www.ajackson.org | Hold Infinity in the palm of your hand | | Houston, Texas

Re: [R] Updating Packages with Personal Library

Test. On Fri, Nov 7, 2008 at 2:40 AM, Alan Lue <[EMAIL PROTECTED]> wrote: > I'm having trouble running `updates.packages()' and installing into a > personal library. > > Setup: > 1. .Renviron file contains: R_LIBS_USER="~/lib/R/%p-library/%v" > 2.

[R] Updating Packages with Personal Library

pdated package was already in my personal library. How can I get R to check packages in my personal library while determining what to update? Alan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.c

[R] Memory limits for large data sets

R for this kind of thing? Cheers, Alan Cohen Alan Cohen Post-doctoral Fellow Centre for Global Health Research 70 Richmond St. East, Suite 202A Toronto, ON M5C 1N8 Canada (416) 854-3121 (cell) (416) 864-6060 ext. 3156 (0ffice) __ R-help@r-project.org mai

[R] Computing Mean Lifetime from Hazard

ted from a well-known distribution? Thanks for your help, Alan -- Alan Cox Director, User Experience iContact, Corp. p 919.459.1038 f 919.287.2475 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read t

[R] Confidence Intervals on Hazard Plots

te curves, haz1 looks very different than haz2, but I'd like to see if it's different enough.  How can I go about creating confidence intervals and plotting them? Thanks, Alan -- Alan Cox Director, User Experience iContact, Corp. p 919.459.1038 f 919.287.2475 _

[R] Hazard Rate Confidence Intervals

te curves, haz1 looks very different than haz2, but I'd like to see if it's different enough.  How can I go about creating confidence intervals and plotting them? Thanks, Alan -- Alan Cox Director, User Experience iContact, Corp. p 919.459.1038 f 919.287.2475 [[alt

[R] Projecting Survival Curve into the Future

x27;ll hit 50%? We started a new program 3.5 months ago, and I believe that this set of accounts behaves differently than the rest of our company's accounts. Thanks very much, Alan -- Alan Cox Director, User Experience iContact, Corp. p 919.459.1038 f 919.287.2475 [[a

[R] Constrained Optimization

Hello, I am trying to run a constrained optimization in R. "constrOptim" is really useful and has helped me a lot, but unfortunately, it doesn't provide the hessian. Is there a solution to this problem? I've tried "optim" with penalty-functions and "L-

Re: [R] Help with hazard plots

I clicked "Send" before making sure I thanked anyone who took the time to help me out.  Sorry about that.  To all who read or respond: thanks!. - Original Message - From: "Alan Cox" <[EMAIL PROTECTED]> To: r-help@r-project.org Sent: Thursday, July 31, 200

  1   2   >