Re: [R] Using R to Compute Covariance

2014-07-26 Thread peter dalgaard
-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. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School

Re: [R] Function assignment

2014-07-26 Thread peter dalgaard
. __ 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. -- Peter Dalgaard, Professor, Center for Statistics

Re: [R] I need help in seeing the code

2014-07-22 Thread peter dalgaard
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, Denmark Phone: (+45)38153501 Email

Re: [R] Maximum likelihood estimation (stats4::mle)

2014-07-22 Thread peter dalgaard
-guide.html and provide commented, minimal, self-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

[R] R 3.1.1 is released

2014-07-10 Thread Peter Dalgaard
for various platforms will appear in due course. Perhaps with some delays due to vacations. In particular, the Mac OS X maintainer is traveling and may be without Internet access for some days yet. For the R Core Team Peter Dalgaard These are the md5sums for the freshly created files, in case you

Re: [R] R Studio v3.0.3 for Windows 32bits is too slow

2014-07-09 Thread peter dalgaard
. __ 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. -- Peter Dalgaard, Professor, Center for Statistics

Re: [R] how does a valid subscript can produce an subscript out of bounds error?

2014-07-05 Thread peter dalgaard
/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, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

Re: [R] Fisher Scoring v/s Coordinate Descent for MLE in R

2014-07-04 Thread peter dalgaard
Vij __ 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. -- Peter Dalgaard

Re: [R] combining data from multiple read.delim() invocations.

2014-07-02 Thread peter dalgaard
than using 136 rbind() calls while expanding the data frame on each iteration. Presumably, most of the time is spent elsewhere, but 8s in the opposite direction is a little surprising. You might want to check the timing of just the df.all - do.call(rbind, df.list) bit. -- Peter Dalgaard

Re: [R] Sample all possible contingency tables both margin fixed

2014-06-24 Thread peter dalgaard
, minimal, self-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r

Re: [R] Reading in a csv - 2 different results

2014-06-22 Thread peter dalgaard
ways, or at the very least waste time trying to do so. As you probably don't need that, use lapply(). -pd -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda

Re: [R] R is converting arg input to scientific notation, which is bad!

2014-06-20 Thread peter dalgaard
-- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] extract a subset of non-contiguous elements of a matrix

2014-06-19 Thread peter dalgaard
, reproducible code. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org mailing list https

Re: [R] prediction based on conditional logistic regression, clogit

2014-06-18 Thread peter dalgaard
to cumulate the information represented by the log-likelihood contributions. I.e. you can do condtional logistic regression with a Cox likelihood because the Cox likelihood _is_ a conditional logistic regression likelihood. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business

Re: [R] prediction based on conditional logistic regression clogit

2014-06-16 Thread peter dalgaard
/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, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

Re: [R] correlation given p value and sample size

2014-06-16 Thread peter dalgaard
://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, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

Re: [R] question about chi values GLM

2014-06-14 Thread peter dalgaard
-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. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School

[R] Release of R 3.1.1 scheduled for July 10

2014-06-13 Thread Peter Dalgaard
Just a quick notice that we intend to have a patch release on July 10. Nickname is still undecided. Due to calendar conflicts, the binaries for some platforms may be delayed by a week or so. -pd -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3

Re: [R] Interval for rnorm command?

2014-06-09 Thread peter dalgaard
, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org mailing list https

Re: [R] R 3.1.0 interpreting large negative seconds since epoch as LMT?

2014-06-06 Thread peter dalgaard
of arc short of GMT+1 (at 15° 0´ E). If you're at about 74W (e.g., New York time), it sort of all figures. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

Re: [R] R 3.1.0 interpreting large negative seconds since epoch as LMT?

2014-06-06 Thread peter dalgaard
On 06 Jun 2014, at 16:46 , Prof Brian Ripley rip...@stats.ox.ac.uk wrote: so that says New York switched to EST in 1883 by 3:58 (not 4:02). As I also would have gotten, had I been able to subtract 44:27 from 48:25 properly Tnx -- Peter Dalgaard, Professor Center for Statistics

Re: [R] detect escape character

2014-06-05 Thread peter dalgaard
: Frede Aakmann Tøgersen Specialist, M.Sc., Ph.D. Plant Performance Modeling I find this very entertaining, because that's exactly what I do -- Plant Performance and Modeling -- but when I do it it involves botany. :) Sarah -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen

Re: [R] seeking latex peace

2014-05-25 Thread peter dalgaard
installing; it is so much easier to get things done right that way than it is to maintain your own texmf tree and/or mucking around in a system-installed one. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45

Re: [R] Inconsistent alphabetisation issue

2014-05-23 Thread peter dalgaard
. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] select randomly from a list

2014-05-22 Thread peter dalgaard
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. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3

Re: [R] How to compute a P-value for a complex mixture of chi-squared distributions in R

2014-05-22 Thread peter dalgaard
://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

Re: [R] Summary to data frame in R!!

2014-05-21 Thread peter dalgaard
On 07 May 2014, at 17:55 , Marc Schwartz marc_schwa...@me.com wrote: as.data.frame(sapply(mtcars[, c(mpg, disp)], summary)) I'd use lapply() rather than sapply() here. That goes directly from a list to a data frame, avoiding the detour via a matrix. -- Peter Dalgaard, Professor Center

Re: [R] approxfun question

2014-05-21 Thread peter dalgaard
-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org mailing

Re: [R] \ escape sequence and windows path

2014-05-20 Thread peter dalgaard
, minimal, self-contained, reproducible code. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r

Re: [R] Second axis on bottom of graph

2014-05-20 Thread peter dalgaard
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. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3

Re: [R] Second axis on bottom of graph

2014-05-20 Thread peter dalgaard
-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. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School

Re: [R] dev-lang/R-3.1.0: biocLite(vsn) removes all files in /

2014-05-18 Thread peter dalgaard
@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. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business

Re: [R] Merging data.table and data.frame

2014-05-13 Thread peter dalgaard
-contained, reproducible code. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org mailing list

Re: [R] File coding problem: how to read a windows-1252 encoded file

2014-05-13 Thread peter dalgaard
-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. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School

Re: [R] as.character(quote(x$y) ) = $ x y not x$y?

2014-05-10 Thread peter dalgaard
it, but the trouble is, is.is(is.is) is FALSE. Duncan Worse, is.is(mybaby) | is.aint(mybaby) comes out TRUE whatever the value of mybaby. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk

Re: [R] Installation problem

2014-05-07 Thread peter dalgaard
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 Statistics, Copenhagen Business School Solbjerg Plads 3, 2000

Re: [R] Lattice Histogram with Normal Curve - Y axis as percentages

2014-05-06 Thread peter dalgaard
and provide commented, minimal, self-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

Re: [R] Problem with products in R ?

2014-05-04 Thread peter dalgaard
, minimal, self-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r

Re: [R] Problem with products in R ?

2014-05-04 Thread peter dalgaard
]] __ 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. -- Peter Dalgaard, Professor, Center

Re: [R] correlation with missing values.. different answers

2014-04-14 Thread peter dalgaard
commented, minimal, self-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r

Re: [R] Selecting variables in a multivariate regression

2014-04-14 Thread peter dalgaard
. __ 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. -- Peter Dalgaard, Professor, Center for Statistics

[R] R 3.0.3 is released

2014-04-10 Thread Peter Dalgaard
for various platforms will appear in due course. (There seems to be a minor mishap with the NEWS file not getting copied to CRAN. I'll attend to that shortly.) For the R Core Team Peter Dalgaard These are the md5sums for the freshly created files, in case you wish to check

[R] R 3.1.0 is released

2014-04-10 Thread Peter Dalgaard
at a CRAN site nearer to you. Binaries for various platforms will appear in due course. (There seems to be a minor mishap with the NEWS file not getting copied to CRAN. I'll attend to that shortly.) For the R Core Team Peter Dalgaard These are the md5sums for the freshly created files, in case you

Re: [R] premature evaluation of symbols. Is that the way to describe this problem?

2014-04-10 Thread peter dalgaard
largely the point of the maneuver, but there is also a slight risk that the object could change or disappear due to side effects when evaluating other arguments. This isn't different from the usual perils of lazy evaluation, though. -pd -- Peter Dalgaard, Professor Center for Statistics

Re: [R] Determining 32- vs. 64-bit Windows

2014-04-10 Thread peter dalgaard
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. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3

Re: [R] Binom.test - hudge difference in p-value for little differences in PD forecast

2014-04-10 Thread peter dalgaard
it is in the second case. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org mailing list

Re: [R] Icelandic Characters in Mac

2014-04-03 Thread peter dalgaard
something like the described effect. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org mailing

Re: [R] A vector of normal distributed values with a sum-to-zero constraint

2014-04-02 Thread peter dalgaard
-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, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

Re: [R] cov.wt gives different results from other (co)variance functions (cov, wtd.var)

2014-03-31 Thread peter dalgaard
; in variance-weighting, I'd expect n/(n-1) (with n=117 here) which is not the same thing, unless the wt are equal. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda

Re: [R] GC overhead limit exceeded

2014-03-29 Thread peter dalgaard
it? For problems within packages, try the package maintainer. He might want to know what exactly you are trying to do. Most people here will not even realize that this is to do with the rrdf package. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000

Re: [R] ask for your help

2014-03-29 Thread peter dalgaard
. A Chinese R user lover. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

Re: [R] Problem with do.call().

2014-03-28 Thread peter dalgaard
(quote(plot), list(quote(x))) : 'what' must be a character string or a function which sees a bit overzealous. I see no reason that it shouldn't work.) -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45

Re: [R] How to import SAS file defining codes by SAS format catalog file

2014-03-22 Thread peter dalgaard
, USA __ 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. -- Peter Dalgaard

Re: [R] help

2014-03-19 Thread peter dalgaard
, then readHTMLTable() is not going to help you. So perhaps first see what is in there; either by opening the URL in a browser and using View Source, or by using con - url(url) ## you might want to reconsider the naming... readLines(con) close(con) -- Peter Dalgaard, Professor Center for Statistics

Re: [R] assigning dataframes in an ifelse statement

2014-03-17 Thread peter dalgaard
, and that they return a value.) -pd -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org

Re: [R] strangely long floating point with write.table()

2014-03-15 Thread peter dalgaard
, 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. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen

Re: [R] strangely long floating point with write.table()

2014-03-15 Thread peter dalgaard
On 15 Mar 2014, at 20:54 , Mike Miller mbmille...@gmail.com wrote: On Sat, 15 Mar 2014, peter dalgaard wrote: I don't think so. I think some of your numbers differ sufficiently from numbers with only a few digits to the right of the decimal that write.table needs to write them

Re: [R] How to get the adjusted R squared from GLS() estimator

2014-03-14 Thread peter dalgaard
unclear what the adj. R-squared should even mean. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r

[R] 3.1.0 on April 10

2014-03-10 Thread Peter Dalgaard
We intend to release R 3.1.0 on April 10. Prerelease versions will appear starting March 13 (for the source version, binaries may come a little later). See developer.r-project.org for details. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000

[R] R 3.0.3 is released

2014-03-06 Thread Peter Dalgaard
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) = cbf6da8f886ccd8d0dda0cc7ffd1b8ec MD5 (COPYING) = eb723b61539feef013de476e68b5c50a MD5

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

2014-03-06 Thread peter dalgaard
, self-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org

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

2014-03-06 Thread peter dalgaard
On 06 Mar 2014, at 19:16 , peter dalgaard pda...@gmail.com wrote: constant. Adding to that, and MA process with a neighbor correlation on the order of .9 is not possible. To paraphrase, you're doing essentially this: Sorry, dyslectic fingers syndrome: [...] _an_ MA process

Re: [R] Is this a mistake in 'An Introduction to R'?

2014-03-05 Thread peter dalgaard
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, Denmark Phone: (+45)38153501 Email: pd

Re: [R] data

2014-02-20 Thread peter dalgaard
-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. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School

Re: [R] Trouble with pmvnorm?

2014-02-10 Thread peter dalgaard
, at 18:04, peter dalgaard wrote: ou almost said it yourself: Your integrand doesn't vectorize. The direct culprit is the following: If x is a vector, what is lower=c(x,x,x,x)? A vector of length 4*length(x). And pmvnorm doesn't vectorize so it wouldn't help to have lower= as a matrix

Re: [R] Trouble with pmvnorm?

2014-02-08 Thread peter dalgaard
-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. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School

Re: [R] proportional weights

2014-02-06 Thread peter dalgaard
. __ 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. -- Peter Dalgaard, Professor

Re: [R] subset and na.rm not really suppressing NA values

2014-01-24 Thread peter dalgaard
://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, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

Re: [R] how to get the numbers of factors in a matrix

2014-01-22 Thread peter dalgaard
. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

Re: [R] how to get the numbers of factors in a matrix

2014-01-21 Thread peter dalgaard
and summarize those values. Actually, just factor(fff) will do. Or droplevels(fff). -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

Re: [R] Type III tests and Cox models

2014-01-20 Thread peter dalgaard
On 20 Jan 2014, at 18:47 , Terry Therneau thern...@mayo.edu wrote: The short summary: I was suspicious before, now I know for certain that it is misguided, and the phreg implementation of the idea is worse. A fortune candidate, if ever I saw one. -pd -- Peter Dalgaard, Professor, Center

Re: [R] Type III tests and Cox models

2014-01-20 Thread peter dalgaard
On 20 Jan 2014, at 23:05 , Göran Broström goran.brost...@umu.se wrote: On 01/20/2014 07:02 PM, peter dalgaard wrote: On 20 Jan 2014, at 18:47 , Terry Therneau thern...@mayo.edu wrote: The short summary: I was suspicious before, now I know for certain that it is misguided, and the phreg

Re: [R] Estimating parameters of 3 parameters lognormal distribution

2014-01-18 Thread peter dalgaard
commented, minimal, self-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help

Re: [R] (no subject)

2014-01-15 Thread peter dalgaard
__ 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. -- Peter Dalgaard, Professor Center for Statistics

Re: [R] SE for R-squared

2014-01-13 Thread peter dalgaard
@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. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School

Re: [R] Exact logistic regression, data frame issue

2014-01-11 Thread peter dalgaard
, but a little debugging suggests that line 3 gets zapped in the calculation of design.frame. Presumably it isn't happy with n=0. Passing dataset=test[-3,] makes something happen. -pd -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg

Re: [R] with() and within() functions inside lapply() not seeing outside of its environment?

2014-01-07 Thread peter dalgaard
. __ 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. -- Peter Dalgaard

Re: [R] p value of trends for odds ratios (or hazard ratios)

2014-01-04 Thread peter dalgaard
. __ 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. -- Peter Dalgaard, Professor, Center

Re: [R] p value of trends for odds ratios (or hazard ratios)

2014-01-04 Thread peter dalgaard
/4 peter dalgaard pda...@gmail.com On 04 Jan 2014, at 12:53 , Uwe Ligges lig...@statistik.tu-dortmund.de wrote: On 04.01.2014 06:39, zhu yao wrote: Dear Sir Many papers calculated the p value of trends for odds ratios of ordered category variables. I have found the tabodds command

Re: [R] Where did lost variables go, with example

2014-01-02 Thread peter dalgaard
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, Denmark Phone: (+45)38153501 Email: pd

Re: [R] The difference between SAS and R

2013-12-22 Thread peter dalgaard
, self-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org

Re: [R] charToRaw(Œ) is not 8C in R console

2013-12-13 Thread peter dalgaard
, self-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org

Re: [R] fisher.test - can I use non-integer expected values?

2013-12-11 Thread peter dalgaard
with simulated p-value (based on 1 replicates) data: c(1, 10, 9) X-squared = 7.3, df = NA, p-value = 0.0252 (The p-values _will_ differ because chi-square critical regions are slightly different from those based on the point probabilities.) -- Peter Dalgaard, Professor Center

Re: [R] Intepreting lm() results with factor

2013-12-03 Thread peter dalgaard
is zero (which is very rarely a relevant hypothesis). -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r

Re: [R] why change days of the week from a factor to an ordered factor?

2013-12-02 Thread peter dalgaard
(dodgers$ordered_day_of_week, levels=c(Monday, Tuesday, Wednesday, Thursday, Friday, Saturday, Sunday), labels=c(Mon, Tue, Wed, Thur, Fri, Sat, Sun)) And BTW, it doesn't (and didn't) create an ordered factor, just a factor with a different level ordering. -- Peter Dalgaard, Professor, Center

Re: [R] error in eval

2013-11-26 Thread peter dalgaard
, reproducible code. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org mailing list https

Re: [R] use of bquote

2013-11-22 Thread peter dalgaard
(pi==.(pi) ~ e==.(exp(1))) )) -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r-project.org

Re: [R] Recode values

2013-11-22 Thread peter dalgaard
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, Denmark Phone: (+45)38153501 Email

Re: [R] Does function read.sas7bdat() have some memory limitations?

2013-11-21 Thread peter dalgaard
__ 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. -- Peter Dalgaard, Professor Center

Re: [R] R- package for Parametric Survival Analysis with Left-censored data

2013-11-19 Thread peter dalgaard
/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, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

Re: [R] Bug in predict.lm?

2013-11-16 Thread peter dalgaard
for global variable 'w'). Apparently, the code checker is not smart enough to detect cases where a local variable is _sometimes_ not defined. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd

Re: [R] Double Pareto Log Normal Distribution

2013-11-12 Thread peter dalgaard
and provide commented, minimal, self-contained, reproducible code. -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com

Re: [R] Optimization failed in fitdistr (Weibull distribution)

2013-10-30 Thread peter dalgaard
On 29 Oct 2013, at 21:35 , Rolf Turner r.tur...@auckland.ac.nz wrote: On 10/29/13 19:44, peter dalgaard wrote: SNIP There really is no substitute for knowledge and understanding! Did it not occur to you that the Windspeed column needs to enter into your analysis? SNIP

Re: [R] Optimization failed in fitdistr (Weibull distribution)

2013-10-29 Thread peter dalgaard
, minimal, self-contained, reproducible code. -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv: pda...@gmail.com __ R-help@r

Re: [R] R vs octave development strategy (and success)

2013-10-29 Thread peter dalgaard
/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 Frederiksberg, Denmark

Re: [R] Nonparametric k-way ANOVA

2013-10-24 Thread peter dalgaard
/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, Denmark Phone: (+45

Re: [R] Function not working as I'd like

2013-10-23 Thread peter dalgaard
://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 Frederiksberg

Re: [R] How to view un-sampled data from a randomly sampled dataset

2013-10-23 Thread peter dalgaard
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, Denmark Phone: (+45)38153501 Email: pd@cbs.dk Priv

Re: [R] Where is element 30?

2013-10-23 Thread peter dalgaard
On Oct 23, 2013, at 20:04 , MacQueen, Don wrote: Follow these examples: grep(5,1:10) [1] 5 grep(3, c(1,5,2,3,6)) [1] 4 Don't: grep(5,1:15) [1] 5 15 -pd -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark

Re: [R] How to view un-sampled data from a randomly sampled dataset

2013-10-23 Thread peter dalgaard
), c(n, N-n))) y - data[ix,] x - data[!ix,] Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of peter dalgaard Sent: Wednesday, October 23, 2013 12:37 PM To: erinu Cc

<    4   5   6   7   8   9   10   11   12   13   >