Re: [R] using subset() in data frame

r goes back to the .csv file which I read into R. For some reason it's propagating a space after macro that's not visible in the file, but R picks up. In future, I'll be careful to note how R expresses variables read from external files. Sincerely, Robert Walters _

Re: [R] using subset() in data frame

Chuck Cleland wrote: > On 2/23/2008 9:13 AM, Robert Walters wrote: >> Chuck Cleland wrote: Chuck, For the record, I might add that that the following two variants for subsetting worked equally well: fit1 <- lm(pore.pct ~ Db, subset(data.b, porosity == "macro ")) fi

[R] Kaiser-Meyer-Olkin

at case I suppose I should find out how to submit it. Robert "Tim" Kopp http://analytic.tripod.com/ [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

[R] can the matrix size limit be increased?

2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [, 25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [, 37] [,38] ... Thanks, Rob Robert W. Leach Scientific Programmer Center for

[R] Trying to "expand" some data - Newbie needs help

Woodley Sepulveda Van.Nuys Woodman Valley.College Laurel.Canyon North.Hollywood 1 9011880 5 2 2 2 2 2 2 22 2 6 4 1 2 9011890 1 1 1 1 1 1 2 1

Re: [R] Trying to "expand" some data - Newbie needs help

t;- apply(CurrentX1Sums, 1, sum) NewTargetX1Sums <- as.matrix(NewTargetData) NewTargetX2Sums <- apply(NewTargetX1Sums, 1, sum) Robert Farley Metro www.Metro.net -Original Message- From: Erik Iverson [mailto:[EMAIL PROTECTED] Sent: Wednesday, October 29, 2008 19:58 To: Farl

[R] Digest setting won't "take"

; box and change the radio button to "On". I click the "Submit My Changes" button which takes me back to the mailing list log-in page. I log back in and the setting is "Off". It looks like all settings are re-set to the default. How do I change these sett

Re: [R] Digest setting won't "take"

No Joy. :-( I've tried with both IE (7.0) and Firefox (3.0.3) and the same thing happens. Firefox has cookies enabled, and I set IE to explicitly allow http://www.r-project.org/ cookies. Robert Farley Metro www.Metro.net -Original Message- From: [EMAIL PROTECTED] [m

[R] Why are these different?

2 2 2 1 2 2 1 23 2 2 1 1 [4,] 9011970 1 1 1 1 2 1 1 26 2 2 2 24 [5,] 9012040 2 2 2 3 2

[R] R 2.8.0 for Ubuntu Intrepid Ibex i386

Hi all Is there going to be a binary of R 2.8.0 for Ubuntu Intrepid Ibex i386 - i.e., not just amd 64bit? Thx, RZ [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do r

[R] How do I read a text (.csv) file to match a matrix/cross tab? (Object confusion??)

lveda Van.Nuys Woodman Valley.College Laurel.Canyon North.Hollywood [1,] 9011880 5 2 2 2 2 2 2 22 2 6 4 1 [2,] 9011890 1 1 1 1 1 1 2

Re: [R] maptools sunrise sunset function

uncertainty. How SHOULD I specify the TZ? I would think that PST is a pretty "standard" ;-) time zone. Yes, I steal a lot of code snippets. I justify it by telling myself it helps me learn R. Robert Farley Metro www.Metro.net -Original Message- From: [EMAIL PROTECT

[R] Lattice: groups and no groups with panel.superpose

on how to "merge" the output from the two above graphs would be greatly appreciated! Kind regards, Robert Buitenwerf Ecologist SAEON South Africa _ [[elided Hotmail spam]] _

Re: [R] Lattice: groups and no groups with panel.superpose

Thanks Gabor and Baptiste, Both your suggestions give me the results I was looking for! Robert > CC: r-help@r-project.org > From: [EMAIL PROTECTED] > To: [EMAIL PROTECTED] > Subject: Re: [R] Lattice: groups and no groups with panel.superpose > Date: Sat, 15 Nov 2008 16:59:3

[R] Estimating the standard error when you have sampling weights.

Hi, Where can I find information ( freely available on the Internet , and also books or other sources ) on how having sampling weights changes the calculation of the standard error (of means and proportions)? How good is R for this type of procedure? And SAS? thanks Robert

[R] Cochran-Armitage

Hello, Which package allows to use Cochrana-Armitage trend test? I tried to search for but I found only package coin in which there is no explicit function. Best, RobMac __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] warnings although R CMD check runs without any problem

nce to a file, which triggers these messages, I > have currently no idea where to take a closer look at. Well, you did not give the whole output, and in it there really are names of files and locations. You can also use an editor to search for 1) any instance of the work exprSet or

Re: [R] Network Construction in R

can use biocViews to see what is available. best wishes Robert Johannes Urpelainen wrote: > Hi, > > I am trying to construct a social network from a data frame with rows > > acra numa acrb numb yearexpab eaboimpab iabo > 10 USA2 CAN 20 1957 4017

Re: [R] linewidth in Rgraphviz

you want... Robert Ben Bolker wrote: > > I'm an Rgraphviz newbie trying to figure out how to do some things/whether > they are possible. > Apparently graphviz does allow you to modify line widths of edges: see > > https://mailman.research.att.com/pipermail/graphviz

Re: [R] igraph and plotting connected components

> > [[alternative HTML version deleted]] > > __ > 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,

[R] calculate spatial distance

oordinates). thanks! Robert Robert Ptacnik Norwegian Institute for Water Research (NIVA) Gaustadalléen 21 NO-0349 Oslo mobile +47 982 277 81 FON NIVA +47 22 18 51 00 FAX +47 2

Re: [R] Splitting up the micEcon package?

so, nor it seems i´ll have it in the near future. On the other hand, judging by the low number of replies it seems to me that the audiece for SSel. in R is still quite limited. My two cents, greetings robert On 10/30/07, Arne Henningsen <[EMAIL PROTECTED]> wrote: > Dear R Users: > >

[R] Integrating a function

Hello everybody! I have problems with integrating my function. My primary function is a “survival function” of the following type: surviv <- exp(-k*x)/(1+exp(alpha*(x-tau))) I would like to integrate this function over a defined range and obtain a vector with all the values from integrate(survi

[R] bug in multcomp?

ny help would be much appreciated. Rob. -- Robert A. Cribbie, PhD Department of Psychology York University Toronto, ONM3J 1P3 Phone: 416-736-2100 x88615 Email: [EMAIL PROTECTED] Fax: 416-736-5814 Website: http://www.psych.yorku.ca/cribbie __ R-help@r-projec

[R] Integrating a two variable function

;- integrate(f1,x[i],x[i+1],int=0) result.1[i] <- integrate(f1,x[i],x[i+1],int=1) } Is there a way to take into account the second variable, i.e. for the integration to take into account int for each x ? Thanks in advance. Robert Schneider

[R] matafor package - categorical moderator interpretation question

uncertain if this is correct, and would appreciate any feedback. Bob ==== Robert Calin-Jageman Professor, Psychology Neuroscience Program Director Dominican University Parmer 210 7900 West Division River Forest, IL 60305 rcalinjage...@dom.edu 708.524.6581 http://calin-jageman.

[R] Environmental Data Connector v1.3

ult: Error in setwd(paste(Sys.getenv("EDC_HOME"), sep = "")) : cannot change working directory During installation EDC_HOME was set to /home/robert/EDC and the directory definitely exists. Does anyone have any suggestions? -- View this message in context: http://r.789695

Re: [R] Environmental Data Connector v1.3

Hi Dan, thanks for your response. The setwd is coded somewhere in the EDC.get function. I guess I could try alter the code but I assume this package should work as is. -- View this message in context: http://r.789695.n4.nabble.com/Environmental-Data-Connector-v1-3-tp4710686p4710701.html Sent

Re: [R] the less-than-minus gotcha

Rolf Turner is right on the money about not mixing-up '=' and '<-' Though this 'gotcha' will always a threat while '<-' is the assignment operator. The old Algol60 syntax of ':=' was less error-prone, but I guess '<-' is too firmly bedded-in to ever change. Meanwhile, spaces around the assi

[R] Solution to communicating with UDP and other interfaces (under Linux) using R

I'd spent hours with Google trying to find how to control a separate software package from R via a UDP interface (I could run the package under a "system" command, but that was too slow). I finally figured out a way to communicate with UDP through R, and it works fine (using the "system" comma

[R] R stops responding/communicating in for loop (lengthy description of issue)

Yeah, this is a bit lengthy, but it's a vexing problem. First, I'm working on learning R, mainly by using it and coming more from a programming aspect (I have the books and have gone through them, but learn best by doing). I have multiple projects going where R is almost necessary. I learned

[R] Solved: Communication in for() loop (Linux version)

Thanks for the advice I've received. FYI - It turns out that the problem was connected to the way R handles pipes and FIFOs, compared to the way the socat command does. (I don't know exactly what, but trying different things solved it!) I found that if you use FIFO() in R AND set up a FIFO c

Re: [R] Null and Alternate hypothesis for Significance test

//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. > [[alternative HTML version deleted]] _

[R] Specifying random effects distribution in glmer()

x27;t use nAGQ>1, not yet implemented summary(fit6) but this ends up with "Error in pgamma(id, shap, scal) : object "shap" not found". My questions are: 1. Can this be done? 2. Am I using the right package and function? 3. What am I doing wrong? Any help would be apprecia

Re: [R] Upgrading R means I lose my packages

point, the 2.7.2 has the same or new packages than 2.7.1, most or all of which will work. 5. I use the "Packages|Update Package ..." to update packages to 2.7.2. 6. Then I delete the 2.7.1 subfolder. You need Administrator rights to do this.

Re: [R] Projecting Survival Curve into the Future

Richie. Mathematical Sciences Unit HSL ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virg

Re: [R] bootstrapped correlation confint lower than -1 ?

E do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.

Re: [R] transformation/link glm in R

ling 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. ======== Robert A. LaBudde, PhD, PAS, Dpl.

Re: [R] online classes or online eduction in statistics? esp. time series analysis and cointegration?

nfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Co

[R] Problem accessing functions in package 'roxygen'

m" I have a similar problem with trim.right(), trim.left() and other functions I've tried. Any ideas? Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timber

Re: [R] Problem accessing functions in package 'roxygen'

' I ended up using trim <- function(x) gsub("^[[:space:]]+|[[:space:]]+$", "", x) instead. At 01:42 PM 9/3/2009, Uwe Ligges wrote: Robert A. LaBudde wrote: I have Vista Home with R-2.9.0, and installed and tried to test the package 'roxygen': > utils:::menu

[R] Confusion on use of permTS() in 'perm'

e-tail p-value of 4/120 = 0.0333, or a 2-tail p-value of 2*4/120 = 0.067. This is not, however, what permTS() returns. The permTS() value of 0.05 appears to correspond to 3 patterns, not 4. I am misunderstanding how to solve this simple problem, or is something going on with permTS() that I

Re: [R] multiple-test correlation

) Then do the pairwise tests without familywise error control. I.e., this is similar to doing the F test in ANOVA before doing LSD testing. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations,

Re: [R] Help with simple bootstrap test

mething to sum all the columns and rows of the table (or a subset should I desire). I've tried writing 'sum' for the function, but this gives me a huge figure which can't possibly be right. Robert A. LaBu

Re: [R] Nonparametric generalization of ANOVA

A search on "bluesky...@gmail.com" shows the user is in Norfolk, VA, USA. At 01:26 PM 3/5/2010, John Sorkin wrote: The sad part of this interchanges is that Blue Sky does not seem to be amiable to suggestion. He, or she, has not taken note, or responded to the fact that a number of people belie

[R] console window always on top

nyone know how to modify this setting? Thanks, -Robert Robert M. Flight, Ph.D. Bioinformatics and Biomedical Computing Laboratory University of Louisville Louisville, KY PH 502-852-0467 EM robert.fli...@louisville.edu EM rfligh...@gmail.com Williams and Holland's Law: If enou

Re: [R] console window always on top

I don't remember inputting anything like that previously, but that worked perfectly. Thank you. -Robert On Thu, Mar 11, 2010 at 13:12, Duncan Murdoch wrote: > On 11/03/2010 11:59 AM, Robert M. Flight wrote: >> >> In R on Windows (ver 7) I have somehow set the preferences

[R] list names

: matrices, R returns $M1 -> listing of numbers in matrix; $M2 -> listing of numbers in matrix, and $M3 -> listing of numbers in the matrix. Is there a way to have it return only the names $M1, $M2, $M3? I know this seems trivial and I'm sure it has been asked before, but I can

[R] Getting identify() to work with lattice::cloud()

ning: "no point within 0.25 inches". Does identify() work with cloud()? If so, what is the correct syntax. Thanks. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.U

Re: [R] How to test frequency independence (in a 2 by 2 table) with many missing values

osting-guide.html and provide commented, minimal, self-contained, reproducible code. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive

Re: [R] How to test frequency independence (in a 2 by 2 table) withmany missing values

tions that shows that the estimation of the fixed effects is unbiased, which assures you that you are estimating the "correct" coefficient. In addition, as Robert says, you should remove all subjects for which both observations are missing. ##First, we simulate some data. ##Note how the data s

[R] reading jpeg images?

software that it was itself unclear how to install. Is there some simpler solution I should be looking at? Are jpeg files so probematic I should be converting them to some other format and using a different package to read that? Thanks Robert Biddle

Re: [R] From THE R BOOK -> Warning: In eval(expr, envir, enclos) : non-integer #successes in a binomial glm!

read chapters 16 and 17 before deciding which example code to run. ======== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive

[R] respecting original matrix dimensions

,1] #tmp gives [,1] [1,]1 [2,]2 [3,]3 I don't see an easy way to do this beyond using matrix and supplying the correct matrix dimensions based on the indexing, or using a conditional to check if only one column/row is being taken out. Any help would be very much appreciated.

Re: [R] Random numbers with PDF of user-defined function

1.694695 ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-467-2947 "Vere scire est pe

Re: [R] Implementation of gamma function for large number

ain myself from asking the obvious question of what possible use gamma(17000) could be to you, and why a simple "infinity" would not work just as well. I'm sure you must have good reason for your request. ==== Rober

Re: [R] [Rd] NetCDF output in R

Hi Nana, This is not a r-devel question. I suspect this should be something like: write.netcdf.time(paste(path,'cam_fore.nc',sep=""), fore[,,times] ,lon,lat,times) Robert On Tue, Nov 10, 2009 at 8:31 AM, nana wrote: > Dear CSAG R users, > > I will be glad if so

Re: [R] random effects in mixed model not that 'random'

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. Robert

Re: [R] Block factor as random or fixed effect?

not, stick with the fixed effects model. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-46

Re: [R] "time series", "longitudinal data" or "trajectories"

e, as opposed to discrete times. "Trajectory" also implies an underlying causal model, as it is a term from kinematics. I hope this helps. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd

Re: [R] "time series", "longitudinal data" or "trajectories"

used. If there are multiple responses all measured at the same times for the subjects, the term "panel data" is used. For controlled experiments, the terms "repeated measures" and "time series" are common. "Longitudinal" could be used, but generally is not. =

Re: [R] Comparing model fits for NLME when models are not nested

lly significant"? ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239

Re: [R] 'stretching' a binomial variable

- any advice much appreciated! nicola Work with events instead of proportions, and use a Poisson model. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcflt

Re: [R] course in ecological data analysis

anyone can help me I will be most grateful. Best regards - J. Capelle Try http://www.statistics.com/ ======== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcflt

Re: [R] Genstat into R - Randomisation test

ing-based approximate method. The name does create confusion in terminology for "randomization" tests for bootstrapping. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Tim

Re: [R] p-values from bootstrap - what am I not understanding?

osting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.

Re: [R] Physical Units in Calculations

already exists. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-

[R] repeated measures anova

Hi I'm having difficulty working out how to get what I think is the appropriate partitioning of variability in a repeated measures setup. I have G=5 treatment-groups, each containing n=6 subjects, and a response is measured on each subject on t=4 occasions. I think the anova degrees of free

[R] repeated measures anova

Hi I'm having difficulty working out how to get what I think is the appropriate partitioning of variability in a repeated measures setup. I have G=5 treatment-groups, each containing n=6 subjects, and a response is measured on each subject on t=4 occasions. I think the anova degrees of freedom

Re: [R] Plotting survival curves after multiple imputation

Can anyone help with this ? On 14/02/2013 14:07, Robert Long wrote: I am working with some survival data with missing values. I am using the mice package to do multiple imputation. I have found code in this thread which handles pooling of the MI results: https://stat.ethz.ch/pipermail/r-help

Re: [R] significance testing for the difference in the ratio of means

__ >> 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. > &

[R] puzzling script bug

Windows XP (SP3) , R 2.15.1 32bit Hi ... I have a script which fails and closes my R session. Unfortunately, it bombs out at a different point each time I run it. I'm guessing that it may be something to do with memory management, or perhaps it's to do with the various .C dll's the script

Re: [R] puzzling script bug

Bob kinley -Original Message- From: Patrick Burns [mailto:pbu...@pburns.seanet.com] Sent: 21 September 2012 09:27 To: Robert Douglas Kinley Subject: Re: [R] puzzling script bug Bob, I haven't noticed if you've had any responses or not. The '.C' calls are almos

[R] Possible error in BCa method for confidence intervals in package 'boot'

argument 'eps' as a default 0.001/n in the calculations makes me suspect the programming in the script. The script apparently works correctly if the number of resamples equals or exceeds the number of original data, but not otherwise. Robert A. LaBudde, PhD, PAS,

[R] Mininum number of resamples required to do BCa bootstrap?

ed from the original sample via jackknife as 1/6 the coefficient of skewness.) The script apparently works correctly if the number of resamples equals or exceeds the number of original data, but not otherwise. Robert A. LaBudd

[R] R / JMP interface

hi useRs I am trying-out the facility to call R code from JMP. details: R 3.1.0 , JMP 11.1.1 , Windows 7 enterprise , all 64 bit. The test-script from the JMP help pages falls over at the first line :- R Init(); giving the error-message :- The installed version of R cannot be used. The entry

Re: [R] R / JMP interface

thanks Michael. That fits in with the response from JMP and with the experience of colleagues. Guess I'll just stick with R :) cheers Bob -Original Message- From: Meyners, Michael [mailto:meyner...@pg.com] Sent: 08 May 2014 07:20 To: Samuel J Gardner; Robert Do

[R] Errors Calculating MVN Likelihood of Time Series with AR(1) Errors

e means the proposed covariance matrix may not be invertible. This said I'm stuck on how to proceed and would be extremely appreciative of any thoughts. Thank you very much, Robert -- Robert Kindman Harvard College, Class of 2014 rkind...@college.harvard.edu +1.91

Re: [R] Errors Calculating MVN Likelihood of Time Series with AR(1) Errors

dhe.ibm.com/infocenter/spssstat/v20r0m0/index.jsp?topic=%2Fcom.ibm.spss.statistics.help%2Fcovariance_structures.htm) appears to be the better reference. I really appreciate your time, insight and help. Robert [[alternative HTML version deleted]]

Re: [R] logistic regression with 50 varaibales

ilman/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. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Lea

[R] Singularity in simple ANCOVA problem

; 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 0.1787 on 6 degrees of freedom Multiple R-squared: 0.7464, Adjusted R-squared: 0.535 F-statistic: 3.532 on 5 and 6 DF, p-value: 0.07811 Why the NA's here? ===

Re: [R] Singularity in simple ANCOVA problem

the groups. Rich Doh! Thanks. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia

Re: [R] Popularity of R, SAS, SPSS, Stata...

ship nearby you can see it happen. Once thing both whales and R users have in common is that, when you sight one, you say "R, Matey! Thar she blows!" ======== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com

Re: [R] Comparing distributions

__ 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] Function to compute the multinomial beta function?

)*Gamma(n3)/Gamma(n1+n2+n3) > beta3<- function (n1, n2, n3) exp(lgamma(n1)+lgamma(n2)+lgamma(n3)-lgamma(n1+n2+n3)) > beta3(5,3,8) [1] 1.850002e-07 ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com L

[R] adding environment variables

runtime. Thanks in advance, -Robert Robert M. Flight, Ph.D. University of Louisville Bioinformatics Laboratory University of Louisville Louisville, KY PH 502-852-0467 EM robert.fli...@louisville.edu EM rfligh...@gmail.com Williams and Holland's Law:        If enough data is collected,

Re: [R] How to do bootstrap for the complex sample design?

7;t be a problem. The fact that you choose a subsample size of 1 means you won't be able to estimate within-region variances unless you make some serious assumptions (e.g., UrbanCommunity effect independent of Region effect). ==== Robe

Re: [R] goodness-of-fit test

-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. ==== Robert A. LaBudde,

[R] Problem reading PostgreSQL data with RODBC

Dear list - I am having a problem using RODBC to access records from tables in a PostgreSQL database. There is no problem establishing the connection using chnl <- odbcConnect (dsn= ... etc. The DSN seems to be properly set up using the PostgreSQL Unicode ODBC driver, and sqlTables(chnl) works

Re: [R] When to use bootstrap confidence intervals?

E do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timb

Re: [R] Finding solution set of system of linear equations.

___ 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. =

[R] How to convert a factor column into a numeric one?

06547 22 -20H1 5.738683 37 -20H1 5.796394 4 14 -20H1 4.413691 50 45H1 6.406547 77 45H1 5.705433 but note that the values for 'Temp' in rows 5 and 7 are 45 and not 4, as expected, although the result is numeric. The

Re: [R] How to convert a factor column into a numeric one?

Exactly! Thanks. At 12:49 AM 6/5/2011, Jorge Ivan Velez wrote: Dr. LaBudde, Perhaps as.numeric(as.character(x)) is what you are looking for. HTH, Jorge On Sun, Jun 5, 2011 at 12:31 AM, Robert A. LaBudde <> wrote: I have a data frame: > head(df) Time Temp Conc ReplLog10 1

Re: [R] How to convert a factor column into a numeric one?

lnames(df)[5]<- 'Log10' At 12:51 AM 6/5/2011, Joshua Wiley wrote: Hi Robert, I would also look into *why* those numeric columns are being stored as factors in the first place. If you are reading the data in with read.table() or one of its wrapper functions (like read.csv), then it wou

Re: [R] How to convert a factor column into a numeric one?

- within(dd, { + a <- as.numeric(as.character(a)) + b <- as.numeric(as.character(b)) +} ) > str(de) 'data.frame': 20 obs. of 3 variables: $ a: num 1 1 1 1 2 2 2 2 3 3 ... $ b: num 1 1 2 2 1 1 2 2 1 1 ... $ y: num 0.6396 1.467 1.8403 -0.0915 0

Re: [R] P values

a, the Z and T distribution are identical. And it is only in Utopia that any P-value less than 0.01 actually corresponds to reality. ======== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.

Re: [R] P values

At 07:10 AM 5/7/2010, Duncan Murdoch wrote: Robert A LaBudde wrote: At 01:40 PM 5/6/2010, Joris Meys wrote: On Thu, May 6, 2010 at 6:09 PM, Greg Snow wrote: Because if you use the sample standard deviation then it is a t test not a z test. I'm doubting that seriously... You calc

Re: [R] Histogram Bin

thz.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. ==== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd

Re: [R] Cubic B-spline, how to numerically integrate?

quot; it. You need to give us the function and the appropriate limits of integration. Can anybody help please? Code and figures here: https://docs.google.com/fileview?id=0B1cQ7z9xYFl2ZTZhMmMyMjAtYTA3Zi00N2QyLTkxNzMtOGYyMjdiOGU2ZWE4&hl=en Thank you, Claudia ===

Re: [R] When normality test violate and sample size is large(n=400) can I use t-test

e simplest distribution-free tests should have near perfect power, so which test you use is not important. ======== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcflt

<    2   3   4   5   6   7   8   9   >