similar to: NLME error model with several responses

Hello all (especially MCLUS users). I'm trying to make use of the MCLUST package by C. Fraley and A. Raftery. My problem is trying to figure out how the (model) identifier (e.g, EII, VII, VVI, etc.) relates to the covariance matrix. The parameterization of the covariance matrix makes use of the method of decomposition in Banfield and Rraftery (1993) and Fraley and Raftery (2002) where

Hi everybody! I want to find a closer neighbourins observation. This is my code: ########################## library(klaR) library(ipred) library(mlbench) data(PimaIndiansDiabetes2) dane=na.omit(PimaIndiansDiabetes2)[,c(2,5,9)] dane[,2]=log(dane[,2]) dane[,1:2]=scale(dane[,1:2]) zbior.uczacy=sample(1:nrow(dane),nrow(dane)/2,F)

Hello, I have a problem, when I run read.fasta.pdb I get this error: pdb/seq: 1 name: dcluster.1 Error: subscript out of bounds My FASTA file is the sequence of the bovine insuline. Thanks in advance, Rg -- View this message in context: http://n4.nabble.com/problem-with-bio3d-package-tp1012116p1012116.html Sent from the R help mailing list archive at Nabble.com.

Hi I want to create a nice picture about my result of k -the nearest neighbours algorithm. Here is my easy code: ################################# library(klaR) library(ipred) library(mlbench) data(PimaIndiansDiabetes2) dane=na.omit(PimaIndiansDiabetes2)[,c(2,5,9)] dane[,2]=log(dane[,2]) dane[,1:2]=scale(dane[,1:2]) zbior.uczacy=sample(1:nrow(dane),nrow(dane)/2,F)

Hi there, >library(faraway) >pima pregnant glucose diastolic triceps insulin bmi diabetes age test 1 6 148 72 35 0 33.6 0.627 50 1 2 1 85 66 29 0 26.6 0.351 31 0 >data(pima) >pima pregnant glucose diastolic triceps insulin bmi diabetes age test 1 6 148 72 35 0 33.6

I just installed 302p1 on a client's machine (I run it on all of mine) and something is definetly wrong, can't figure out what. ssh localhost ssh_exchange_identification: Connection closed by remote host stock config files except ssh allows X11 forwarding any suggestions? Michael at Insulin-Pumpers.org

Unsubscribe: http://mail4.haemoundheilung.me/unsubscribe.php?M=229682&C=740f45a9d08a410ee15403710005cd65&L=1&N=2 Forget anything you`ve ever been told about Diabetes... Because this `average` man`s discovery has not only CURED over 20.000 people of their diabetes in just two weeks, but it has the potential to wipe out the $245 billion diabetes pharma industry FOR GOOD. And get this

Hi, I'm having a bit of trouble with the axes in my plots. I don't like the way R does not have them cross in the origin. Is there another plot/axis function? i tried using abline as suggested by someone from this list, but in my case this gives no satisfactory result, as the abline does sometimes lie on top of the y axis and sometimes not, depending on how i scale the image... #bild16

Hi all, I have some difficulties with the lme function and so this is my problem. Supoose i have the following model y_(ijk)=beta_j + e_i + epsilon_(ijk) where beta_j are fixed effects, e_i is a random effect and epsilon_(ijk) is the error. If i want to estimate a such model, i execute >lme(y~vec.J , random~1 |vec .I ) where y is the vector of my data, vec.J is a factor object

Hello, I have a simulated data structure in which students are nested within teachers, and with each student are associated two test scores. There are 20 classrooms and 25 students per classroom, for a total of 500 students and two scores per student. Here are the first 10 lines of my dataframe "d": studid tchid Y time 1 1 1 -1.0833222 0 2 1 1

Hi R-helpers, I'm new to time series modelling, but my requirement seems to fall just outside the capabilities of the arima function in R. I'd like to fit an ARMA model where the variance of the disturbances is a function of some exogenous variable. So something like: Y_t = a_0 + a_1 * Y_(t-1) +...+ a_p * Y_(t-p) + b_1 * e_(t-1) +...+ b_q * e_(t-q) + e_t, where e_t ~ N(0, sigma^2_t),

Folks Has anyone any experience with these? http://www.automation.siemens.co.uk/main/business%20groups/power/ups/products/ups%20systems/single%20phase/ http://www.automation.siemens.co.uk/main/business%20groups/power/ups/products/ups%20systems/single%20phase/Siemens%20leaflet%20UPS%20Series%20A.pdf They look a very solid machine, I've got a chance to pick one up for NZ$50 with an extra

Dear R-users Coming from a proc mixed (SAS) background I am trying to get into the use of (n)lme. In this connection, I have some (presumably stupid) questions which I am sure someone out there can answer: 1) With proc mixed it is easy to get a hold on the estimated variance parameters as they can be put out into a SAS data set. How do I do the same with lme-objects? For example, I can see the

Dear All, I have a model to predict time series data for example: data(LakeHuron) Lake.fit <- arima(LakeHuron,order=c(1,0,1)) then the function predict() can be used for predicting future data with the model: LakeH.pred <- predict(Lake.fit,n.ahead=5) I can see the result LakeH.pred$pred and LakeH.pred$se but I did not see residual in predict function. If I have a model: [\ Z_t =

Dear List, I'm just teaching myself semi-parametric techniques. Apologies in advance for the long post. I've got observational data and a longitudinal, semi-parametric model that I want to fit in GAM (or potentially something equivalent), and I'm not sure how to do it. I'm posting this to ask whether it is possible to do what I want to do using "canned" commands

---------- Forwarded message ---------- From: ma yuchao <ma.yuchao@gmail.com> Date: 2006-5-20 ÉÏÎç4:01 Subject: hello, everyone To: R-help@stat.math.ethz.ch Hello, R people: I have a question in using fSeries package--the funciton garchFit and garchOxFit if adding a regression to the mean formula, how to estimate the model in R? using garchFit or garchOxFit? For example,

>From what I understand, the problem is due to people's disagreement about what the "correct" behavior should be. I'm pretty sure that the following is the correct behavior from running rsh and ssh often (both fsecure and openssh). Lets say you have a stupid script that does while 1 do sleep 1 done Called foreverSleep on your remote host: rsh remotehost

Dear all I need to run an anova from a factorial model y_{ijk}=\alpha_i+\beta_j+(\alpha\beta)_{ij}+e_{ijk} and calculate type II and III sums of square, but I have an empty cells, so anova function from package car fail. (I believe) y<-c(7,13,6,10,8,11,8,3,7,5,65) a<-as.factor(c(1,1,2,2,3,3,3,1,1,1,2)) b<-as.factor( c(rep(1,7),rep(2,4)) ) table(b,a) # cell (2,3) is empty

> Ady wrote: > > This is a quote from Wonko/Jaclaz, who has also been investigation this > > issue and deserves credit for it: > > Ok, not that I understand the code, but what *somehow* happens is that > > the "embedded" checksum in Isolinux.bin (starting from 4.00) is the > > checksum of the WHOLE file (i.e. starting from offset 0 instead of 64) >

> > E_ Under FreeDOS, SYSLINUX.COM doesn't seem to install ldlinux.* > > (specially when used with "--directory"), and there are no ([error]) > > messages; not even when executed with no parameters, or with > > "--help", or with "--version". The messages disappeared _after_ > > 5.00-pre7. > > Right, as Peter pointed out on