similar to: pure statistical question

Hello! One note/question hier about specification of control-parameters in the lme(...,control=list(...)) function call: i tried to specify tne number of iteration needed via lme(....,control=list(maxIter=..., niterEM=...,msVerbose=TRUE)) but every time i change the defualt values maxIter (e.g. maxIter=1, niterEM=0) on ones specified by me, the call returns all the iterations needed until

Hi, I have been trying to work with error terms given back from aov to make confidence intervals. However, the numbers seem to be incorrect. If there is more than one term in the ANOVA then the error terms can be inflated by the number of factors in the extra terms. The F's are correct so it is right back to the SS. I was wondering if this is standard practice for stats programs or

Dear all, I would be glad if someone could tell me if Tukey HSD (function of stats library) accept the NA values and could correct the unbalanced design (different number of sample in each group) Thanks for help Alice Le Bars

I''m using a plugin which adds a method to ActiveRecord::Base http://similetimelinerailshelper.googlecode.com/svn/trunk/simile_timeline/ This adds an "acts_as" type declaration to the Models, via the file in vendor/plugins/simile_timeline/lib/simile_timeline.rb Rails is working with these declarations, and they are being used successfully in ''regular

The help for the round function states: round rounds the values in its first argument to the specified number of decimal places (default 0). Note that for rounding off a 5, the IEEE standard is used, ``go to the even digit''. Therefore round(0.5) is 0 and round(-1.5) is -2. Could someone provide a reference for this standard? Thanks, Anne ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Anne E.

You must create a output filter. Edit /etc/printcap and add your client printer like this: ... PC-PRINTER:\ :sd=/var/spool/lpd/PP01:\ :mx#0:\ :lp=/dev/null:\ :if=/var/spool/lpd/Filter/PP01:\ :sh: ... Adjust the printcap parameters as needed (spool directory, directory to find the filter script). Make sure the lpd can access the directory. Consult your manual. Now create the filter

Hi, I am trying to recreate the right graph on page 524 of Gelman's 2006 paper "Prior distributions for variance parameters in hierarchical models" in Bayesian Analysis, 3, 515-533. I am only interested, however, in recreating the portion of the graph for the overlain prior density for the half-Cauchy with scale 25 and not the posterior distribution. However, when I try:

Dear colleague: Actually that is what is done. When using the z-test between proportions in two different groups, or using chi-squared to test the null hypothesis of equal proportions of two or more groups, the null hypothesis is that H0: p1=p2=p3.....=p where this p is the pooled proportion, the proportion in the summed groups = n1+n2+n3+..../N1+N2+N3.... in the z-test between two

Hi, i am working on a data set with EDA. That includes QQ-Plots of residuals vs expected normal distribution. What puzzles me is that the range of ordinate and abscissae is so different: while the theoretical quantiles range from [-2, 2] the sample quantiles on the ordinate do extent from [-20, 50]. Quite obviously some kind of transformation is done. Although i intensively RTFM i could not

Is there a "canned" function in R for finding the standard error of the mean? I have tried > sem <- function(x) c(mean =mean(x), + SEM = stdev(x)/sqrt(length(x))) > sem(pnet.lai) Error in sem(pnet.lai) : couldn't find function "stdev" It looks like there is no stdev function in R Thanks, Kirk Kirk R. Wythers email: kwythers at umn.edu University of

Hello all, I am currently using NIS for authentication with clear text passwords. We would like to move to password encryption. Any ideas whether NIS would support encrypted password on samba 3.X ? if not what would be my other options. Sheikji Nazirudeen IT Analyst Syracuse University 315-443-1207

Hi I try to implement a simple 2-factorial repeated-measure anova in the lme framework and would be grateful for a short feedback -my dependent var is a reaction-time (rt), -as dependent var I have -the age-group (0/1) the subject belongs to (so this is a between-subject factor), and -two WITHIN experimental conditions, one (angle) having 5, the other 3 (hands) factor-levels;

Hi, I am getting some strange results using round - it seems that it depends if the number before the decimal point is odd or even .... For example: > round(1.5)[1] 2> round(2.5)[1] 2 While i would expect that round(2.5) be 3 and not 2. Do you have any explanation for that? I really appreciate your input, Monica

Full_Name: Jörg Polzehl Version: 1.5.1 OS: Windows 2000 Submission from: (NULL) (193.175.148.198) When calculating MLE's in a variance component model using constrained optimization, i.e. optim(...,method="L-BFGS-B",...) I observed an inproper behaviour in cases where the likelihood function was evalueted at the constraint. Parameters and value of the function at the constraint

I have question (curiosity) regarding returned values of R's qcauchy () function, for nonexceedance probability (F). It seems the ideal returned range of cauchy distribution should be [-Inf,Inf]. For F=0 > qcauchy(0) [1] -Inf but for F=1 > qcauchy(1) [1] 8.16562e+15 It seems to me that the proper return value should be Inf??? For default (location=0,scale=1) quantile function of

What is the current state of the R-art for ordered probit models, and more esoterically is there any available R strategy for ordered cauchit models, i.e. ordered multinomial alternatives with a cauchy link function. MCMC is an option, obviously, but for a univariate latent variable model this seems to be overkill... standard mle methods should be preferable. (??) Googling reveals that spss

Dear all, I have the following problem. Given an expression object 'expr' containing a certain set of symbols (say 'a', 'b', 'c'), I would like to translate the expression object in an R function of, say, 'a', programmatically. Here an example of what I mean. Given: > expr <- expression(a+b+c) a call like: > asFunctionOf(expr, 'a',

I just created a new smbpasswd file, all users with unset passwords. Problem is, users can't set their passwords, only I can as root. Password is not set / blank, but they still get: >test> smbpasswd >Old SMB password: <enter> >New SMB password: Samba >Retype new SMB password: Samba >machine 127.0.0.1 rejected the password change: Error was : RAP86: The >specified

I am in a situation where I have to fit a distrution, such as cauchy or normal, to an empirical dataset. Well and good, that is easy. But I wanted to assess just how good the fit is, using ks.test. I am concerned about the following note in the docs (about the example provided): "Note that the distribution theory is not valid here as we have estimated the parameters of the normal

Dear friends, Hope you are all doing well and having a great weekend. I have data that was collected on specific gravity and spectrophotometer analysis for 26 mixtures of NG (nitroglycerine), TA (triacetin), and 2 NDPA (2 - nitrodiphenylamine). In the dataset, x1 = %NG, x2 = %TA, and x3 = %2 NDPA. The response variable is the specific gravity, and the rest of the variables are the predictors.