Re: [R] Very slow optim(): solved

2011-07-14 Thread Mike Marchywka
Date: Thu, 14 Jul 2011 12:44:18 -0800 From: toshihide.hamaz...@alaska.gov To: r-h...@stat.math.ethz.ch Subject: Re: [R] Very slow optim(): solved After Googling and trial and errors, the major cause of optimization was not functions, but data

Re: [R] Very slow optim()

2011-07-13 Thread Mike Marchywka
Hamazaki, Hamachan (DFG toshihide.hamazaki at alaska.gov writes: Dear list, I am using optim() function to MLE ~55 parameters, but it is very slow to converge (~ 25 min), whereas I can do the same in ~1 sec. using ADMB, and ~10 sec using MS EXCEL Solver. Are there any tricks

Re: [R] Using t tests

2011-07-10 Thread Mike Marchywka
( after getting confirmation of lack of posts try again, LOL ) From: marchy...@hotmail.com To: r-help@r-project.org Subject: RE: [R] Using t tests Date: Sun, 10 Jul 2011 10:13:51 -0400 ( sorry if this is a repost but I meant to post to list and

Re: [R] problem loading rgdal with Rapache, problem solved due to libexpat with apache.

2011-07-08 Thread Mike Marchywka
to know your intentions, platform, and the versions of the software you are using. Roger Mike Marchywka wrote: Has anyone had problems with Rapache that don't show up on command line execution of R? I just ran into this loading rgdal in Rapache page and having a problem with loading

Re: [R] superimposing network graphs

2011-07-06 Thread Mike Marchywka
Date: Wed, 6 Jul 2011 13:10:14 +0200 To: r-help@r-project.org CC: but...@uci.edu Subject: [R] superimposing network graphs Dear all, I have a undirected network (g), representing all the sexual relationships that ever existed in a model community. I also have a directed edgelist (e)

[R] problem loading rgdal with Rapache, problem solved due to libexpat with apache.

2011-07-06 Thread Mike Marchywka
Has anyone had problems with Rapache that don't show up on command line execution of R? I just ran into this loading rgdal in Rapache page and having a problem with loading shared object. The final complaint was that Stop_XMLParser was undefined- this was surprising since grep -l showed it in

Re: [R] wavelets

2011-07-04 Thread Mike Marchywka
From: jdnew...@dcn.davis.ca.us Date: Mon, 4 Jul 2011 00:45:41 -0700 To: tyagi...@gmail.com; r-help@r-project.org Subject: Re: [R] wavelets Study the topic more carefully, I suppose. My understanding is that wavelets do not in themselves compress anything, but because they sort out the

Re: [R] Protecting R code

2011-07-04 Thread Mike Marchywka
Put it on rapache or otherwise server but this seems like a waste depending on what you are doing Server side is only good way but making c++ may be interesting test Sent from my Verizon Wireless BlackBerry -Original Message- From: Vaishali Sadaphal vaishali.sadap...@tcs.com Date: Mon,

Re: [R] How to fit ARMA model

2011-07-01 Thread Mike Marchywka
From: bbol...@gmail.com Date: Fri, 1 Jul 2011 13:16:45 + Subject: Re: [R] How to fit ARMA model UnitRoot akhussanov at gmail.com writes: Hello, I am having some problems with fitting an ARMA model to my time series data (randomly generated numbers). The thing is I have

Re: [R] Running R from windows command prompt

2011-06-28 Thread Mike Marchywka
To: lig...@statistik.tu-dortmund.de CC: r-help@r-project.org Subject: Re: [R] Running R from windows command prompt Thanks for your help. I tried the way you mentioned for my first question. But I am not getting any results. Can you please explain in detail the process through which

Re: [R] Time-series analysis with treatment effects - statistical approach

2011-06-25 Thread Mike Marchywka
Date: Thu, 23 Jun 2011 15:41:25 -0700 From: jmo...@student.canterbury.ac.nz To: r-help@r-project.org Subject: Re: [R] Time-series analysis with treatment effects - statistical approach Mike Marchywka wrote: I discovered a way to do repetitive tasks that can be concisely

Re: [R] Access R functions from web

2011-06-25 Thread Mike Marchywka
From: orvaq...@gmail.com To: r-help@r-project.org Subject: [R] Access R functions from web I need a way to send R objects and call R functions from web. Is there any project close or similar to that? I want to be able to send an HTTP rquest from an existing application with some

Re: [R] Need help on a R script part

2011-06-25 Thread Mike Miller
this list the other day. Mike On Fri, 24 Jun 2011, Lusk Aris wrote: Hi all, I need all your help on this. I have the next part of code:e1=x1-mean(x1)e2=x2-mean(x2)n1=length(x1)n2=length(x2N=(n1 + n2)nu2=sum( c( ( x1 -mean(x1) )^2 , ( x2-mean(x2) )^2 ) )/Nss=c(e1,e2) b3= N*sum(ss^4

Re: [R] Is R a solution to my problem?

2011-06-24 Thread Mike Marchywka
From: jdnew...@dcn.davis.ca.us Date: Fri, 24 Jun 2011 02:13:29 -0700 To: ikoro...@gmail.com; R-help@r-project.org Subject: Re: [R] Is R a solution to my problem? Almost any tool can be bent to perform almost any task, but that combination sounds like a particularly difficult way to

Re: [R] R-help

2011-06-24 Thread Mike Marchywka
I'm surprised the pdf made it through mail list but many will not open these and I have just cleaned malware off my 'dohs system. Can you explain what you want to do in simple text? Subject: [R] R-help Hi Please assist me to code the attached pdf in R. Your help will be greatly

Re: [R] gaoshan

2011-06-24 Thread Mike Miller
On Fri, 24 Jun 2011, wang peter wrote: aa file is: x 1 NA 2 y 1 NA 3 and r program is aa-read.table(aa,row.names=1) bb-cor(t(aa),method = pearson,use=pairwise.complete.obs) bb x y x 1 1 y 1 1 i am confused why the pearson correlation coefficients between x and y is 1 You have two

Re: [R] Time-series analysis with treatment effects - statistical approach

2011-06-23 Thread Mike Marchywka
the crop then that is what you want to focus the analysis on etc etc etc. From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] on behalf of Mike Marchywka [marchy...@hotmail.com] Sent: Wednesday, June 22, 2011 9:31 PM To: jmo

Re: [R] Time-series analysis with treatment effects - statistical approach

2011-06-22 Thread Mike Marchywka
Subject: [R] Time-series analysis with treatment effects - statistical approach Hello all R listers, I'm struggling to select an appropriate statistical method for my data set. I have collected soil moisture measurements every hour for 2 years. There are 75 sensors taking these

Re: [R] several messages

2011-06-22 Thread Mike Miller
. That's still inefficient compared to the PLINK system, but it is way better than using doubles. Best, Mike -- Michael B. Miller, Ph.D. Minnesota Center for Twin and Family Research Department of Psychology University of Minnesota __ R-help@r-project.org

Re: [R] Re; Getting SNPS from PLINK to R

2011-06-22 Thread Mike Miller
integers, using 3 as the missing value. That's still inefficient compared to the PLINK system, but it is way better than using doubles. Best, Mike -- Michael B. Miller, Ph.D. Minnesota Center for Twin and Family Research Department of Psychology University of Minnesota

Re: [R] Time-series analysis with treatment effects - statistical approach

2011-06-22 Thread Mike Marchywka
Date: Wed, 22 Jun 2011 17:21:52 -0700 From: jmo...@student.canterbury.ac.nz To: r-help@r-project.org Subject: Re: [R] Time-series analysis with treatment effects - statistical approach Hi Mike, here's a sample of my data so that you get an idea what I'm working with. Thanks, data

Re: [R] Hardy Weinberg

2011-06-22 Thread Mike Miller
are in linkage equilibrium (an unfortunate term meaning that the genotypes for pairs of markers are not associated). Mike -- Michael B. Miller, Ph.D. Minnesota Center for Twin and Family Research Department of Psychology University of Minnesota __ R-help

Re: [R] Hardy Weinberg

2011-06-22 Thread Mike Miller
probably understand our questions and answer them. You might consider storing a 2x3 table as one row of length 6 in a matrix with 10,000 rows. Store that as one file -- it's only 60,000 elements, all of them integers. Store 1,000 files. You aren't being clear about what you are doing. Mike

Re: [R] different results from nls in 2.10.1 and 2.11.1

2011-06-20 Thread Mike Marchywka
Date: Mon, 20 Jun 2011 11:25:54 +0200 From: lig...@statistik.tu-dortmund.de To: p.e.b...@dunelm.org.uk CC: r-help@r-project.org; pb...@astro.uni-bonn.de Subject: Re: [R] different results from nls in 2.10.1 and 2.11.1 Since one is a 32-bit and the other one a 64-bit, and therefore the

Re: [R] Multivariate HPD credible volume -- is it computable in R?

2011-06-19 Thread Mike Marchywka
Date: Sun, 19 Jun 2011 19:06:20 +1200 From: agw1...@gmail.com To: r-help@r-project.org Subject: [R] Multivariate HPD credible volume -- is it computable in R? Hi all, I'm new to the list and am hoping to get some advice. I have a set of multivariate data and would like to find the

Re: [R] Server question

2011-06-18 Thread Mike Marchywka
From: oliver.jo...@digred.com To: r-help@r-project.org Date: Fri, 17 Jun 2011 15:35:38 +0100 CC: rebeccahard...@deltaeconomics.com Subject: [R] Server question Hi A client of mine has asked me to investigate the installation of R-software. Could anyone tell me

Re: [R] About 'hazard ratio', urgent~~~

2011-06-14 Thread Mike Marchywka
Date: Mon, 13 Jun 2011 19:44:15 -0700 From: dr.jz...@gmail.com To: r-help@r-project.org Subject: [R] About 'hazard ratio', urgent~~~ Hi, I am new to R. My question is: how to get the 'hazard ratio' using the 'coxph' function in

Re: [R] Need script to create new waypoint

2011-06-14 Thread Mike Marchywka
Now I need to change this data set into one with waypoints at regular intervals: for example 2 I guess your question is about algorithm ideas for making up data due to unspecified percrived need. Anything you can specify completely should be easy to code in R, maybe a bash script, or

Re: [R] About 'hazard ratio', urgent~~~

2011-06-14 Thread Mike Marchywka
Date: Tue, 14 Jun 2011 09:11:25 -0700 From: dr.jz...@gmail.com To: r-help@r-project.org Subject: Re: [R] About 'hazard ratio', urgent~~~ Thanks a lot for the great help~ well, if you are referring to my post, as I indicated I made a lot of

Re: [R] smoothScatter function (color density question) and adding a legend

2011-06-13 Thread Mike Marchywka
Date: Sat, 11 Jun 2011 21:52:21 -0400 From: ccoll...@purdue.edu To: r-help@r-project.org Subject: [R] smoothScatter function (color density question) and adding a legend Dear R experts, I am resending my questions below one more time

Re: [R] RES: Linear multivariate regression with Robust error

2011-06-11 Thread Mike Marchywka
Date: Fri, 10 Jun 2011 16:50:24 -0300 From: filipe.bote...@vpar.com.br To: frien...@yorku.ca; bkkoc...@gmail.com CC: r-help@r-project.org Subject: [R] RES: Linear multivariate regression with Robust error --Forwarded Message Attachment--

Re: [R] Amazon AWS, RGenoud, Parallel Computing

2011-06-11 Thread Mike Marchywka
Date: Sat, 11 Jun 2011 13:03:10 +0200 From: lui.r.proj...@googlemail.com To: r-help@r-project.org Subject: [R] Amazon AWS, RGenoud, Parallel Computing Dear R group, [...] I am a little bit puzzled now about what I could do... It seems like

Re: [R] Amazon AWS, RGenoud, Parallel Computing

2011-06-11 Thread Mike Marchywka
Date: Sat, 11 Jun 2011 19:57:47 +0200 Subject: Re: [R] Amazon AWS, RGenoud, Parallel Computing From: lui.r.proj...@googlemail.com To: marchy...@hotmail.com CC: r-help@r-project.org Hello Mike, [[elided Hotmail spam]] Best to my knowledge

Re: [R] error with geomap in googleVis

2011-06-10 Thread Mike Marchywka
Subject: RE: [R] error with geomap in googleVis Date: Fri, 10 Jun 2011 11:06:47 +0100 From: markus.gesm...@lloyds.com To: marchy...@hotmail.com; r-h...@stat.math.ethz.ch Hi Mike, I believe you are trying to put two charts on the same page

Re: [R] Linear multivariate regression with Robust error

2011-06-10 Thread Mike Marchywka
Date: Fri, 10 Jun 2011 09:53:20 +0530 From: bkkoc...@gmail.com To: r-help@r-project.org Subject: [R] Linear multivariate regression with Robust error Dear all, i am doing linear regression with robust error to know the effect of a (x)

Re: [R] Resources for utilizing multiple processors

2011-06-09 Thread Mike Marchywka
From: rjeffr...@ucla.edu Date: Wed, 8 Jun 2011 20:54:45 -0700 To: r-help@r-project.org Subject: [R] Resources for utilizing multiple processors Hello, I know of some various methods out there to utilize multiple processors but am not sure

Re: [R] error with geomap in googleVis

2011-06-09 Thread Mike Marchywka
To: r-h...@stat.math.ethz.ch From: mjphi...@tpg.com.au Date: Wed, 8 Jun 2011 10:14:01 + Subject: Re: [R] error with geomap in googleVis SNV Krishna primps.com.sg writes: Hi All, I am unable to get the plot geomap in googleVis package.

Re: [R] any documents

2011-06-09 Thread Mike Marchywka
Date: Thu, 9 Jun 2011 02:21:21 -0700 From: wa7@gmail.com To: r-help@r-project.org Subject: [R] any documents Hi, I'm doing a textual analysis of several articles discussing the evolution of prices in order to give a forecast. if

Re: [R] error with geomap in googleVis

2011-06-09 Thread Mike Marchywka
] On Behalf Of Mike Marchywka Sent: 09 June 2011 11:19 To: mjphi...@tpg.com.au; r-h...@stat.math.ethz.ch Subject: Re: [R] error with geomap in googleVis To: r-h...@stat.math.ethz.ch From: mjphi...@tpg.com.au Date: Wed, 8 Jun 2011 10:14:01 +

Re: [R] Can we prepare a questionaire in R

2011-06-08 Thread Mike Marchywka
Date: Wed, 8 Jun 2011 12:37:33 +0530 From: ammasamri...@gmail.com To: r-help@r-project.org Subject: [R] Can we prepare a questionaire in R Is there a way to prepare a questionnaire in R like html forms whose data can be directly populated into

Re: [R] RgoogleMaps Axes

2011-06-08 Thread Mike Marchywka
Date: Tue, 7 Jun 2011 09:50:10 -0700 From: egregory2...@yahoo.com To: r-help@r-project.org Subject: [R] RgoogleMaps Axes R Help, I posted a question on StackOverflow yesterday regarding an issue I've been having with the RgoogleMaps

Re: [R] Logistic Regression

2011-06-07 Thread Mike Marchywka
Date: Tue, 7 Jun 2011 01:38:32 -0700 From: farah.farid@student.aku.edu To: r-help@r-project.org Subject: [R] Logistic Regression I am working on my thesis in which i have couple of independent variables that are categorical in nature and the

Re: [R] Populating values from html

2011-06-07 Thread Mike Marchywka
Date: Tue, 7 Jun 2011 03:35:46 -0700 From: ammasamri...@gmail.com To: r-help@r-project.org Subject: [R] Populating values from html can we populate values into an excel sheet from html forms that has to be used in R for data analysis Can we

[R] predict with eha package

2011-06-02 Thread Mike Harwood
generated successfully, others output error f.ph.fit - predict(f.ph, newdata=mort, type='lp') f.ph.eha.fit - predict(f.ph.eha, newdata=mort, type='lp') f.aft.fit - predict(f.aft, newdata=mort, type='lp') Mike __ R-help@r-project.org mailing

Re: [R] error in model specification for cfa with lavaan-package

2011-06-01 Thread Mike Cheung
Dear Alain, There were 16 variables with 10 cases with missing values. The sample covariance matrix is not positive definite. It has nothing to do with lavaan. You need more cases before you can fit a CFA with 16 variables. Regards, Mike

Re: [R] Identifying sequences

2011-06-01 Thread Mike Marchywka
Date: Wed, 1 Jun 2011 17:12:29 +0200 From: cjp...@gmail.com To: r-help@r-project.org Subject: Re: [R] Identifying sequences Thanks to David, Thierry and Jonathan for your help. I have been able to put this function together a=1:10 b=20:30

Re: [R] error in model specification for cfa with lavaan-package

2011-06-01 Thread Mike Cheung
and/or the data are missing at random. By the way, if your questions are directly related to SEM, you may get more responses from SEMNET (http://alabamamaps.ua.edu/archives/semnet.html). Hope it helps. Regards, Mike -- - Mike W.L. Cheung

Re: [R] Compiling C-code in Windows

2011-05-31 Thread Mike Marchywka
Date: Tue, 31 May 2011 11:37:56 -0400 From: murdoch.dun...@gmail.com To: tom.osb...@iinet.net.au CC: r-h...@stat.math.ethz.ch; pmi...@ff.uns.ac.rs On 31/05/2011 7:50 AM, Tom Osborn wrote: You could use cygwin's cc/gcc, or the Watcom opensource

Re: [R] newbie: fourier series for time series data

2011-05-31 Thread Mike Marchywka
( hotmail won't mark text so I'm top posting... ) Can you post the data? Personally I'd just plot abs(fft(x)) and see what you see, as well as looking at Im(fft(x))/Re(fft(x)) or phase spectrum. Now, presumably you are nominally looking for something with a period of 1 year, that part could be

Re: [R] Value of 'pi'

2011-05-30 Thread Mike Marchywka
Date: Sun, 29 May 2011 23:09:47 -0700 From: jwiley.ps...@gmail.com To: vincy_p...@yahoo.ca CC: r-help@r-project.org Subject: Re: [R] Value of 'pi' Dear Vincy, I hope that in school you also learned that 22/7 is an approximation. Please

Re: [R] Nested design

2011-05-28 Thread Mike Marchywka
Date: Sat, 28 May 2011 09:33:03 -0700 From: jwiley.ps...@gmail.com To: bjorn.robr...@gmail.com CC: r-help@r-project.org Subject: Re: [R] Nested design Hi, If you are not asking for stats help, then do you understand the model and are

Re: [R] predictive accuracy

2011-05-26 Thread Mike Marchywka
Date: Thu, 26 May 2011 13:50:15 -0700 From: gunter.ber...@gene.com To: ahmed.el-taht...@pfizer.com CC: r-help@r-project.org Subject: Re: [R] predictive accuracy 1. This is not about R, and should be taken off list. Well, depending on what

Re: [R] Processing large datasets/ non answer but Q on writing data frame derivative.

2011-05-25 Thread Mike Marchywka
Date: Wed, 25 May 2011 09:49:00 -0400 From: ro...@bestroman.com To: biomathjda...@gmail.com CC: r-help@r-project.org Subject: Re: [R] Processing large datasets Thanks Jonathan. I'm already using RMySQL to load data for couple of days. I

Re: [R] Processing large datasets

2011-05-25 Thread Mike Marchywka
Date: Wed, 25 May 2011 10:18:48 -0400 From: ro...@bestroman.com To: mailinglist.honey...@gmail.com CC: r-help@r-project.org Subject: Re: [R] Processing large datasets Hi, If your datasets are *really* huge, check out some packages listed

Re: [R] Processing large datasets

2011-05-25 Thread Mike Marchywka
Date: Wed, 25 May 2011 12:32:37 -0400 Subject: Re: [R] Processing large datasets From: mailinglist.honey...@gmail.com To: marchy...@hotmail.com CC: ro...@bestroman.com; r-help@r-project.org Hi, On Wed, May 25, 2011 at 11:00 AM, Mike

[R] so current status of interactive PDF creation isl? trying to explore econometriic convergence thread

2011-05-23 Thread Mike Marchywka
I think I just saw a thread go by on how do I make interactive PDF from rgl output but I ignored it at the time and can't seem to find consensus result on google. Essentially I wanted to create pdf or other output to share a 3D plot and let viewers interact with it but it still wasn't clear on

Re: [R] RGL package installation problem on Centos

2011-05-23 Thread Mike Marchywka
Date: Mon, 23 May 2011 14:43:59 +0100 From: arraystrugg...@gmail.com To: r-help@r-project.org Subject: [R] RGL package installation problem on Centos Dear R users, I have installed the latest version of R from source on Centos (using configure

Re: [R] maximum likelihood convergence reproducing Anderson Blundell 1982 Econometrica R vs Stata

2011-05-22 Thread Mike Marchywka
with bobyqa in this example. In short, even with scaling and exact gradients, this optimization problem is recalcitrant. Best, Ravi. From: Mike Marchywka [marchy...@hotmail.com] Sent: Thursday, May 12, 2011 8:30 AM To: Ravi Varadhan; pda

[R] predict 'expected' with eha package

2011-05-21 Thread Mike Harwood
'. Since my 'real' data is left-truncated and right-censored I cannot use survreg, and I wanted to investigate the output from eha. Thanks in advance! Mike data(mort) aftreg(Surv(enter, exit, event) ~ ses, data = mort) Call: aftreg(formula = Surv(enter, exit, event) ~ ses, data = mort) Covariate

Re: [R] maximum likelihood convergence reproducing Anderson Blundell 1982 Econometrica R vs Stata

2011-05-21 Thread Mike Marchywka
. From: Mike Marchywka [marchy...@hotmail.com] Sent: Thursday, May 12, 2011 8:30 AM To: Ravi Varadhan; pda...@gmail.com; alex.ols...@gmail.com Cc: r-help@r-project.org Subject: RE: [R] maximum likelihood convergence reproducing Anderson Blundell 1982

Re: [R] text mining analysis and word visualization of pdfs

2011-05-19 Thread Mike Marchywka
Date: Wed, 18 May 2011 15:24:49 +0530 From: ashimkap...@gmail.com To: k...@huftis.org CC: r-h...@stat.math.ethz.ch Subject: Re: [R] text mining analysis and word visualization of pdfs On Wed, May 18, 2011 at 1:44 PM, Karl Ove Hufthammer wrote:

Re: [R] Help, please

2011-05-19 Thread Mike Marchywka
From: dwinsem...@comcast.net To: julio.flo...@spss.com.mx Date: Thu, 19 May 2011 10:40:08 -0400 CC: r-help@r-project.org Subject: Re: [R] Help, please On May 18, 2011, at 6:29 PM, Julio César Flores Castro wrote: Hi, I am using R

Re: [R] Power Spectrum from STFT (e1071)?

2011-05-13 Thread Mike Marchywka
Date: Tue, 10 May 2011 11:55:34 -0700 From: ajf...@psu.edu To: r-help@r-project.org Subject: [R] Power Spectrum from STFT (e1071)? Hello. Does anyone know how to generate a power spectrum from the STFT function in package e1071? The output for this function supplies the Fourier

Re: [R] Powerful PC to run R

2011-05-13 Thread Mike Marchywka
Date: Fri, 13 May 2011 12:38:51 +0200 From: haenl...@escpeurope.eu To: r-help@r-project.org Subject: [R] Powerful PC to run R Dear all, I'm currently running R on my laptop -- a Lenovo Thinkpad X201 (Intel Core i7 CPU, M620, 2.67 Ghz, 8 GB RAM). The problem is that some of my

Re: [R] How to extract information from the following dataset?

2011-05-12 Thread Mike Marchywka
Date: Thu, 12 May 2011 10:43:59 +0200 From: jose-marcio.mart...@mines-paristech.fr To: xzhan...@ucr.edu CC: r-help@r-project.org Subject: Re: [R] How to extract information from the following dataset? Xin Zhang wrote: Hi all, I have

Re: [R] maximum likelihood convergence reproducing Anderson Blundell 1982 Econometrica R vs Stata

2011-05-12 Thread Mike Marchywka
So what was the final verdict on this discussion? I kind of lost track if anyone has a minute to summarize and critique my summary below. Apparently there were two issues, the comparison between R and Stata was one issue and the optimum solution another. As I understand it, there was some

Re: [R] maximum likelihood convergence reproducing Anderson Blundell 1982 Econometrica R vs Stata

2011-05-09 Thread Mike Marchywka
Date: Mon, 9 May 2011 22:06:38 +1200 From: alex.ols...@gmail.com To: pda...@gmail.com CC: r-help@r-project.org; da...@otter-rsch.com Subject: Re: [R] maximum likelihood convergence reproducing Anderson Blundell 1982 Econometrica R vs Stata

[R] Syntax for iter.max in rms

2011-05-08 Thread Mike Harwood
you in advance, and all insights and criticisms are appreciated. Mike library(rms) f - bj(Surv(ftime, stroke) ~ rcs(age,5) + hospital,x=TRUE, y=TRUE) f - bj(Surv(ftime, stroke) ~ rcs(age,5) + hospital, link='identity',x=TRUE, y=TRUE) No convergence in 50 steps Failure in bj.fit f - bj(Surv

Re: [R] Confidence intervals and polynomial fits

2011-05-08 Thread Mike Marchywka
From: pda...@gmail.com Date: Sun, 8 May 2011 09:33:23 +0200 To: rh...@sticksoftware.com CC: r-help@r-project.org Subject: Re: [R] Confidence intervals and polynomial fits On May 7, 2011, at 16:15 , Ben Haller wrote: On May 6, 2011, at 4:27 PM, David Winsemius wrote: On May 6,

Re: [R] nls problem with R

2011-05-05 Thread Mike Marchywka
Date: Thu, 5 May 2011 01:20:33 -0700 From: sterles...@hotmail.com To: r-help@r-project.org Subject: Re: [R] nls problem with R ID1 ID2 t V(t) 1 1 0 6.053078443 2 1 0.3403 5.56937391 3 1 0.4181 5.45484486 4 1 0.4986 5.193124598 5 1 0.7451

Re: [R] ID parameter in model

2011-05-04 Thread Mike Harwood
(= 2.2.0), survival, graphics Maintainer: Göran Broström g...@stat.umu.se Packaged: 2011-03-01 14:56:12 UTC; gb Repository: CRAN Date/Publication: 2011-03-01 15:50:52 Built: R 2.13.0; i386-pc-mingw32; 2011-04-15 08:22:36 UTC; windows Mike On Mon, May 2, 2011 at 10:38 AM, Mike Harwood harwood

Re: [R] nls problem with R

2011-05-04 Thread Mike Marchywka
Date: Wed, 4 May 2011 07:07:44 -0700 From: sterles...@hotmail.com To: r-help@r-project.org Subject: Re: [R] nls problem with R Thanks Andrew. I am sorry for some typos that I omit some numbers of T2. Based on your suggestion,I think the problem is in the initial values. And I will read

Re: [R] UNIX-like cut command in R

2011-05-03 Thread Mike Miller
to correspond to cut -c, not to cut -f. Can it work with delimiters or only with character counts? Mike __ 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

Re: [R] UNIX-like cut command in R

2011-05-03 Thread Mike Miller
very close to what I'm looking for, but it seems to correspond to cut -c, not to cut -f. Can it work with delimiters or only with character counts? Mike x - this is a string unlist(strsplit(x, ))[c(1,4)] Thanks. I did figure that one out a couple of messages back, but to get it do behave

[R] ID parameter in model

2011-05-02 Thread Mike Harwood
, with the results at the bottom. Thanks in advance! Mike head(mort) ## data clearly contains multiple entries for some of the dataframe ids no.id.aft - aftreg(Surv(enter, exit, event) ~ ses, data = mort) ## Inital model id.aft - aftreg(Surv(enter, exit, event) ~ ses, data = mort, id=id

[R] UNIX-like cut command in R

2011-05-02 Thread Mike Miller
? I'm also guessing that there is a web page somewhere that will tell me how to do a lot of common GNU/UNIX/Linux text util commmand-line kinds of things in R. By that I mean by using R functions, not by making system calls. Does anyone know of such a web page? Thanks in advance. Mike

[R] List of Data Frames

2011-05-02 Thread Mike Smith
on each week, but have an average running up to that point in the season.  Thus the list would be indexed by weeks, and then there's a data frame of the game and all relevant statistics.  Thank You, Mike Smith __ R-help@r-project.org mailing list https

Re: [R] UNIX-like cut command in R

2011-05-02 Thread Mike Miller
with a vector of strings. Mike __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] UNIX-like cut command in R

2011-05-02 Thread Mike Miller
On Mon, 2 May 2011, Gabor Grothendieck wrote: On Mon, May 2, 2011 at 10:32 PM, Mike Miller mbmille...@gmail.com wrote: On Tue, 3 May 2011, Andrew Robinson wrote: try substr() OK.  Apparently, it allows things like this... substr(abcdef,2,4) [1] bcd ...which is like this: echo abcdef

Re: [R] how to generate a normal distribution with mean=1, min=0.2, max=0.8

2011-04-29 Thread Mike Miller
the correct expected max and min? That's another legitimate question. Mike -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Mao Jianfeng Sent: Thursday, April 28, 2011 12:02 PM To: r-help@r-project.org Subject: [R] how

Re: [R] how to generate a normal distribution with mean=1, min=0.2, max=0.8

2011-04-29 Thread Mike Miller
On Fri, 29 Apr 2011, David Winsemius wrote: On Apr 29, 2011, at 1:29 PM, Mike Miller wrote: On Fri, 29 Apr 2011, Giovanni Petris wrote: Well, but the original poster also refers to 0.2 and 0.8 as expected min and max, in which case we are back to a joke... Well, he is a lot better

Re: [R] fisher exact for 2x2 table

2011-04-29 Thread Mike Miller
is exceptionally clear on this point, which may be true, but what is the point? The point is that he prefers a different *name* for the test than the rest of us. Everyone agrees on the math/stat. Mike -- Michael B. Miller, Ph.D. Minnesota Center for Twin and Family Research Department of Psychology

Re: [R] regular expression in gsub() for strings with leading backslash

2011-04-29 Thread Mike Miller
set. Mike -- Michael B. Miller, Ph.D. Minnesota Center for Twin and Family Research Department of Psychology University of Minnesota __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] fisher exact for 2x2 table

2011-04-28 Thread Mike Miller
tables and larger sample sizes the computational cost becomes quite enormous. Mike __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

Re: [R] fisher exact for 2x2 table

2011-04-28 Thread Mike Miller
On Fri, 29 Apr 2011, Thomas Lumley wrote: On Fri, Apr 29, 2011 at 8:01 AM, Mike Miller mbmille...@gmail.com wrote: On Thu, 28 Apr 2011, viostorm wrote: I'm using fisher.exact on a 4x2 table and it seems to work. Does anyone know exactly what is going on?  I thought fisher.exact is only

Re: [R] fisher exact for 2x2 table

2011-04-28 Thread Mike Miller
been adamant when s/he was wrong. Show your biostatistician the MathWorld page. Mike __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

Re: [R] how to generate a normal distribution with mean=1, min=0.2, max=0.8

2011-04-28 Thread Mike Miller
. Mike -- Michael B. Miller, Ph.D. Bioinformatics Specialist Minnesota Center for Twin and Family Research Department of Psychology University of Minnesota On Thu, 28 Apr 2011, Carl Witthoft wrote: That method (creating lots of samples and throwing most of them away) is usually frowned upon

Re: [R] how to generate a normal distribution with mean=1, min=0.2, max=0.8

2011-04-28 Thread Mike Miller
- qnorm(runif(n, p_L, p_U), mean=m, sd=s) Or it could be written on one line: x - qnorm(runif(n, pnorm(L, mean=m, sd=s), pnorm(U, mean=m, sd=s)), mean=m, sd=s) Mike On Thu, 28 Apr 2011, Mike Miller wrote: Good point. It would be absurdly inefficient if the upper and lower limits

Re: [R] Speed up plotting to MSWindows graphics window

2011-04-27 Thread Mike Marchywka
Date: Wed, 27 Apr 2011 11:16:26 +0200 From: jonat...@k-m-p.nl To: r-help@r-project.org Subject: [R] Speed up plotting to MSWindows graphics window Hello, I am working on a project analysing the performance of motor-vehicles through messages logged over a CAN bus. I am using R 2.12 on

Re: [R] Speed up plotting to MSWindows graphics window

2011-04-27 Thread Mike Marchywka
Date: Wed, 27 Apr 2011 14:40:23 +0200 From: jonat...@k-m-p.nl To: r-help@r-project.org Subject: Re: [R] Speed up plotting to MSWindows graphics window On 27/04/2011 13:18, Mike Marchywka wrote: Date: Wed, 27 Apr 2011 11:16:26 +0200

Re: [R] blank space escape sequence in R?

2011-04-25 Thread Mike Miller
to evoke the hexadecimal ascii form of blank? That works in perl as does \040 for the octal form. Mike __ 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

Re: [R] Survival analysis: same subject with multiple treatments and experience multiple events

2011-04-23 Thread Mike Marchywka
Date: Fri, 22 Apr 2011 10:00:31 -0400 From: littleduc...@gmail.com To: r-help@r-project.org Subject: [R] Survival analysis: same subject with multiple treatments and experience multiple events Hi there, I need some help to figure out what

Re: [R] Using Java methods in R

2011-04-23 Thread Mike Marchywka
Date: Sat, 23 Apr 2011 05:32:59 -0700 From: hill0...@umn.edu To: r-help@r-project.org Subject: Re: [R] Using Java methods in R No answer to my post, so let's try a simpler question. Am I doing this correctly? I have the RGui with R

Re: [R] Using Java methods in R

2011-04-23 Thread Mike Marchywka
From: marchy...@hotmail.com To: hill0...@umn.edu; r-help@r-project.org Date: Sat, 23 Apr 2011 15:12:30 -0400 Subject: Re: [R] Using Java methods in R Date: Sat, 23 Apr 2011

Re: [R] How to answer the question about transitive correlation?

2011-04-22 Thread Mike Marchywka
Date: Fri, 22 Apr 2011 11:42:59 +0800 From: mailzhu...@gmail.com To: r-help@r-project.org Subject: [R] How to answer the question about transitive correlation? Hi, everyone. I know it may be a basic statistical question. But I can't find a

Re: [R] survexp with weights

2011-04-21 Thread Mike Harwood
system i386, mingw32 status major 2 minor 12.2 year 2011 month 02 day25 svn rev54585 language R version.string R version 2.12.2 (2011-02-25) package ‘survival’ version 2.36-5 On Apr 20, 11:03 am, Mike Harwood harwood

Re: [R] Extrapolating data points for individuals who lack them

2011-04-20 Thread Mike Marchywka
From: de...@exeter.ac.uk To: r-help@r-project.org Date: Wed, 20 Apr 2011 12:41:29 +0100 Subject: [R] Extrapolating data points for individuals who lack them Hi, We have an experiment where individuals responses were measured over 5 days. Some

[R] survexp with weights

2011-04-20 Thread Mike Harwood
Hello, I probably have a syntax error in trying to generate an expected survival curve from a weighted cox model, but I can't see it. I used the help sample code to generate a weighted model, with the addition of a weights=albumin argument (I only chose albumin because it had no missing values,

Re: [R] regression and lmer

2011-04-18 Thread Mike Marchywka
Date: Mon, 18 Apr 2011 03:27:40 -0700 From: lampria...@yahoo.com To: r-help@r-project.org Subject: Re: [R] regression and lmer Dear all, I hope this is the right place to ask this question.   ( hotmail not marking your text sorry I can';t find

Re: [R] Rsquared for anova

2011-04-17 Thread Mike Marchywka
( did this msg make it through the lists as rich text? hotmail didn't seem to think it was plain text?) Anyway, having come in in the middle of this it isn't clear if your issues are with R or stats or both. Usually the hard core stats people punt the stats questions to other places but both

<    1   2   3   4   5   6   7   8   9   >