Displaying 20 results from an estimated 500 matches similar to: "repeated measures ANOVA using a cov matrix"
2009 Feb 23
1
why results from regression tree (rpart) are totally inconsistent with ordinary regression
Hi,
In my analysis of impacts of insecticide-treated bednets on malaria, I
look at the relationship between malaria incidence and mosquito
behaviors. The condensed data set is copied here. Ordinary regression
(lm) shows that Incidence was negatively related to Mortality. This
makes sense because the latter reflected the strength of killing
mosquitoes by insecticide-treated nets. Since the
2008 Feb 19
4
[LLVMdev] 2008-01-25-ByValReadNone.c Failure
Hi all,
I'm seeing this failure on my PPC G4 box running TOT with llvm-gcc
4.2. Is anyone else seeing this? I'm sure it's related to the byval
stuff that's recently gone into LLVM. I'm attaching the output of
this command:
$ llvm-gcc -emit-llvm -O3 -S -o - -emit-llvm /Users/wendling/llvm/
llvm.src/test/CFrontend/2008-01-25-ByValReadNone.c
As you can see in it, there
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
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
2012 Aug 03
1
Multiple Comparisons-Kruskal-Wallis-Test: kruskal{agricolae} and kruskalmc{pgirmess} don't yield the same results although they should do (?)
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
2018 May 16
0
Systemfit
Sadly you failed to set your email program to send plain text and the data is corrupted at my end.
I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis.
I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need
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
2008 Mar 25
1
Subset of matrix
Dear R users
I have a big matrix like
6021 1188 790 290 1174 1015 1990 6613 6288
100714
6021 1 0.658 0.688 0.474 0.262 0.163 0.137 0.32
0.252 0.206
1188 0.658 1 0.917 0.245 0.331 0.122 0.148 0.194
0.168 0.171
790 0.688 0.917 1 0.243 0.31 0.122 0.15 0.19
0.171 0.174
290 0.474
2000 Dec 06
3
write.table
Good morning,
suppose the following:
m <- round(matrix(rnorm(16), ncol=4), 3)
a <- rev(c(0.01, 0.025, 0.05, 0.1))
rownames(m) <- a
colnames(m) <- c("0.25,0.75", "0.4,0.6", "0.1,0.9", "0.4,0.9")
m
0.25,0.75 0.4,0.6 0.1,0.9 0.4,0.9
0.1 1.034 -0.119 -1.213 0.619
0.05 0.035 1.074 0.525 1.671
0.025 -1.687 0.960
2010 Feb 17
2
extract the data that match
Hi r-users,
I would like to extract the data that match. Attached is my data:
I'm interested in matchind the value in column 'intg' with value in column 'rand_no'
> cbind(z=z,intg=dd,rand_no = rr)
z intg rand_no
[1,] 0.00 0.000 0.001
[2,] 0.01 0.000 0.002
[3,] 0.02 0.000 0.002
[4,] 0.03 0.000 0.003
[5,] 0.04 0.000 0.003
[6,]
2007 Aug 10
7
Help wit matrices
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
2007 Aug 09
2
Systematically biased count data regression model
Dear all,
I am attempting to explain patterns of arthropod family richness
(count data) using a regression model. It seems to be able to do a
pretty good job as an explanatory model (i.e. demonstrating
relationships between dependent and independent variables), but it has
systematic problems as a predictive model: It is biased high at low
observed values of family richness and biased low at
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
2011 Nov 01
3
factor level issue after subsetting
Dear list,
I cannot figure out why, after sub-setting my data, that particular item
which I don't want to plot is still in the newly created subset (please
see example below). R somehow remembers what was in the original data
set. A work around is exporting and importing the new subset. Then it's
all fine; but I don't like this idea and was wondering what am I missing
here?
Thanks!
2012 Apr 20
1
predictOMatic for regression. Please try and advise me
I'm pasting below a working R file featuring a function I'd like to polish up.
I'm teaching regression this semester and every time I come to
something that is very difficult to explain in class, I try to
simplify it by writing an R function (eventually into my package
"rockchalk"). Students have a difficult time with predict and newdata
objects, so right now I'm
2008 Aug 25
3
lmer4 and variable selection
Dear list,
I am currently working with a rather large data set on body temperature
regulation in wintering birds. My original model contains quite a few
dependent variables, but I do not (of course) wish to keep them all in my
final model. I've fitted the following model to the data:
>
2012 Aug 21
1
make check fails two tests on RHEL 6 build
I am installing R 2.15.1 onto RHEL 6, using gcc 4.7.0 with Intel MKL
10.3.7 and the following environment:
export BLAS_LIBS="-Wl,--start-group
/usr/caen/intel-12.1/mkl/lib/intel64/libmkl_gf_lp64.a
/usr/caen/intel-12.1/mkl/lib/intel64/libmkl_sequential.a
/usr/caen/intel-12.1/mkl/lib/intel64/libmkl_core.a -Wl,--end-group
-lpthread"
export LAPACK_LIBS="-Wl,--start-group
2009 Apr 17
4
Separating variables in read.table
If I have a table (we'll call it, "test") containing two columns (as below):
i x1 x2 x3 x4 x5 y
0 1.125 0.232 7.160 0.0859 8.905 1.5563
7 0.920 0.268 8.804 0.0865 7.388 0.8976
15 0.835 0.271 8.108 0.0852 5.348 0.7482
22 1.000 0.237 6.370 0.0838 8.056 0.7160
29 1.150 0.192 6.441 0.0821 6.960 0.3130
37 0.990 0.202 5.154 0.0792 5.690 0.3617
44 0.840 0.184 5.896 0.0812 6.932 0.1139
58
2003 Apr 24
1
write.table problem
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
2007 Aug 13
1
GML with tweedie: AIC=NA
Dear Catarina,
I prefer to leave the AIC value as NA for the tweedie GLM family
because it takes extra time to compute and is only occasionally
wanted. It's easy to compute the AIC yourself using the dtweedie()
function of the tweedie package.
Best wishes
Gordon
At 03:05 AM 14/08/2007, Catarina Miranda wrote:
>Dear Gordon;
>
>I have also sent this email to R help mailing list,