search for: 0.321

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 =

Hi there, I am doing multiple comparisons for data that is not normally distributed. For this purpose I tried both functions kruskal{agricolae} and kruskalmc{pgirmess}. It confuses me that these functions do not yield the same results although they are doing the same thing, don't they? Can anyone tell my why this happens and which function I can trust? kruskalmc() tells me that there are no

Hello R users, Can any one please help me to find a way to build phylogenetic tree by R program from any distance matrix. Suppose I have a data like : MATRIX [1] '1' 0.0 [2] '2' 0.071 0.0 [3] '3' 0.036 0.286 0.0 [4] '4' 0.429 0.75 0.714 0.0 [5] '5' 0.679 0.179 0.214 0.536 0.0 [6] '6' 0.893 0.929 0.964 0.464 0.357 0.0 [7]

Dear subscribers, I am looking for a function which would allow me to model the dependent variable as the number of successes in a series of Bernoulli trials. My data looks like this ID TRIALS SUCCESSESS INDEP1 INDEP2 INDEP3 1 4444 0 0.273 0.055 0.156 2 98170 74 0.123 0.456 0.789 3 145486 30 0.124

Dear friends of nlme, Running quinModel (Pinheiro/Bates page 380) on R (current release, windows) gives: Nonlinear mixed-effects model fit by maximum likelihood Model: conc ~ quinModel(Subject, time, conc, dose, interval, lV, lKa, lCl) Data: Quinidine Log-likelihood: -497 Fixed: lV + lKa + lCl ~ 1 lV lKa lCl 5.382 -0.273 2.470 Random effects: Formula: list(lV ~ 1, lCl ~ 1)

Hi all, i have a data set cintaining 14 columns and 11 rows. Rows represent single point and columns represent the parameter measured. I wiuld like to compare the data to see which are more alike. I used the cluster analysis, but now i ma wondering if there are some other methods, since the cluster analysis did not give me the result i like. I tried to use factanal() for factor analysis

Hi, For further calculations I need to extract indivdual Variances of different random effects from a fitted model. I found out how to extract the correlations (VarCorr(m1)@reSumry$group1) but I was not able to find a way to extract the other components individually. To extract the Residuals I tried: (ranef(m1)@ stdErr) which unfortunately did not work. Thank you very much for your help!

Dear R-helpers, To get curves for a pseudo cohort other than the one centered at the mean of the covariates, I have been trying to use the newdata argument to survfit with no success. Here is my model statement, the newdata and the ensuing error. What am I doing wrong? > summary(fit) Call: coxph(formula = Surv(Start, Stop, Event, type = "counting") ~ Week + LagAOO + Prior.f +

Hi all, I have a dataframe of behavioral observations from 360 fish, each with 241 observation points(rows), which looks like this: > head(d) fish treatment tank trial video tid pid ang.chg abs.ac t len vel d2p x y 1 1 3 1 1 1 1 1 NA NA 0.0 0.000 NA NA 5.169 9.617 2

Dear all, can somebody please help me how to calculate "McFadden R-square" for a LOGIT model? Corresponding definition can be found here: http://publib.boulder.ibm.com/infocenter/spssstat/v20r0m0/index.jsp?topic=%2Fcom.ibm.spss.statistics.help%2Falg_plum_statistics_rsq_mcfadden.htm Here is my data: Data <- structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1,

Hello all, I am working with a 1000x1000 matrix, and I would like to return a 1000x1000 matrix that tells me which value in the matrix is greater than a theshold value (1 or 0 indicator). i have tried mat2<-as.matrix(as.numeric(mat1>0.25)) but that returns a 1:100000 matrix. I have also tried for loops, but they are grossly inefficient. THanks for all your help in advance. Lanre

I have the manual for S+ 6 and I'm trying to use R for the Principle Component Analysis example and I'm getting a few interesting answers... The log is as follows: R : Copyright 2001, The R Development Core Team Version 1.3.0 (2001-06-22) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type `license()' or

Hi, I have a Community (COM) composed of 6 species: A, B, C, D, E & F. The density of my Community is thus (Eq.1): dCOM = dA + dB + dC + dE + dF I would like to calculate and plot a linear regression between the density of each of my species and the density of the whole community (illustrating how the density of each species varies with variations of the whole community). For example, I would

Hi there, This problem has been dogging me for a bit, and I'm trying to figure out why. When running the the subsets function in the leaps library, R is giving me the following error message > lvodsub <- subsets(pred, resp$LVOD) Warning message: XHAUST returned error code -999 in: leaps.exhaustive(a, really.big = really.big) but this still happens if I add the really.big option:

I am trying to get confidence bands for a non-linear power function (y=mx^b). I thought I should be able to figure it out, but can't. Are there any R consultant? I would be willing to pay some amount of money, but not sure such consultants exist. I fit power functions to lots of data, and this would be very useful. I would ideally like to have confidence bands for the mean function and a

Dear all, I am fitting a LOGIT model on this Data........... Data <- structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 0, 1, 1, 0, 1, 0, 47, 58, 82, 100, 222, 164, 161, 70, 219, 81, 209, 182, 185, 104, 126, 192, 95, 245, 97, 177, 125, 56, 85, 199, 298, 145, 78, 144, 178, 146, 132, 98, 120, 148, 123, 282, 79, 34, 104,

Hello, I've been unsuccessfully trying to reproduce a sem from Grace et al. (2010) published in Ecological Monographs: http://www.esajournals.org/doi/pdf/10.1890/09-0464.1 The model in question is presented in Figure 8, page 81. The errors that I've been getting are: 1. Using a correlation matrix: res.grace <- sem(grace.model, S = grace, N = 190) Warning message: In sem.default(ram

Dear R helpers, I have been using the loadings function from the multiv library and I get the typical output (see below). When I try to export these results to a file using a write.table() I get the following error message "Error in as.data.frame.default(x[[i]], optional = TRUE) : can't coerce loadings into a data.frame" Any idea why write.table is doing that and any

I have a simple regression using lm(). If I just want to check the coefficient, I can use summary(lm())$coef; if I need the standard error, I can use summary(lm())$s, if I need the residuals, I can use summary(lm())$res. OK. How can I get the R-squares and Adjusted R-squares using $...? Is there a function, like objects(), that can show all the references for values? Thanks a lot! -- View this

#I am trying to understand how R fits models for contrasts in a #simple one-way anova. This is an example, I am not stupid enough to want #to simultaneously apply all of these contrasts to real data. With a few #exceptions, the tests that I would compute by hand (or by other software) #will give the same t or F statistics. It is the contrast estimates that R produces #that I can't seem to