Displaying 20 results from an estimated 100 matches similar to: "nlme Output"
2006 May 01
2
Pasting data into scan()
The file TENSILE.DAT from the Hand et al "Handbook of Small Data Sets"
looks like this:
0.023 0.032 0.054 0.069 0.081 0.094
0.105 0.127 0.148 0.169 0.188 0.216
0.255 0.277 0.311 0.361 0.376 0.395
0.432 0.463 0.481 0.519 0.529 0.567
0.642 0.674 0.752 0.823 0.887 0.926
except that my mail client has replaced the tab separators by blanks. If
I paste this data into R 2.2.1 what I get is
2009 Jan 14
1
Can't find files after install package (Windows)
I am a total newbie at R but experienced with computers. If this is not the
right forum for this question, please let me know one that is. I searched in
the R manuals in the Help section, the Nabble help and Rseek but no luck.
I am trying to install the package gcl on my recently installed R 2.8.1
I searched on cran for the windows version and found a zip file for the
package
In RGui, I ran
2017 Oct 21
2
Help_urgent_how to calculate mean and sd in biomod 2
Hello
I am new in R. I am trying to implement Biomod 2 package.
However, I have a doubt. I want to calculate the mean and sd of
"Testing.data"
(ROC and TSS)
> # let's print the ROC scores of all selected models
> myBiomodModelEval_55["ROC","Testing.data",,,]
RUN1 RUN2 RUN3 RUN4 RUN5 RUN6 RUN7 RUN8 RUN9 RUN10
0.938 0.938 0.926 0.931 0.939
2005 Nov 17
1
Principal Components Analysis (PR#8320)
Full_Name: Sahotra Sarkar
Version: 2.2.0
OS: Windows XP Professional
Submission from: (NULL) (146.6.130.180)
The following two commands should give the same results for the eigenvectors but
do not (there is a sign reversal for the first one):
> summary(princomp(bumpus),loading = TRUE)
Importance of components:
Comp.1 Comp.2 Comp.3 Comp.4 Comp.5
2011 Jun 08
2
Results of CFA with Lavaan
I've just found the lavaan package, and I really appreciate it, as it
seems to succeed with models that were failing in sem::sem. I need
some clarification, however, in the output, and I was hoping the list
could help me.
I'll go with the standard example from the help documentation, as my
problem is much larger but no more complicated than that.
My question is, why is there one latent
2011 Apr 05
6
simple save question
Hi,
When I run the survfit function, I want to get the restricted mean
value and the standard error also. I found out using the "print"
function to do so, as shown below,
print(km.fit,print.rmean=TRUE)
Call: survfit(formula = Surv(diff, status) ~ 1, type = "kaplan-meier")
records n.max n.start events *rmean *se(rmean) median
200.000
2012 Oct 22
1
glm.nb - theta, dispersion, and errors
I am running 9 negative binomial regressions with count data.
The nine models use 9 different dependent variables - items of a clinical
screening instrument - and use the same set of 5 predictors. Goal is to
find out whether these predictors have differential effects on the items.
Due to various reasons, one being that I want to avoid overfitting models,
I need to employ identical types of
2012 Jul 02
1
How to get prediction for a variable in WinBUGS?
Dear all,I am a new user of WinBUGS and need your help. After running the following code, I got parameters of beta0 through beta4 (stats, density), but I don't know how to get the prediction of the last value of h, the variable I set to NA and want to model it using the following code.Does anyone can given me a hint? Any advice would be greatly appreciated.Best
2012 Jun 02
2
mgcv (bam) very large standard error difference between versions 1.7-11 and 1.7-17, bug?
Dear useRs,
I reran an analysis with bam (mgcv, version 1.7-17) originally
conducted using an older version of bam (mgcv, version 1.7-11) and
this resulted in the same estimates, but much lower standard errors
(in some cases 20 times as low) and lower p-values. This obviously
results in a larger set of significant predictors. Is this result
expected given the improvements in the new version? Or
2008 Mar 04
2
Asking, are simple effects different from 0
Hello, R-i-zens. I'm working on an data set with a factorial ANOVA
that has a significant interaction. I'm interested in seeing whether
the simple effects are different from 0, and I'm pondering how to do
this. So, I have
my.anova<-lm(response ~ trtA*trtB)
The output for which gives me a number of coefficients and whether
they are different from 0. However, I want the
2011 Jul 25
4
ggplot question: changing the label for the Y axis on a histogram
Some help with how to re-label the vertical axis in a histogram would be appreciated.
qplot(off.sc,weight=rel.freq,binwidth=.29,main="test Figure"+ylab("New from inside"))+ylab("New from outside")+
xlab("off.sc\nAggregated frequency plots for 17 equal intervals.")
The code
2004 Jan 15
1
nlme vs aov with Error() for an ANCOVA
Hi
I compouted a multiple linear regression with repeated measures on one
explanatory variable:
BOLD peak (blood oxygenation) as dependent variable,
and as independent variables I have:
-age.group (binaray:young(0)/old(1))
-and task-difficulty measured by means of the reaction-time 'rt'. For
'rt' I have repeated measurements, since each subject did 12 different
tasks.
-> so
2010 May 28
4
vlookup in R?
Hi R-users,
I would like to search for the values of seq that match my rand values. In excel I will use =VLOOKUP(G2,$E$2:$F$32,2). For example, for rand=.262 it will give me approximately seq=120 and rand=0.964293344, seq=460 and etc.
E F G
cdf seq rand
0.00E+00 0 0.262123478
1.56E-03 20 0.964293344
1.55E-02 40 0.494827113
5.30E-02 60
2006 Apr 06
1
interpreting anova summary tables - newbie
Hello,
Apologies if this is the wrong list, I am a first-time poster here. I
have an experiment in which an output is measured in response to 42
different categories.
I am only interested which of the categories is significantly different
from a reference category.
Here is the summary of the results:
summary(simple.fit)
Call:
lm(formula = as.numeric(as.vector(TNFa)) ~ Mutant.ID, data =
2006 Mar 08
1
RES: survival
Dear Thomas,
The head of my dataset
> head(wsuv)
parcel sp time censo treatment
species
1 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1 1
2 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1 1
3 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1 1
4 S8 Poecilanthe effusa ( Hub. ) Ducke. 1 1 1
2011 Sep 26
3
survival analysis: interval censored data
hello:
my data looks like:
time1 time2 event catagoria
2004 2006 1 C
2004 2005 0 C
2005 2010 1 E
2007 2009 1 C
2006 2007 0 E
2008 2010 0 C
2008 2010 1 E
...
and the census interval is 1 year
I have tried this
2018 May 15
0
Systemfit
... and the mailing list is picky about attachments... whatever you attached did not conform to the stringent requirements mentioned in the Posting Guide. Pasting the code right into the email is usually safest, though you DO have to post using plain text (as the Posting Guide indicates) or your code may get mangled by the automatic html format removal.
On May 15, 2018 7:04:31 AM PDT, Bert Gunter
2018 May 15
2
Systemfit
OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel.
#Reproducible script
Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv")
View(Empdata)
install.packages("systemfit")
2018 May 15
1
Systemfit
Unless there is good reason not to, always cc the list -- there are lots of
smarter folks than I on it who can help.
I may or may not have time to look at this. Hopefully someone else will.
-- Bert
Bert Gunter
"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip
2005 Jan 30
1
x86_64 timer resolution in ping results
Hi,
There seems to be a problem with the timer resolution on x86_64. This is
clear in ping results - I don't know if it's just a ping problem or if
it manifests elsewhere as well. It does also affect mtr and traceroute
times.
I guess it's more a kernel issue than a CentOS problem, but is anyone
else seeing the same?
Ping results have a resolution of 10ms - so a ping <10ms