similar to: linear mixed model with nested factors

Does anyone else miss email-friendly character graphics such as the following example, produced using Minitab? Histogram of C6 N = 478 N* = 21 Each * represents 2 observation(s) Midpoint Count -12 16 ******** -11 53 *************************** -10 63 ******************************** -9 83

Dear R-Friends, Sorry to bother everyone with my earlier question: "Do we have similar R function to work like Minitab's INDICATOR?" I found the way to make things work: e.g., > x <- c(2,2,5,3,6,5,NA) > model.matrix(~ factor(x) - 1) factor(x)2 factor(x)3 factor(x)5 factor(x)6 1 1 0 0 0 2 1 0 0 0

hi folks, Bit of a newbie, but I've spent a fair bit of time looking for an answer on this, with no joy. Can anyone help me? Dataset: A single column of values in a csv file (eg. 52, 53, 54, 85, etc) Goal: In Minitab, you have what they call a dot plot. It's a histogram, where a single dot represents a set of identical values (eg. 57, 57, 57 would be one dot). Multiple dots are

Dear R experts, Can R make an usual dotplot just like Minitab and other softwares? I have the following data, and can use dotchart to graph a dotplot: y=c(2.873438152e-01, -8.732895642e-01, 4.579001889e-01, 1.047395204e+00, 8.491182299e-02 , -1.938007105e+00, -1.273708343e+00, 9.848010588e-05, 7.238490734e-01, -1.490552717e+00) dotchart(y, xlab="10 observations from

Hi, list, suppose i have such a data frame: trash <- data.frame(cbind(seq(1:5),c('a','a','b','a','b'),c('b','a','b','b','a'))) names(trash) <- c('age','typeI','typeII') and I want to change all 'a's to be 0 and 'b's to be 1. temp <- as.matrix(trash)

Dear R-help, Using R 1.2.1 on Windows98 to run a factor analysis on a 64x150 matrix of data generated from a simulation model, factanal() reported that it failed to find a solution. Looking at the factanal code, I see the immediate condition that triggered the result: if (best == Inf) stop("Unable to optimize from these starting value(s)") So I am sure factanal() is giving

Hi, I have found a minor problem with leaps(). In 1.3.1 under Windows 2000 I seem to only be able to obtain values for one statistic at a time. That is choosing method=c("Cp","adjr2","r2") just gives Cp values. To mimic the output of Minitab's MTB > BReg 'Fertility' 'Agriculture'-'Infant.Mortality' ; SUBC> NVars 1 5; SUBC>

Hello all I'm trying to learn how to fit a nested model in R. I found a toy example on internet where a dataset that have?3 areas and 4 sites within these areas. When I use Minitab to fit a nested model to this data, this is the ANOVA table that I got: Nested ANOVA: y versus areas, sites Analysis of Variance for y Source DF SS MS F P areas 2 4.5000 2.2500

Dear R-Friends, Do we have a similar R function to work like Minitab's INDICATOR? INDICATOR C2, C11-C15 C2 C11 C12 C13 C14 C15 2 1 0 0 0 0 2 1 0 0 0 0 5 ------> 0 0 0 1 0 3 0 1 0 0 0 6 0 0 0 0 1 5 0 0 0 1 0 * * * * *

I have written a function to emulate minitab's QQ plotting output (with SW test and AD test results on the graph): mtab.norm<-function(x) { library(nortest) library(lattice) x<-as.numeric(x) x<-as.vector(x) plot.ht<-4.6 plot.wd<-4.6 pt.ht=plot.ht/5 txt.sz<-(plot.ht/7.5) X11(width=plot.wd, height=plot.ht, bg='gray96') qqplot(x, pch=16, cex=pt.ht,

I'm having trouble using the read.mtp function in the foreign package to import datasets from minitab (.mtw) format. Specifically, each file I try to import fails to load any data beyond a header row stating the version of Minitab that saved the dataset. I get this error: incomplete final line found by readtableHeader on 'income.mtw' The dataset appears to be complete (ie, using

Hi I am performing an analysis of variance with two factors, each with two levels. I have differing numbers of observations in each of the four combinations, but all four combinations *are* present (2 of the factor combinations have 3 observations, 1 has 4 and 1 has 5) I have used both anova(aov(...)) and anova(lm(...)) in R and it gave the same result - as expected. I then plugged this into

Dear all, I'm trying to get a Brownian bridge kernel (kernelbb) for each combination of two consecutive animal locations (see commands below) and put them, with a loop, inside a list. It works well at the beginning but after 42 runs, it appears the following warning : >Error in seq.default(yli[1], yli[2], by = diff(xg[1:2])) :  >  invalid (to - from)/by in seq(.) I looked at the

Dear R list, I have been trying to do a linear model, extracting the effect of a covariate.... and the results do not match, when I do it with other programs (e.g. minitab).... so it is obvious that I was doing something wrong. Whan I do it with minitab, I have this results: (sector is a factor and depth is the covariate): Source DF Seq SS Adj SS Adj MS F P

Dear All, Would like to ask the inconsistency in the autocorrelation from R with SPSS/Minitab. I have tried a dataset x with 20 data (1-20) and ask R to give the autocorrelation of different lags using the command < acf(x, lag.max=100, type = "correlation"), However while SPSS and Minitab give the same answers (0.85 for lag1), R gives 0.3688 which is much smaller. Obviously, the

Hello, Given a probability function: p(x) = 12 / (25(x+1)) , x=0, 1, 2, 3 we generate the following values: C1 C2 0 0.48 1 0.24 2 0.16 3 0.12 Now, I'm supposed to create 50 random values using this table. In MiniTab, I simply selected Calc -> Random Data -> Discrete, and selected the columns, and it created 50 random values in a new column.[1] How do I do the

Hi, I ran a small data set from a factorial experiment through R, Minitab and JMP... the result from R is significantly different from what Minitab or JMP give... The data set is at the following link: http://www.personal.psu.edu/nug107/Uploads/2x3_16repsANOVA.txt The first 5 columns are the factors and the next three are responses. In particular, for the response beta11MSE, two of the

I recently downloaded R for Windows, running on Win XP. I'm trying to create a perspective plot but not having any luck after reading the R manual and several examples found on the Internet. I have a 100 x 100 matrix of Z data as a tab-delimited text file exported from Minitab. I read this in to R using read.delim; this seemed to go ok. I created X and Y using seq() to get 100 divisions

Dear all, Does R or S-plus or any of their packages provide any command to form any of the following procedures to find Best Regression Equation - 1. 'All Possible Regressions Procedures' (is there any automated command to perform 2^p regressions and ordering according to criteria R2(adj), mallows Cp, s2- by not setting all the regression models manually), 2. 'Backward

Hello - I have two files which play back in VLC correctly. They were encoded using this command: ffmpeg -y -i '/tmp/2e905b981d30324f9a0664cbabb3f59a20170406-21128-r0neh9.wav' -acodec libvorbis '/tmp/2e905b981d30324f9a0664cbabb3f59a20170406-21128-r0neh920170406-21128-146yyex.ogg However ov_pcm_total reports a zero length for one of the two files. Is this a problem with my use of the