similar to: by function does not separate output from function with mulliple parts

Displaying 20 results from an estimated 100 matches similar to: "by function does not separate output from function with mulliple parts"

2010 Apr 05
1
new to R, analysis of latency data
Hi, I'd like to move from excel to R because our dataset are so large. Here's what my data looks like: Transaction Rate Run# Transaction Type Location Latency in Seconds 10 1 Order A 0 10 1 Order B
2003 Aug 19
3
On the Use of the nnet Library
Dear List, I am trying to solve a problem by the neural network method(library: nnet). The problem is to express Weight in terms of Age , Sex and Height for twenty people. The data frame consists of 20 observations with four variables: Sex, Age, Height and Weight. Sex is treated as a factor, Age and Weight are variables normalized to unity, as usual. I wanted to construct a neural network, and so
2009 Apr 15
2
AICs from lmer different with summary and anova
Dear R Helpers, I have noticed that when I use lmer to analyse data, the summary function gives different values for the AIC, BIC and log-likelihood compared with the anova function. Here is a sample program #make some data set.seed(1); datx=data.frame(array(runif(720),c(240,3),dimnames=list(NULL,c('x1','x2','y' )))) id=rep(1:120,2); datx=cbind(id,datx) #give x1 a
2012 Jun 02
8
Error
Hello list, after somes days I have this error when I run git pull in samba4 error: Unable to find 3b6a88a1a49fdc23cd956c6583e9ec56e541eeae under http://gitweb.samba.org/samba.git Cannot obtain needed blob 3b6a88a1a49fdc23cd956c6583e9ec56e541eeae while processing commit 1a596076697624545f415eb558c4647267cfe46c. error: Fetch failed. root at jupiter:/home/run# error: Unable to find
2019 Apr 24
1
Bug in "stats4" package - "confint" method
Dear R developers, I noticed a bug in the stats4 package, specifically in the confint method applied to ?mle? objects. In particular, when some ?fixed? parameters define the log likelihood, these parameters are stored within the mle object but they are not used by the ?confint" method, which retrieves their value from the global environment (whenever they still exist). Sample code: >
2018 Jan 17
1
Assessing calibration of Cox model with time-dependent coefficients
I am trying to find methods for testing and visualizing calibration to Cox models with time-depended coefficients. I have read this nice article <http://journals.sagepub.com/doi/10.1177/0962280213497434>. In this paper, we can fit three models: fit0 <- coxph(Surv(futime, status) ~ x1 + x2 + x3, data = data0) p <- log(predict(fit0, newdata = data1, type = "expected")) lp
2011 Oct 06
1
anova.rq {quantreg) - Why do different level of nesting changes the P values?!
Hello dear R help members. I am trying to understand the anova.rq, and I am finding something which I can not explain (is it a bug?!): The example is for when we have 3 nested models. I run the anova once on the two models, and again on the three models. I expect that the p.value for the comparison of model 1 and model 2 would remain the same, whether or not I add a third model to be compared
2012 Nov 15
1
Step-wise method for large dimension
Hi , I want to apply the following code fo my data with 400 predictors. I was wondering if there ia an alternative way instead of typing 400 predictors for the following code. I really appreciate your help. fit0<-lm(Y~1, data= mydata) fit.final<- lm(Y~X1+X2+X3+.....+X400, data=mydata) ??? step(fit0, scope=list(lower=fit0, upper=fit.final), data=mydata, direction="forward")
2018 Jan 18
1
Time-dependent coefficients in a Cox model with categorical variants
First, as others have said please obey the mailing list rules and turn of First, as others have said please obey the mailing list rules and turn off html, not everyone uses an html email client. Here is your code, formatted and with line numbers added. I also fixed one error: "y" should be "status". 1. fit0 <- coxph(Surv(futime, status) ~ x1 + x2 + x3, data = data0) 2. p
2009 May 11
1
Warning trying to plot -log(log(survival))
windows xp R 2.8.1 I am trying to plot the -log(log(survival)) to visually test the proportional hazards assumption of a Cox regression. The plot, which should give two lines (one for each treatment) gives only one line and a warning message. I would appreciate help getting two lines, and an explanation of the warning message. My problem may the that I have very few events in one of my strata,
2008 Nov 26
1
Request for Assistance in R with NonMem
Hi I am having some problems running a covariate analysis with my colleage using R with the NonMem program we are using for a graduate school project. R and NonMem run fine without adding in the covariates, but the program is giving us a problem when the covariate analysis is added. We think the problem is with the R code to run the covariate data analysis. We have the control stream, R code
2011 May 08
1
anova.lm fails with test="Cp"
Here is an example, modified from the help page to use test="Cp": -------------------------------------------------------------------------------- > fit0 <- lm(sr ~ 1, data = LifeCycleSavings) > fit1 <- update(fit0, . ~ . + pop15) > fit2 <- update(fit1, . ~ . + pop75) > anova(fit0, fit1, fit2, test="Cp") Error in `[.data.frame`(table, , "Resid.
2004 Apr 27
0
Extracting labels for residuals from lme
Dear R-helpers, I want to try to extract residuals from a multi-level linear mixed effects model, to correlate with another variable. I need to know which residuals relate to which experimental units in the lme. I can show the labels that relate to the experimental units via the command ranef(fit0)$resid which gives: 604/1/0 -1.276971e-05 604/1/1 -1.078644e-03 606/1/0 -7.391706e-03 606/1/1
2009 Jun 15
2
coxph and robust variance estimation
Hello, I would like to compare two different models in the framework of Cox proportional hazards regression models. On Rsitesearch and google I don't find a clear answer to my question. My R-Code (R version 2.9.0) coxph.fit0 <- coxph(y ~ z2_ + cluster(as.factor(keys))+ strata(stratvar_), method="breslow" ,robust=T ) coxph.fit1 <- coxph(y ~ z_ +
2009 May 10
2
plot(survfit(fitCox)) graph shows one line - should show two
R 2.8.1 Windows XP I am trying to plot the results of a coxph using plot(survfit()). The plot should, I believe, show two lines one for survival in each of two treatment (Drug) groups, however my plot shows only one line. What am I doing wrong? My code is reproduced below, my figure is attached to this EMail message. John > #Create simple survival object >
2019 Dec 27
0
"simulate" does not include variability in parameter estimation
On 26/12/2019 11:14 p.m., Spencer Graves wrote: > Hello, All: > > > ????? The default "simulate" method for lm and glm seems to ignore the > sampling variance of the parameter estimates;? see the trivial lm and > glm examples below.? Both these examples estimate a mean with formula = > x~1.? In both cases, the variance of the estimated mean is 1. That's how
2006 Feb 27
1
help with step()
Folks: I'm having trouble doing a forward variable selection using step() First, I fit an initial model: fit0 <- glm ( est~1 , data=all, subset=c(n>=25) ) then I invoke step(): fit1 <- step( fit0 , scope=list(upper=est~ pcped + pchosp + pfarm ,lower=est~1)) I get the error message: Error in eval(expr, envir, enclos) : invalid 'envir' argument I looked at the
2006 Oct 18
1
MARS help?
I'm trying to use mars{mda} to model functions that look fairly close to a sequence of straight line segments. Unfortunately, 'mars' seems to totally miss the obvious places for the knots in the apparent first order spline model, and I wonder if someone can suggest a better way to do this. The following example consists of a slight downward trend followed by a jump up after
2019 Dec 27
1
"simulate" does not include variability in parameter estimation
On 2019-12-27 04:34, Duncan Murdoch wrote: > On 26/12/2019 11:14 p.m., Spencer Graves wrote: >> Hello, All: >> >> >> ? ????? The default "simulate" method for lm and glm seems to ignore the >> sampling variance of the parameter estimates;? see the trivial lm and >> glm examples below.? Both these examples estimate a mean with formula = >>
2007 Jul 14
0
memory problem
Hi, My computer has 2GB of ram and I also request 2GB of virtual ram from c drive, therefore totally I have 4GB of ram. Before I open R workshop, I also add "C:\Program Files\R\R-2.5.0\bin\Rgui.exe" --max-mem-size=3000Mb--max-vsize=3000Mb" into the target of R by right clicking the R icon-properties. I am running a program (DPglmm), which is built in DPpackage. The program is