similar to: Fitting t-Student(mu, sigma, nu)

Hello everybody, Can anyone tell me how I could analyze data that are at a contigency table form? I already found function cfa in the cfa package but I still don't understand how I could use this function in order to elaborate a contigency table. Every answer is welcome! --------------------------------- ÁðïêôÞóôå ôçí äùñåÜí óáò@yahoo.gr [[alternative HTML version deleted]]

Hello, I have 4 questions in R and I am looking for 4 answers! the questions are the below: 1) is there a function in R for the Least Significant Difference method for multiple comparisons? 2) Is there a function in R that gives me the studentized deleted residuals and the leverage points? 3) is there a way that I could create a link in a program of mine, so as to select an option when my program

Hello everybody, I've been practicing with some data in R and SPSS and I noticed that there are some differences in ANOVA results. For example with : y<-c(1,2,34,2,3,45,2,1,67,3,2,67,2,2,98,4,4,23,1,1,23,2,3,45) and x<-rep(c(1,2,3),8) I get in R ( with summary(aov(y~x)) ) MSres=350.7 df=22 while in SPSS I get MSres=221.9 df=21 Can enyone explain me what is the problem or

Hello R users, I have three questions and I would be grateful if someone could give me an answer to each of these. 1) I have constracted a function that returns an output, which runs in a while( condition ){ run function } loop. I would like to know if there is a way to get the outputs in different windows, every time the function runs, so as to compare easier the results. 2) In my

hello, I tested all config with my ups and psp all seems to be OK now, I had change ups name and user name in ups.conf, upsd.users, upsmon.conf etc etc ... but I have quickly these events in my syslog without doing anyting what can I do ? Nov 15 11:09:22 mezenc upsmon[3865]: Poll UPS [pm2200 at localhost] failed - Data stale thanks a lot for your help Sabine GOUDARD Ecole Nationale

Greetings, I am in the process of moving several large sites to a Tyan+BQ+CentOS box (built from the Nuonce ISO) and I see that special characters are not showing up in the browser properly but do show up in the source "view source" properly. -- very strange. Example: on the web page --> Heizungsbau M?ller in the source --> Heizungsbau M?ller When I view this on the existing

#bug 1: besselI() for nu<0 and expon.scaled=TRUE #tested with R-devel (2007-06-17 r41981) x <- 2.3 nu <- -0.4 print(paste(besselI(x, nu, TRUE), "=", exp(-x)*besselI(x, nu, FALSE))) #fix: #$ diff bessel_i_old.c bessel_i_new.c #57c57 #< bessel_k(x, -alpha, expo) * ((ize == 1)? 2. : 2.*exp(-x))/M_PI #--- #> bessel_k(x, -alpha, expo) * ((ize == 1)? 2. :

Hi, I am trying to do a PCA on my data but I keep getting the error message svd(x, nu=0) infinite or missing values >From the messages posted on the subject, I understand that the NAs in my data might be the problem, but I thought na.omit would take care of that. Less than 5% of my cells are missing data. However, the NAs are not regularly distributed across my matrix: certain cases and

I'm trying to get started with maximum likelihood estimation with a simple regression equivalent out of Bolker (Ecological Models and Data in R, p302). With this code: #Basic example regression library(bbmle) RegData<-data.frame(c(0.3,0.9,0.6),c(1.7,1.1,1.5)) names(RegData)<-c("x", "y") linregfun = function(a,b,sigma) { Y.pred = a+b*x

Dear List, I am trying to label a plot with the symbol +/- sigma. Using something like - expression (2*sigma) gives me the symbol 2ó. However, adding +/- to it beats me. The code I am using is: plot(x,y,type="l",main=" expression(paste("±", plain(2*ó)),sep=""). Any suggestion will be appreciated. Best Ogbos [[alternative HTML version deleted]]

Hi, I'm playing around with gamlss and don't entirely understand the sigma result from an attempted lognormal fit. In the example below, I've created lognormal data with mu=10 and sigma=2. When I try a gamlss fit, I get an estimated mu=9.947 and sigma=0.69 The mu estimate seems in the ballpark, but sigma is very low. I get similar results on repeated trials and with Normal and

Dear all &nbsp; I wish to carry out sigma- and beta-convergence analysis in respect of panel data [wherein current value of one of the variables needs be regressed upon suitably transformed lagged values of another variable(s)]. I am quite new to the R-language and am not very much aware of the availbaility of suitable package(s)/ code in the language. Can any one help me in letting me know of

Again, Hi All, Thanks for your replies. I probably should've mentioned that I'm not subscribed to the list but I found your answers on the archives. I just wanted to go on record and mention that I do in fact have all the original CD's and that by re-encoding I meant repeating the entire process of ripping the original CD's (with cdparanoia (of course ;)) and then encoding them

Dear Sir, I have data, coming from tests, consisting of 300 values. Is there a way in R with which I can confirm this data to 68-95-99.8 rule or three-sigma rule? I need to look around percentile ranks and prediction intervals for this data. I, however, used SixSigma package and used ss.ci() function, which produced 95% confidence intervals. I still am not certain about percentile ranks

I'm proposing to add something like this to the stats package : ---------------------------------------------------------- ### "The" sigma in lm/nls - "like" models: sigma <- function(object, ...) UseMethod("sigma") ## works whenever deviance(), nobs() and coef() do fine: sigma.default <- function (object, use.fallback=TRUE, ...)

Dear R community, I have a data matrix (531X314), and would like to apply the prcomp. However, I got this error Lapack message. I am using R3.2.2. I googled a bit and found that it might be related to converge issue. ?Just wonder if there is a way to get around it? Thank you very much! Ace On Thursday, December 29, 2016 11:44 AM, Ista Zahn <istazahn at gmail.com> wrote: Use

How might you fit a generalized linear model (glm) with variance = mu+theta*mu^2 (where mu = mean of the exponential family random variable and theta is a parameter to be estimated)? This appears in Table 2.7 of Fahrmeir and Tutz (2001) Multivariate Statisticial Modeling Based on Generalized Linear Models, 2nd ed. (Springer, p. 60), where they compare "log-linear model fits to

Dear list, I was trying to create a VarClass for nlme to work with Fay-Herriot (FH) models. The idea was to create a modification of VarComb that instead of multiplying the variance functions made their sum (I called it varSum). After some fails etc... I found that the I was not getting the expected results because I needed to make sigma fixed. Trying to find how to make sigma fixed I run into

Hello I have run an anova analysis for the fallowing model: H_obs=mu+REGION+MANAGEMENT + e When I run it in ASRelm I get the p-value for mu, and, of course also for the two dependent variables (REGION and MANAGEMENT) When I run it in R, I do not get the pvalue for mu. Can some one help me to understand why? and if it is possible to estimate the pvalue for mu in anova() in R? I attach the

Hi Folks, I am interested in extending the repertoire of link functions in glm(Y~X, family=binomial(link=...)) to include a "tan" link: eta = (4/pi)*tan(mu) i.e. this link bears the same relation to the Cauchy distribution as the probit link bears to the Gaussian. I'm interested in sage advice about this from people who know their way aroung glm. >From the surface, it looks