similar to: lme in loop help

Sorry I had typo in previous email, this typo corrected version: Dear R experts I have simple question, please execuse me: #example data, the real data consists of 20000 pairs of variables K1 <- c(1,2,1, 1, 1,1); K2 <- c(1, 1,2,2, 1,2); K3 <- c(3, 1, 3, 3, 1, 3) M1a <- rep( K1, 100); M1b <- rep(K2, 100) M2a <- rep(K1, 100); M2b <- rep(K1, 100) M3a <- rep(K1, 100); M3b

Hi there, I'm using R1.0.1 Windows 98. This file contains some inputs and an aov function code. Can someone check it for me? Somehow I got completely different answer when typing them in R and in Splus. Splus gives me this: > summary( Turnip.aov ) Error: Blocks Df Sum of Sq Mean Sq F Value Pr(F) Residuals 3 163.7367 54.57891 Error: Plots %in% Blocks

I would like to fit regression models of the form y ~ ns(x, df=splineDF) where splineDF is passed as an argument to a wrapper function. This works fine if the regression function is lm(). But with lme(), I get two different errors, depending on how I handle splineDF inside the wrapper function. A workaround is to turn the lme() command, along with the appropriate value of splineDF, into a text

I would like to fit regression models of the form y ~ ns(x, df=splineDF) where splineDF is passed as an argument to a wrapper function. This works fine if the regression function is lm(). But with lme(), I get two different errors, depending on how I handle splineDF inside the wrapper function. A workaround is to turn the lme() command, along with the appropriate value of splineDF, into a text

A maybe trivial and stupid question: In the case of a lm or glm fit, it is quite informative (to me) to have a look to the null deviance and the residual deviance of a model. This is generally provided in the print method or the summary, eg: Null Deviance: 658.8 Residual Deviance: 507.3 and (a bit simpled minded) I like to think that the proportion of deviance 'explained' by the

Hi everybody, I have an experiment examining risky choice behavior where two groups of subjects were unevenly divided across two different MRI scanners while they performed a task. Each subject's data was recorded once and only once on a particular scanner. The table describing the distribution of subjects across the scanner (3TE and 3TW) and groups is below. 3TE 3TW Group1 10

Hi everybody, I've written a function to run an LME model on data derived from functional magnetic resonance images. When I run the function with contrasts included I get the following error Error in eval(expr, envir, enclos) : object 'inModelFormula' not found I think it has something do do with the way contrast evaluates arguments, but I've got no idea how to fix it. The code

Hello, Dr. Viechtbauer. I am trying to perform a meta-analyis on a group of before-after studies using Metafor. I read your webpage including your correspondence with Dr. Dewey (https://stat.ethz.ch/pipermail/r-help/2012-April/308946.html), who also conducted a similar study. These information is very hepful, but I have one additonal question which I wonder if you can give me some instruction.

Dear all My dataframe has > 80,000 variables which I can not everytime load into R using *.txt files (read.table option), cost me time and sometime computer decomes not responsive. So I need a way to save my dataframe in my workdirectory as such. Execuse me if the problem is too simple. I tried the following, I do not know what is wrong with it: Example data: x <- 1:10 y <- 21:30 z

Dear list I need your help: Execuse me for my limited R knowledge. #example data set set.seed (134) lm=c(1:4) block = c(rep(lm,6)) gen <- c(rep(1, 4), rep(2, 4), rep(3, 4), rep(4, 4),rep(5, 4),rep(6, 4)) X1 = c( rnorm (4, 10, 4), rnorm (4, 12, 6), rnorm (4, 10, 7),rnorm (4, 5, 2), rnorm (4, 8, 4), rnorm (4,7, 2)) X2 = X1 + rnorm(length(X1), 0,3) yvar <- c(X1, X2) X <- c(rep( 1,

[ I'm teaching log-linear models at present, and looking at Agresti's book. ] Often when producing analyis of variance or analysis of deviance tables it is important to get effects fitted in a specific order. At present in S / R this is done as follows: anova( glm( terms(formula, keep.order=T), ...) ) Can anyone think of a reason why keep.order should not become an argument to

Dear all, I have a dataframe with multiple observations and the levels as the last column, as in: d <- data.frame(A=sample(1:100,12),B=sample(1:100,12),levels=c(rep('A',4),rep('B',4),rep('C',4))) > d A B levels 1 77 40 A 2 14 18 A 3 56 7 A 4 46 27 A 5 63 35 B 6 80 21 B 7 3 54 B 8 93 76 B 9 5 46 C 10 16 53

Dear R-experts Execuse me for an easy question, but I need help, sorry for that. >From days I have been working with a large dataset, where operations are needed within a component of dataset. Here is my question: I have big dataset where x1:.....x1000 or so. What I need to do is to work on 4 consequite variables to calculate a statistics and output. So far so good. There are more vector

Hi all, I am having memory allocation problem with my R 2.2.1 for Mac OS. The following is the error message that I get. I do not get this message if I break down the large dataset in to sub datasets. I think breaking up the dataset is not a sustainable solution in the long run. The data that I am analysing is essentially big, and it would be reasonable to do the analyis on the whole dataset

My objective is to look at differences in two species of fish from morphometric measurements. My morphometric measurements are head length, eye diameter, snout length, and measurements from tail to each fin. I want to use discrimanant function analyis to determine if there are differences between the two species. I am familiar with R but new to discrimannt function analysis. I want to learn

Hi: Execuse me if this has been mentioned before, I''ve looked around the archive but did not find anything about the problem I met. I''m trying to use tcindex to classify the the packet outgoing to one interface. I do as the following : tc q a dev eth0 root handle 10:0 cbq bandwidth 10Mbit avpkt 1500 tc c a dev eth0 parent 10:0 classid 10:1 cbq bandwidth 10Mbit rate

Hello, I have installed Ubuntu 9.10 and it has ext4 filesystem. I have one very elementary question. When I say my filesystem is "ext4", which directories are part of it. I mean from the root I can see some directories such as /proc, /tmp, /dev etc. Are they store on the disk which have formatted with ext4, of certain files resides somewhere else. I am only sure only about /home

Hi! I''m from Russia. It''s my first post in this forum. Execuse for my bad language, but i hope you''re help me! I try to deploy rails app on apache + fastcgi hosting. .htaccess file contents: AddHandler fcgid-script .fcgi AddHandler cgi-script .cgi Options +FollowSymLinks +ExecCGI RewriteEngine On RewriteRule ^$ index.html [QSA] RewriteRule ^([^.]+)$ $1.html [QSA]

Hi, extractRangeChecksFromBranch uses BranchProbabilityInfo to decide whether its worth trying the InductiveRangeCheckElimination transformation. For the following example: void split() { for (int i = 0; i < 100; ++i) { if (i < 99) do_something() else do_something_else() } } But the reported BPI is reported as 50/50 to whether do_something will be called, but we

Adding Maxim On Apr 9, 2018, at 10:06 AM, Sam Parker via llvm-dev <llvm-dev at lists.llvm.org<mailto:llvm-dev at lists.llvm.org>> wrote: Hi, extractRangeChecksFromBranch uses BranchProbabilityInfo to decide whether its worth trying the InductiveRangeCheckElimination transformation. For the following example: void split() { for (int i = 0; i < 100; ++i) { if (i < 99)