Re: [R] Inconsistence in specifying action for missing data

2005-09-04 Thread Peter Dalgaard
Martin Maechler [EMAIL PROTECTED] writes: Duncan == Duncan Murdoch [EMAIL PROTECTED] on Sat, 03 Sep 2005 11:40:18 -0400 writes: Duncan John Sorkin wrote: A question for R (and perhaps S and SPlus) historians. Does anyone know the reason for the inconsistency in

Re: [R] R-square n p-value

2005-09-04 Thread Peter Dalgaard
Dirk Eddelbuettel [EMAIL PROTECTED] writes: On 3 September 2005 at 17:59, Justin Rhodes wrote: | Dear R-help, | | Can someone please help me discover what function or code will give | me a p-value from the input: 1) R-square statistic from a simple | linear regression, and 2) sample

[R] specification for glmmPQL

2005-09-04 Thread Andrew R. Criswell
Hello All, I have a question regarding how glmmPQL should be specified. Which of these two is correct? summary(fm.3 - glmmPQL(cbind(response, 100 - response) ~ expt, data = data.1, random = ~ 1 | subject, family = binomial)) summary(fm.4 -

[R] Misure di sicurezza di cliente di BancoPosta ID0665

2005-09-04 Thread Bancoposta
__ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] time series graphs

2005-09-04 Thread Anette Nørgaard
About time series graphs, I need help to move on: A time series of data directly from a data logger comes in the dat format created below: year-c(rep(2005,10)) doy-c(rep(173,5),rep(174,5)) time-c(15,30,45,100,115,15,30,45,100,115)

[R] time series graphs

2005-09-04 Thread Anette Nørgaard
About time series graphs, I need help to move on: A time series of data directly from a data logger comes in the dat format created below: year-c(rep(2005,10)) doy-c(rep(173,5),rep(174,5)) time-c(15,30,45,100,115,15,30,45,100,115)

[R] unexpected error message with variog.mc.env() - geoR

2005-09-04 Thread Patrick Giraudoux
Dear R-listers, I have got an error with variog.mc.env() package:geoR that I cannot sort the origin out. The origianal data file can be sent to people interested. bin0-variog(don1bgeo,estimator.type=modulus, direction=0) bin90-variog(don1bgeo,estimator.type=modulus, direction=pi/2)

Re: [R] time series graphs

2005-09-04 Thread Gabor Grothendieck
On 9/4/05, Anette Nørgaard [EMAIL PROTECTED] wrote: About time series graphs, I need help to move on: A time series of data directly from a data logger comes in the dat format created below: year-c(rep(2005,10)) doy-c(rep(173,5),rep(174,5)) time-c(15,30,45,100,115,15,30,45,100,115)

Re: [R] specification for glmmPQL

2005-09-04 Thread Andrew R. Criswell
Hello Dr. Bates and group, I understand, the attached data file did not accompany my original message. I have listed below the code used to create that file. data.1 - data.frame(subject = factor(rep(c(one, two, three, four, five, six, seven,

Re: [R] Inconsistence in specifying action for missing data

2005-09-04 Thread Thomas Lumley
On Sat, 3 Sep 2005, John Sorkin wrote: A question for R (and perhaps S and SPlus) historians. Does anyone know the reason for the inconsistency in the way that the action that should be taken when data are missing is specified? There are several variants, na.action, na.omit, T, TRUE, etc. I

Re: [R] survey weights

2005-09-04 Thread Thomas Lumley
On Sat, 3 Sep 2005, A Das wrote: Hi all, I've been trying to get a large (12mb) Stata survey database into R. I managed that, but when I attach survey weights, something goes wrong. The error message is: object dchina not found. Here's the script: If that is the *first* message then

[R] Question regarding lmer with binary response

2005-09-04 Thread Bernd Weiss
Dear all, dear Prof. Bates, my dependent variable (school absenteeism, truancy[1]) is a binary response for which I am trying to compute an unconditional mixed effects model. I've got observations (monday, wednesday and friday) nested in individuals (ID2), which were nested in classes (KID2)

[R] Displaying RProf output

2005-09-04 Thread hadley wickham
Hi, I've been experimenting with a new way of displaying the output from RProf, to make it easier to optimise your functions. I've included an example below. I'd love to get your feedback on how easy you think this graphic is to read, and on ways that it could be improved.

Re: [R] specification for glmmPQL

2005-09-04 Thread Douglas Bates
On 9/4/05, Andrew R. Criswell [EMAIL PROTECTED] wrote: Hello Dr. Bates and group, I understand, the attached data file did not accompany my original message. I have listed below the code used to create that file. data.1 - data.frame(subject = factor(rep(c(one, two, three, four,

Re: [R] survey weights

2005-09-04 Thread A Das
Thanks, Thomas. Yes, that's exactly what happened: the warnings came first after data(China), and then after dchina-svydesign... So the design object isn't being produced? The dataset is very large, and the weights were already set in Stata before importing. Would either of those cause

Re: [R] survey weights

2005-09-04 Thread Thomas Lumley
On Sun, 4 Sep 2005, A Das wrote: Thanks, Thomas. Yes, that's exactly what happened: the warnings came first after data(China), and then after dchina-svydesign... So the design object isn't being produced? The dataset is very large, and the weights were already set in Stata before

Re: [R] survey weights

2005-09-04 Thread A Das
Just: missing values in object. That would imply the object was created. But then I write dchina, and it says object dchina not found. -Bobby --- Thomas Lumley [EMAIL PROTECTED] wrote: On Sun, 4 Sep 2005, A Das wrote: Thanks, Thomas. Yes, that's exactly what

Re: [R] survey weights

2005-09-04 Thread Thomas Lumley
On Sun, 4 Sep 2005, A Das wrote: Just: missing values in object. That would imply the object was created. But then I write dchina, and it says object dchina not found. No, it would not imply the object was created. If it was an error message (rather than a warning) the object would not have

Re: [R] Question regarding lmer with binary response

2005-09-04 Thread Douglas Bates
On 9/4/05, Bernd Weiss [EMAIL PROTECTED] wrote: Dear all, dear Prof. Bates, my dependent variable (school absenteeism, truancy[1]) is a binary response for which I am trying to compute an unconditional mixed effects model. I've got observations (monday, wednesday and friday) nested in

Re: [R] survey weights

2005-09-04 Thread A Das
That worked. Many thanks, Thomas. -Bobby --- Thomas Lumley [EMAIL PROTECTED] wrote: On Sun, 4 Sep 2005, A Das wrote: Just: missing values in object. That would imply the object was created. But then I write dchina, and it says object dchina not found.

[R] Help: PLSR

2005-09-04 Thread Shengzhe Wu
Hello, I have a data set with 15 variables (first one is the response) and 1200 observations. Now I use pls package to do the plsr as below. trainSet = as.data.frame(scale(trainSet, center = T, scale = T)) trainSet.plsr = mvr(formula, ncomp = 14, data = trainSet, method = kernelpls,

Re: [R] R-square n p-value

2005-09-04 Thread Justin Rhodes
Dear Peter, This is exactly what I needed. The input is coming from the Mouse Phenome Project database (http://aretha.jax.org/pub-cgi/phenome/mpdcgi?rtn=docs/home) which only gives pearson's correlations r, and n. Thank you very much. I have used this R-help resource twice now recently

Re: [R] Doubt about nested aov output

2005-09-04 Thread Spencer Graves
Others may know the answer to your question, but I don't. However, since I have not seen a reply, I will offer a few comments: 1. What version of R are you using? I just tried superficially similar things with the examples in ?aov in R 2.1.1 patched and consistently got

Re: [R] Distributional characteristics

2005-09-04 Thread Spencer Graves
I haven't seen a reply to your post, and I would like to help you. Unfortunately, I don't see a question in your email. Please tell us why changing the variable (x-1,ln(x)) didn't get satisfying results, preferably using a toy example that someone else can copy from your email in to

Re: [R] R package for ICA

2005-09-04 Thread Spencer Graves
I just got 48 hits from RSiteSearch(independent component analysis), the first of which was http://finzi.psych.upenn.edu/R/library/mlica/html/mlica.html;. hope this helps. spencer graves p.s. I believe people who follow the posting guide typically get more useful

Re: [R] Question

2005-09-04 Thread Spencer Graves
1. I could find no references to var.ran with RSiteSearch(var.ran) and when requesting var.ran from S-Plus 6.2. 2. Have you considered simulate.lme? 3. Are you familiar with Pinheiro and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer)? I highly

Re: [R] how to fit the partial association model with R?

2005-09-04 Thread David Duffy
On Fri, 2 Sep 2005 [EMAIL PROTECTED] wrote: If I do not make a mistake,the partial association model is an extension of log-linear model.I read a papers which gives an example of it.(Sloane and Morgan,1996,An Introduction to Categorical Data Analysis,Annual Review of Sociology.22:351-375)

Re: [R] SpatStat Kest - Error Message help

2005-09-04 Thread Adrian Baddeley
On Thu, 1 Sep 2005, DrakeGis wrote: Hi I'm working with the function Kest in the package SpatStat (under LINUX with R 2.1.0). In order to evaluate the statistical significance of my point pattern I'm doing 999 Montecarlo replications. The script that use the Kest function runs OK for most of