Re: [R] Plotting Landscape in R-Studio

LEASE do read the posting guide >> http://www.R-project.org/posting-guide.html >> and provide commented, minimal, self-contained, reproducible code. > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > ht

Re: [R] Differential equations

ps://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. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Fre

[R] calculating the cumulative distribution function for a Von Mises Mixture model?

would be grateful for any suggestions, please? How is the direction for the lower limit for the CDF of a Von Mises Mixture model determined? Kind Regards Peter Mills University of Strathclyde Level 4, Technology and Innovation Centre (TIC), 99 George Street, Glasgow, G1 1RD

Re: [R] Marg.fct function

a) > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commente

Re: [R] Change (merge) vector by row names.

h(fullnames),1) rownames(extended) <- fullnames extended[rownames(reduced),] <- reduced where, presumably, you can find fullnames fromm the fit of a full model. -pd -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg

Re: [R] Error: long vectors (argument 1) are not supported in .Fortran

few weeks (possibly months), I do have a WGCNA update in the works that __should__ work on blocks larger than 46300. Best, Peter On Fri, Feb 3, 2017 at 7:31 AM, Bert Gunter <bgunter.4...@gmail.com> wrote: > Probably wrong list. Try the Bioconductor list instead. > > Cheers,

Re: [R] RQDA issue

That would be http://rqda.r-forge.r-project.org? It appears that the entire r-forge site is offline just now, but it is hosted at WU-Vienna, so they are the ones that have the hard facts about what has happened. Peter D. On 30 Jan 2017, at 08:09 , Miya Draga <dr...@ualberta.ca> wrote: &

Re: [R] Running Omega in R

>> >> >> [[alternative HTML version deleted]] >> >> __ >> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >> https://stat.ethz.ch/mailman/listinfo/r-help >> PLEASE do read the

Re: [R] Match ISO 8601 week-of-year numbers to month-of-year numbers on Windows with German locale

>>> LC_CTYPE=German_Germany.1252 >> LC_MONETARY=German_Germany.1252 >>>> [4] LC_NUMERIC=C LC_TIME=English_United >>>> States.1252 >>>> >>>> attached base packages:[1] stats graphics grDevices utils >

Re: [R] Tobit Regression with unbalanced Panel Data

e that the procedure has converged. Also, a logSigma of -5.4 suggests that you are working with small numbers -- it sometimes helps to scale things by a factor of 100 or 1000. -pd -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiks

Re: [R] if zlib version >= 1.2.5... no

This was fixed (independently?) in r-devel today: Peter-Dalgaards-MacBook-Air:R pd$ svn log m4/R.m4 | head r71889 | lawrence | 2017-01-04 04:57:31 +0100 (Wed, 04 Jan 2017) | 4 lines R_ZLIB macro tests ZLIB_VERNUM to handle

Re: [R] R in raspberry Pi

I can see the file under this link: http://www.floppybunny.org/robin/web/rbook/online_chapters/r_and_the_raspberry_pi.pdf Make sure the (English) words are not split - my first attempt contained raspber_ry and thus it failed. Peter On Mon, Jan 2, 2017 at 5:10 PM, John Sorkin <j

Re: [R] [FORGED] Export R output in Excel

rtment of Statistics > University of Auckland > Phone: +64-9-373-7599 ext. 88276 > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guid

Re: [R] histogram first bar wrong position

>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >> 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 c

Re: [R] data manipulation

//stat.ethz.ch/mailman/listinfo/r-help> >>> PLEASE do read the posting guide >>> http://www.R-project.org/posting-guide.html >>> <http://www.R-project.org/posting-guide.html> >>> and provide commented, minimal, self-contained, re

Re: [R] Infelicity in print output with matrix indexing of `[.data.frame`

re not coerced by being innocent bystanders of >> a dataframe colum that was not relevant. >> >> as.matrix( x[ min( i[ , 1]):max( i[ , 1]), min( i[ ,2 ]):max(i[ , 2]) >> ])[ >> cbind( i[,1]-min( i[ , 1]) +1, i[,2]- min( i[ ,2 ]) +1) ] > > ____

Re: [R] Log plus one transformation in R

ng guide http://www.R-project.org/ >> posting-guide.html >> and provide commented, minimal, self-contained, reproducible code. >> >> *Confidentiality Statement:* >> >> This email message, including any attachments, is for ...{{dropped:10}} > > _

Re: [R] Function implemented in R returns the wrong value

> > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-gu

Re: [R] Help!

___ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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. -- Peter Dalgaa

Re: [R] LaTeX errors when creating PDF version.

___ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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. -- Peter Dalgaar

Re: [R] create list of vectors in for loop

nts conversion to 10x10 matrix) or lapply(1:10, function(i) rnorm(10)) -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Email: pd@cbs.dk Priv: pda...@gmail.com __

Re: [R] i am trying to teach myself R

imal, self-contained, reproducible code. > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-g

Re: [R] errors when installing packages

> __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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-

Re: [R] GAM with the negative binomial distribution: why do predictions no match with original values?

at any > more, so I guess at some point I should remove this feature... > Fortune candidate... -pd > best, > Simon > > > > > On 23/11/16 00:24, Marine Regis wrote: >> Thanks a lot for your answers. >> Peter: sorry, here is the missing informati

Re: [R] Chi2 algorithm - R

___ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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

Re: [R] GAM with the negative binomial distribution: why do predictions no match with original values?

on deleted]] >> >> __ >> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >> https://stat.ethz.ch/mailman/listinfo/r-help >> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html >

Re: [R] using pmax in presence of NAs

__ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > 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. -- Peter Dalgaard, P

Re: [R] Gobbling up a repeating, irregular list of data

You will then have figure out how to deal with the empty cells in R. Peter On Thu, Nov 10, 2016 at 8:26 PM, Morway, Eric <emor...@usgs.gov> wrote: > What would be the sophisticated R method for reading the data shown below > into a list? The data is output from a numerical model. Pasti

Re: [R] Helmert contrasts and weighted means

erently, you would have to modify the entire design matrix if you removed a few rows. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Email: pd@cbs.dk Priv

Re: [R] creating lists of random matrices

Add a simplify = FALSE to the call to replicate, and you'll get a list. replicate(5, matrix(rnorm(4), 2, 2), simplify = FALSE) Peter On Wed, Nov 9, 2016 at 10:45 AM, Rui Barradas <ruipbarra...@sapo.pt> wrote: > Hello, > > I also thought of replicate() but it creates an 2x2x5 ar

Re: [R] Alternative to apply in base R

uide >> http://www.R-project.org/posting-guide.html >> and provide commented, minimal, self-contained, reproducible code. > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinf

Re: [R] problem installing R 2.5

lman/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. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, D

Re: [R] Sorting of character vectors

le("LC_COLLATE", "C"). That should work as long as you stay in good ol' ASCII. (2) Figure out (Don't look at me!) how to diddle the ICU settings for your system, icuSetCollate() is claimed to be your friend. -pd > Thanks for your help > > Pascal > >

Re: [R] multiplying a matrix by a vector

R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >>>> 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. &g

Re: [R] Updated package: survival_2.40-1

that ``we have found a bug in your code: > it gives a different answer than SAS''. (This is an actual quote.) Fortune nomination. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Em

[R] R 3.3.2 is released

rer to you. Binaries for various platforms will appear in due course. For the R Core Team, Peter Dalgaard These are the md5sums for the freshly created files, in case you wish to check that they are uncorrupted: MD5 (AUTHORS) = eb97a5cd38acb1cfc6408988bffef765 MD

Re: [R] cdf in R give probability of random variable

> On 21 Oct 2016, at 20:23 , peter dalgaard <pda...@gmail.com> wrote: > > In both cases it works out nicer if you do > > names(px) <- x > barplot(px) > Um, unless of course you want the cdf as a step function, in which case check the help page for plot f

Re: [R] cdf in R give probability of random variable

ed]] Notice that plain text posting is strongly preferred in here. HTML messes up code badly sometimes. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Email: pd@cbs.dk Priv: pda.

Re: [R] findInterval() surprising behavior

It's a bug (left.open=FALSE code gets executed in some cases). Hoping to have a fix tested and in place before 3.3.2. -pd > On 17 Oct 2016, at 21:48 , Rui Barradas wrote: > > Hello, > > Same on Windows 7. > > > findInterval(x=c(6, 1, 1, 1), vec=c(0, 1, 3, 5, 10),

Re: [R] nls.lm

f you just want the minimum or if you need to see the cov. matrix in order to see which parameters are unidentifiable. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Email: pd@cb

Re: [R] nls.lm

This would seem to apply to the add-on package minpack.lm. That package has a maintainer... Offhand, I would expect that this is a sanity check that, broadly speaking, prevents you from trying to solve a system of equations with more unknowns than equations. This is not a sufficient condition:

Re: [R] sort()ing strings

dardization. Names don't necessarily mean the same thing cross-platform. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Email: pd@cbs.dk

Re: [R] (no subject)

8 tools_3.2.5 Rcpp_0.12.7 > nlme_3.1-128grid_3.2.5 knitr_1.14 > [8] tensorA_0.36lattice_0.20-34 > > Can anyone explain what's going on? > > Bob > > -- > Bob O'Hara > > Biodiversity and Climate Research Centre > Senckenberganlage 25 >

Re: [R] Recoding lists of categories of a variable

** > This email and any attachments are confidential. Any use...{{dropped:8}} > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guid

Re: [R] Finding starting values for the parameters using nls() or nls2()

s://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. >> >> >> >> -- >> Andrew Robinson >> Deputy Di

Re: [R] how to work ttkspinbox ? win BUG ??

> On 09 Oct 2016, at 15:27 , Fox, John <j...@mcmaster.ca> wrote: > > Dear Peter, > > Thanks for the explanation. > > Is there a reason that there's no (unthemed) tkspinbox() in tclctk, defined > as tkspinbox <- function(parent, ...) tkwidget(parent, &

Re: [R] how to work ttkspinbox ? win BUG ??

displays a spinbox in an X-window: >>> >>> t <- tktoplevel() >>> s = ttkspinbox(t, from = 1.0, to = 100.0, textvariable='spinval') >>> tkpack(s) >>> >>> David Winsemius >>> Alameda, CA, USA >>> >>> &g

Re: [R] prcomp(): How do I multiply two matrices

to transpose the rotation matrix first. - Peter D. > On 07 Oct 2016, at 20:24 , T.Riedle <tr...@kent.ac.uk> wrote: > > Dear R-users, > > I am trying to do a principal components analysis using the attached data. My > code looks as follows. I want to calculate the time

[R] R-3.3.2 on October 31

Just a quick heads-up, mostly for those who want to keep their packages up to date with respect to updates of R: We intend to have a patch release on October 31. Nickname and detailed schedule will be made available on developer.r-project.org in due course. For the R Core Team Peter Dalgaard

[R] R-3.3.2 on October 31

Just a quick heads-up, mostly for those who want to keep their packages up to date with respect to updates of R: We intend to have a patch release on October 31. Nickname and detailed schedule will be made available on developer.r-project.org in due course. For the R Core Team Peter Dalgaard

Re: [R] What does lm() output coefficient mean when it's been given a categorical predictor of string values?

interpreted as 1 and 'poika' as 0? Short answer: Yes. Long answer: Yes, if treatment contrast parametrization is being used. See help(contrasts) for a lead-in to an even longer answer. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederik

Re: [R] Multivariable Wald to test equality of multinomial coefficients

[,1] [,2] [,3] [,4] [,5] [,6] [1,]100 -100 [2,]0100 -10 [3,]00100 -1 > A %*% beta [,1] [1,] 50.2709610 [2,] 22.5015565 [3,] -0.3689079 > t(A %*% beta) %*% solve(A %*% vcov(mod) %*% t(A), A %*% beta) [,1] [1,] 3.592326

Re: [R] Bootstrapping in R

s in coef(lm(df$y~df$x, data=d[f])). -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-

Re: [R] using read.csv2()

> On 29 Sep 2016, at 11:40 , Duncan Murdoch wrote: > > > It's a "factor". read.csv2() defaults to a decimal separator of "," rather > than ".", so the last column doesn't look like numbers, and they're being > read as character strings, and then automatically

Re: [R] Memory not release when an environment is created

>>> [[alternative HTML version deleted]] >>> __ >>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >>> https://stat.ethz.ch/mailman/listinfo/r-help >>> PLEASE do read the posting guide http:/

Re: [R] separate commands by semicolon

ts...) The comment char could also appear in a character string, where it does not mean the start of a comment... Not sure how to accomplish that using strsplit (or in general using just regular expressions). Peter __ R-help@r-project.org mailing list -- T

Re: [R] stfrtime function not returning proper results through sqldf package in R

On 16 Sep 2016, at 15:23 , PIKAL Petr <petr.pi...@precheza.cz> wrote: > Hi Peter > > The devil is in detail > > Data from OP had different format and was transferred to Date object by > as.Date, which results in incorrect values (and NA if not transferred) > df <

Re: [R] stfrtime function not returning proper results through sqldf package in R

ransfer of the email. > > In case that this e-mail forms part of business dealings: > - the sender reserves the right to end negotiations about entering into a > contract in any time, for any reason, and without stating any reasoning. > - if the e-mail contains an offer, the recip

Re: [R] Visualizing and clustering one half of a symmetric matrix

interpret the matrix as a distance matrix. Read carefully through the help on pheatmap to make sure the function plots what you want it to plot. HTH, Peter On Thu, Sep 15, 2016 at 7:38 PM, Khan, Saad M. (MU-Student) <smk...@mail.missouri.edu> wrote: > Hi all, > > I have a dist

[R] Creating dataframe with subtotals by all fields and totals of subtotals

an "All" salesmen subtotal and rbind it with the original dataset. I could do that for "Drink" and "Region", then also do combinations of salesmen, drink, and region subtotals. However, I'm hoping somebody out there is more clever than I am. Thanks! Peter dat <

Re: [R] Help with strftime error "character string is not in a standard unambiguous format"

anks to all who supplied > the answer and all who supplied additional useful information. Usually happens because you think it is p for print, so the other one must be f. (It's really _p_arse and _f_ormat.) -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solb

Re: [R] Help with strftime error "character string is not in a standard unambiguous format"

> On 12 Sep 2016, at 07:57 , Chris Evans <chrish...@psyctc.org> wrote: > >> strftime("Thu, 25 Aug 2016 6:34 PM",format="%a, %d %b %Y %I:%M %p") strptime, not strftime... -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business Schoo

Re: [R] understanding with

le expression. A little more than that, witness the difference between these two: > plot(airquality$Ozone) > with(airquality, plot(Ozone)) -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Offic

Re: [R] Using apply on a three dimensional matrix and passing multiple arguments to user defined function

would extend to using my mask function. Anyway, thanks again for both of your help. Cheers, Justin -- Justin Peter Research Fellow International Centre for Applied Climate Sciences, University of Southern Queensland West St, Toowoomba, QLD, 4350 Australia Email: justin.pe...@usq.

[R] Using apply on a three dimensional matrix and passing multiple arguments to user defined function

ssing? Thanks in advance for you help. Kind regards, Justin -- Justin Peter Research Fellow International Centre for Applied Climate Sciences, University of Southern Queensland West St, Toowoomba, QLD, 4350 Australia Email: justin.pe...@usq.edu.au<mailto:justin.pe...@usq.edu.au> Ph: +61

Re: [R] save R object into a remote directory

t; Stellenbosch University > South Africa > > Tel : +33 - (0)9 53 10 27 44 > Cell: +33 - (0)6 85 62 59 98 > Fax : +33 - (0)9 58 10 27 44 > > Fax (D):+49 - (0)3 21 21 25 22 44 > > email: rai...@krugs.de > > Skype: RMkrug > &

Re: [R] readBin documentation error

ossible character strings have a corresponding function. It is not clear why we allow both "int" and "integer" here, but there is no reason to expect int() to exist. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business Sc

Re: [R] Need advice on linear programming in R

t.org/web/views/Optimization.html#MathematicalProgrammingSolvers ? -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Email: pd@cbs.dk Priv: pda...@gmail.com __

Re: [R] GLM output problem

ept)x I(x^2) -0.44608 0.26767 -0.01458 -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Office: A 4.23 Email: pd@cbs.dk Priv: pda...@gmail.com __

Re: [R] Estimated Effects Not Balanced

<-rnorm(12)* > *mod<-aov(y~A+B)* As an aside, notice that this is one case where posting in HTML goes against your own interests (let alone the posting guide). The above cannot be cut+pasted, so fewer people will be willing to try it. -- Peter Dalgaard, Professor, Center for Statistics, Co

Re: [R] I had 1 problem with r

I think this is error. > Please give me answer, sir. > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE d

Re: [R] Efficiently parallelize across columns of a data.table

would not worry about copying the table to every slave node, since the copy doesn't really happen unless the worker processes modify the table. HTH, Peter On Fri, Aug 19, 2016 at 11:22 AM, Rebecca Payne <rebeccapa...@gmail.com> wrote: > I am trying to parallelize a task acros

[R] lavaan: Modification indices and intercepts

Hello, I’m examining CFA model invariance using the lavaan package in R. I’ve run an analysis to test “strong” model invariance (e.g., fixing loadings and intercepts to be equivalent across groups) but cannot seem to retrieve modification indices for freeing intercept parameters. Here is

Re: [R] need some help with date

t;MD", "NY", >"VA", "NY", "VT", "MD", "VT", "MD", "", "NY", "MD", "MD", >"NY", "VT", "VT", "NY", "OH", "CT", &

Re: [R] Help with non standard evaluation and require function

;- y + substitute (x + y) + } > f(z) x + 3.14159265358979 It is usually a better idea to handle these issues with bquote(), though. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg

Re: [R] Freq in dataframe

84643 > E-mail: m.trip...@aifa.gov.it > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-gu

Re: [R] Odds ratios in logistic regression models with interaction

d(coef(fit3.16), d) > se.ln.or <- sqrt(t(d) %*% vcov(fit3.16) %*% d) > exp(est.ln.or) > # [1,] 3.899427 > exp(est.ln.or - 1.96 * se.ln.or) > # [1,] 1.712885 > exp(est.ln.or + 1.96 * se.ln.or) > # [1,] 8.877148 > > -Original Message- > From: peter dalgaard [mai

Re: [R] Odds ratios in logistic regression models with interaction

(0.29, 3.79) > # 20 2.01 (0.91, 4.44) > # 25 3.90 (1.71, 8.88) > # 30 7.55 (1.95, 29.19) > > > > > [[alternative HTML version deleted]] > > ______ > R-help@r-project.org mailing list -- To UNSUBSCRIBE a

Re: [R] what happened to inside-r? [possibly OT]

mal, self-contained, reproducible code. >> > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.o

Re: [R] Aggregate matrix in a 2 by 2 manor

5.5 13.5 21.5 29.5 > aggregate.nx.ny.array.aperm(tst.small,FUN=colSums) [,1] [,2] [,3] [,4] [1,] 14 46 78 110 [2,] 22 54 86 118 > aggregate.nx.ny.forloop(tst.small,FUN=sum) [,1] [,2] [,3] [,4] [1,] 14 46 78 110 [2,] 22 54 86 118 cheers Peter > O

Re: [R] Aggregate matrix in a 2 by 2 manor

Hi all, thanks for the suggestions, I did some timing tests, see below. Unfortunately the aggregate.nx.ny.array.apply, does not produce the expected result. So the fastest seems to be the aggregate.nx.ny.expand.grid, though the double for loop is not that much slower. many thanks Peter >

Re: [R] Please advice me about nls.lm

I think you need someone who speaks your own language. Have a look at https://www.facebook.com/groups/KoreaRUsers/ - Peter D. > On 28 Jul 2016, at 08:54 , 박희국 <bbag...@gmail.com> wrote: > > Dear Sir/Madam > > > Dear > > Hello, I am a student in South Korea

Re: [R] Reduce woes

If you have a simple list of vectors (call it lst), use lengths = sapply(lst, length) In general, you may want to look at functions lapply and sapply which apply a function over a list, in this case the function length(). Peter On Wed, Jul 27, 2016 at 8:20 AM, Stefan Kruger <stefan.

[R] Aggregate matrix in a 2 by 2 manor

y over the new dimension, but that does not work, since it casts the data along the row. > matrix(apply(array(tst,dim=c(nx,ny,8)),3,mean),nrow=nrow(tst)/ny) [,1] [,2] [,3] [,4] [1,] 2.5 10.5 18.5 26.5 [2,] 6.5 14.5 22.5 30.5 cheers Peter __

Re: [R] removing all non-numeric characters from a string, but not "."

3.04" "12.5" -pd > > Thanks a lot! > > -- > Dimitri Liakhovitski > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read t

Re: [R] lm() silently drops NAs

ponse, Martin. > > Hadley > > -- > http://hadley.nz > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html >

[R] Regression with ARMA residual

n? Please advice Kind regards Peter 082 456 4669 -Original Message- From: Ismail SEZEN [mailto:sezenism...@gmail.com] Sent: 13 July 2016 04:55 PM To: Mangalani Peter Makananisa Cc: r-help@r-project.org Subject: Re: [R] Dates in R (Year Month) You can not convert numeric vectors directly

Re: [R] newbee Caret load question

read the posting guide http://www.R-project.org/posting-guide.html >> and provide commented, minimal, self-contained, reproducible code. > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > h

Re: [R] Dates in R (Year Month)

Thanks. -Original Message- From: David L Carlson [mailto:dcarl...@tamu.edu] Sent: 13 July 2016 05:05 PM To: Mangalani Peter Makananisa; r-help@r-project.org Subject: RE: Dates in R (Year Month) You need to look at the examples on the manual pages for ?yearmon and ?strptime: >

Re: [R] Reference for aov()

ect.org mailing list -- To UNSUBSCRIBE and more, see > 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. -- Peter Dalgaard, Professor, Center for Statis

[R] Dates in R (Year Month)

) library(chron, zoo) Z = as.yearmon(X) # it is not working please assist Kind regards Peter Please Note: This email and its contents are subject to our email legal notice which can be viewed at http://www.sars.gov.za/Pages/Email-disclaimer.aspx [[alternative HTML version deleted

[R] Fwd: How to make the "apply" faster

Forgot to cc the list... -- Forwarded message -- From: Peter Langfelder <peter.langfel...@gmail.com> Date: Sat, Jul 9, 2016 at 1:32 PM Subject: Re: [R] How to make the "apply" faster To: Debasish Pai Mazumder <pai1...@gmail.com> You could try the follo

Re: [R] 'ref' must be an existing level

Hmm, and levels(m11$prog)? (There's a chance that a space has sneaked in, but your HTML mail makes it hard to see column alignment) -pd > On 09 Jul 2016, at 12:50 , Shivi Bhatia <shivipm...@gmail.com> wrote: > > Hi Peter, > It gives me breakdown of all categories with the

Re: [R] 'ref' must be an existing level

lass of the var and is > a factor variable. > > [[alternative HTML version deleted]] > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-projec

Re: [R] tclArray() --- "active"?!

nd more, see > 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. -- Peter Dalgaard, Professor, Cente

Re: [R] BCa Bootstrapped regression coefficients from lmrob function not working

_ > De : Bert Gunter <bgunter.4...@gmail.com> > À : varin sacha <varinsa...@yahoo.fr> > Cc : peter dalgaard <pda...@gmail.com>; R-help Mailing List > <r-help@r-project.org> > Envoyé le : Mercredi 6 juillet 2016 1h19 > Objet : Re: [R] BCa Bootstrappe

Re: [R] BCa Bootstrapped regression coefficients from lmrob function not working

> boot.ci(results, type= "bca",index=2) > > > Any help would be highly appreciated, > S > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] strange behavior of lchoose in combinatorics problem

> On 25 Jun 2016, at 16:13 , Gonçalo Ferraz <gferra...@gmail.com> wrote: > > PROBLEM: the sum of the probabilities in lpvec should be <=1, but it is not. > The sum is something on the order of 1.48e-13. Um, in which sense is 1.48e-13 not <=1 ??? -- Peter Dal

Re: [R] [FORGED] Re: Generate list if sequence form two vector element

gt;> >>>> b<-c(10,7,20,2) >>>> >>>> >>>> Expected outcome : >>>> >>>> d<-list(1:10,3:7,6:20,2:9) > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE an

[R] R 3.3.1 is released

rer to you. Binaries for various platforms will appear in due course. For the R Core Team, Peter Dalgaard These are the md5sums for the freshly created files, in case you wish to check that they are uncorrupted: MD5 (AUTHORS) = eb97a5cd38acb1cfc6408988bffef765 MD5 (COPYING) = eb723b61539feef013de476e68

<    1   2   3   4   5   6   7   8   9   10   >