Displaying 20 results from an estimated 200 matches similar to: "problems with for loop"
2007 Jan 29
1
Fwd: Re: LSD multiple comparison test
I am returning this to the R-help list. Please keep followups on the list.
Yes, it can be done. It is not currently easy because multcomp doesn't
have the syntax yet. Making this easy is on Torsten's to-do list for the
multcomp package.
See the MMC.WoodEnergy example in the HH package. The current version on CRAN
is HH_1.17. Please see the discussion of this example in R-help:
2007 Jan 29
3
Multiple comparisons when interacction
In the model:
lm.1 <- lm(variable ~ BLOC + TIL * YEAR , data=selvanera)
I found TIL*YEAR interaction significant. Then I am trying to compare
means of the different levels of TIL inside every YEAR using:
mc.2 <- glht(lm.1, linfct = mcp(TIL*YEAR="Tukey"))
summary(mc.2, test = univariate())
but it does not work.
There is any way of doing this, like the SLICE option in
2009 May 11
1
predict function the other way around
Dear List,
Consider the following example
x=c(1,2,3,4,5)
y=c(2,4,6,8,10)
linearmodel=lm(y~x)
To predict a y-value if you know the corresponding x value is very easy
with the command predict.
predict(linearmodel, newdata=(x=1.5))
The other way around, to predict an x-value with a corresponding
y-value, doesn't work unfortunately.
Is there another function that can do that, or do I need to
2009 Jul 15
2
storing lm() results and other objects in a list
to clean up some code I would like to make a list of arbitrary length
to store?various objects for use in a loop
sample code:
############ BEGIN SAMPLE ##############
# You can see the need for a loop already
linearModel1=lm(modelSource ~ .,mcReg)
linearModel2=step(linearModel1)
linearModel3=lm(modelSource ~ .-1,mcReg)
linearModel4=step(linearModel3)
#custom
linearModel5=lm(modelSource ~ .
2009 Nov 09
1
Getting Sphericity Tests for Within Subject Repeated Measure Anova (using "car" package)
Hello everyone,
I am trying to do within subjects repeated measures anova followed by the
test of sphericity (sample dataset below).
I am able to get either mixed model or linear model anova and TukeyHSD, but
have no luck with Repeated-Measures Assuming Sphericity or Separate
Sphericity Tests.
I am trying to follow example from "car" package, but it seems that I am not
getting something
2009 Jul 21
2
Odd coefficent behavior
Why are my coefficients getting appended with a 1? It borks a match I
do later against the original list that doesn't have the random 1
added to the end.
> linearModel[[1]]
Call:
lm(formula = modelSource ~ +UNITBUILD + UNITDB + ITBUILD + ITDB +
UATBUILD + UATDB + HOGANCODE + RCF + ReleaseST1 + ReleaseST2 +
ReleaseBLA + Small.Bank.Acquisitions + HLY.NewYear + HLY.MLK +
HLY.PRES +
2009 Nov 09
1
Getting Sphericity Tests for Within Subject Repeated Measure Anova (using "car" package) (Adjusted Dataset)
[corrected dataset below]
Hello everyone,
I am trying to do within subjects repeated measures anova followed by the
test of sphericity (sample dataset below).
I am able to get either mixed model or linear model anova and TukeyHSD, but
have no luck with Repeated-Measures Assuming Sphericity or Separate
Sphericity Tests.
I am trying to follow example from "car" package, but it seems
2009 Dec 03
0
Problem with predict() and factors
I am working on a script that takes numeric performance indicators and runs
them against a series of regressors (dummy regressors, yes\no stuff via 0
and 1, e.g. Was is Christmas this week 0=no, 1=yes).
The script is as follows (Written as a function):
-- Begin Script --
doEnv <- function(HOUR,ENVNAME,REPORTNAME) {
library(RODBC)
library(forecast)
library("geneplotter")
2016 Apr 04
0
Test for Homoscedesticity in R Without BP Test
Hi Deepak,
In econometrics there is another test very often used : the white test.
The white test is based on the comparison of the estimated variances of residuals when the model is estimated by OLS under the assumption of homoscedasticity and when the model is estimated by OLS under the assumption of heteroscedastic.
The White test with R
install.packages("bstats")
library(bstats)
2016 Apr 04
1
Test for Homoscedesticity in R Without BP Test
On Mon, 4 Apr 2016, varin sacha via R-help wrote:
> Hi Deepak,
>
> In econometrics there is another test very often used : the white test.
> The white test is based on the comparison of the estimated variances of
> residuals when the model is estimated by OLS under the assumption of
> homoscedasticity and when the model is estimated by OLS under the
> assumption of
2003 Aug 07
0
new version of Rcmdr package
Dear list members,
I've uploaded a new version of the Rcmdr package to CRAN. There are many
additions and (I hope) improvements, as indicated in the relevant portion
of the CHANGES file, reproduced below. As usual, comments, suggestions, and
bug reports are appreciated.
John
-------------------------------------------
Version 0.8-4
o Minor bug fixes and additional input-error
2006 Mar 11
1
A question on plot
I’m trying to obtain a plot of given function, say the density of a
normal, in which I could indicate the area under the curve between two
points (a,b).
I.e. the P(a<X<b). Could anyone indicate me which is the instruction for
obtain such a shadow area under a curve?
I guess the solution is not difficult, but I can not find the procedure.
Thanks
Albert Sorribas
Departament de Ciències
2006 May 02
1
Is there a bug in CrossTable (gmodels)
Library gmodels include a function CrossTable that is useful for
crosstabulation. In the help, it is indicated that one can call this
function as CrossTable(data), were data is a matrix. However, when I try
to use this option, it doesn't help. Any idea? Is there a bug?
Thanks for your help.
Prof. Albert Sorribas
Grup de Biomatem?tica i Bioestad?stica
Departament de Ci?ncies M?diques B?siques
2006 Apr 26
3
Were to find appropriate functions for a given task in R
This is a generic request concerning were to look for finding
appropriate information on a precise procedure in R.
I’m using R for teaching introductory statistics and my students are
learning how to deal with it. However, I find it difficult to locate
some of the procedures. For instance, for basic crosstabulation, it is
obvious that basic functions as table, ftable, and prop.table can be
used.
2008 Aug 14
1
AMI and extensions.conf
Hello
I'm looking for a wayy to modify extensions.conf
It seems that PutConfig AMI command is not supposed to work on extensionsq.conf
Any ideas?
Thanks
Vadim
1997 Nov 29
1
Workgroup popups with samba
Hello,
I'd like to know if there's a way to receive Workgroup's winpopup messages
with samba under linux.
If it's impossible, what are the technical problems caused by such a request ?
Thank you
--
,,,
o +
+-----------------------ooO---(_)---Ooo---------------------------+
| "Xylo", a.k.a
2007 Aug 02
1
how to plot a differential equation?
Hi,
I would like to plot the following equation:
dF(x)/dx=(k1+k2F(x))(1-F(x))
where k1 and k2 are parameters that I have estimated already.
How can I plot the curve in R?
Thanks!
Montserrat Rue
Universitat de Lleida (Spain)
[[alternative HTML version deleted]]
2016 Apr 04
4
Test for Homoscedesticity in R Without BP Test
Respected Sir,
I am doing a project on multiple linear model fitting and in that project I
have to test Homoscedesticity of errors I have google for the same and
found bptest for the same but in R version 3.2.4 bp test is not available.
So please suggest me a test on homoscedesticity ASAP as we have to submit
our report on 7-04-2016.
P.S. : I have plotted residuals against fitted values and it is
2008 Apr 09
2
How to estimate a hazard ratio using an external hazard function
Hi,
I would like to compare the hazard functions of two samples using the
Cox proportional hazards model. For sample 1 I have individual time-to-
event data. For sample 2 I don't have individual data, but grouped
data that allows to obtain a hazard function.
I am wondering if there is an R function that allows to obtain a
hazard ratio of the two hazard funtions (under the
2006 Jan 26
3
Samba daemons hang trying to lock locking.tdb
Hi,
A day or so after starting samba, some daemons (diferent forks) begin to
hang. Then, the WinXP clients hang too completely.
When I try to figure out what is happen, I see that smbd daemons hangs
always in a fcntl64() call:
# strace -p 6414
Process 6414 attached - interrupt to quit
fcntl64(14, F_SETLKW64, {type=F_WRLCK, whence=SEEK_SET, start=3684, len=1} <unfinished ...>
Process 6414