similar to: sem error message

Dear R community, I need to plot the results of some simulations I did using QTL Cartographer. I am plotting LOD scores over three chromosomes. The three plot have to be one next to the other. The procedure I am using is: par(mfrow=c(1,3)) plot(x$x, x$y, ylim=c(0,35), type="l", col="blue", las=1, xaxs="i", yaxs="i", xlab="X Chromosome",

Hello, I'm working on some dialplan rules to pull multiple users into a conference call. I have some fairly straightforward rules which start up a new MeetMe conference, allow escape with the * key to invite more users, then use a features.conf sequence to bring the new user into the conference with ChannelRedirect. The problem I'm running into is the time in the MeetMe conference

I am wondering how to interpret the parameter estimates that lm() reports in this sort of situation: y = round(rnorm(n=24,mean=5,sd=2),2) A = gl(3,2,24,labels=c("one","two","three")) B = gl(4,6,24,labels=c("i","ii","iii","iv")) # Make both observations for A=1, B=4 missing y[19] = NA y[20] = NA data.frame(y,A,B) nonadd = lm(y ~

Hi all: Can someone advice me on how to hold the residuals Mean sq value on a string so it can be used in other calculations. I was trying something like this: Msquare<-dfr$Mean sq but fails..Thanks dfr <- read.table(textConnection("percentQ Efficiency 1.565 0.0125 1.94 0.0213 0.876 0.003736 1.027 0.006 1.536 0.0148 1.536 0.0162 2.607 0.02 1.456 0.0157 2.16 0.0103

I am having trouble using back ticks with the R extractor function 'predict' and an lm() model.? I'm trying too construct some nice vectors that can be used for plotting the two types of regression intervals.? I think it works with normal column heading names but it fails when I have "special" back-tick names.? Can anyone help with how I would reference these?? Short of

https://bugs.freedesktop.org/show_bug.cgi?id=85570 Bug ID: 85570 Summary: DPMS does not turn off LCD backlight on G73 (NV4B) [GeForce 7600 GS] Product: xorg Version: unspecified Hardware: Other OS: Linux (All) Status: NEW Severity: normal Priority: medium Component:

Dear forum members, How can I force a negative slope in a linear regression even though the slope might be positive? I will need it for the purpose of determining the trend due reasons other than biological because the biological (genetic) trend is not positive for these data. Thanks. Julia Example of the data: [1] 1.254 1.235 1.261 0.952 1.202 1.152 0.801 0.424 0.330 0.251 0.229

?s 17:38 de 30/11/2023, Robert Baer escreveu: > I am having trouble using back ticks with the R extractor function > 'predict' and an lm() model.? I'm trying too construct some nice vectors > that can be used for plotting the two types of regression intervals.? I > think it works with normal column heading names but it fails when I have > "special"

Hi, I've stumbled upon a problem which I can't solve easily with the available options in tinc - at least as far as I see. If enlightenment is all I need, I'll happily accept pointers ;) I try to establish a connection between two hosts. Each host has multiple addresses assigned to it's internet interface. A stripped down list would be: Host 1: 2001:780:0:1e::1

Dear all, I am unable to understand how Hline() works in tabular(). I've read the vignette and the help page, and here this example compiles perfectly fine: latex( tabular( Species + Hline() + 1 ~ Heading()*mean*All(iris), data=iris) ) However, if I try it on my own data it fails. Consider this: set.seed(1) Xa <- data.frame(p=rep(c("First group","Second

Dear All, I am interested in using R to summarize data so that I can THEN do some additional analyses using behavioral data. I study birds, mostly, and often behavioral data that includes timing. So, for example, each line may be the occurrence of some behavior, and it would include identifiers for species, individual, date, time of observation and the various behavioral observations that may

I have a stock price dataset a snippet of which is: > plcm60[1:15, c(1,3,4,5,6,7)] DATE BIDLO ASKHI PRC VOL RET 1 1/2/03 9.450 9.79 9.700 1531819 0.018907 2 1/3/03 9.670 9.94 9.940 1582192 0.024742 3 1/6/03 9.830 10.05 9.960 1843298 0.002012 4 1/7/03 9.835 10.38 10.350 1412441 0.039157 5 1/8/03 10.220 10.67 10.260 961400 -0.008696 6 1/9/03 10.280

Hi, I am trying to figure out how to lable the second y-axis after the following codes: plot(x, y, xlab="time", ylab="pay1" ) par(new=TRUE) plot(x,y2, ann=FALSE, xaxt="n", yaxt="n", pch=7 ) axis(4) Then, I want to label the second y axis "pay2". I tried "title(ylab="pay2")", but it put this lable on the first

Hi all, I attempted to add an EL5.1 client to our puppet server (EL5), and after signing the client cert, got the error "Certificates were not trusted: hostname not match with the server certificate" I found the mailing list discussion and the relevant page: http://www.reductivelabs.com/trac/puppet/wiki/RubySSL-2007-006 As far as I can tell, my puppermaster''s cert CN matches

Good afternoon, We are running callgrind and cg_annotate version 3.6.1 on Centos Linux Version 5.5 x86_32. One month ago Mr. Josef Weidenorfer issued a special patch that fixed callgrind on Centos Linux Version 5.5 x86_32. We can now profile complex C++ programs which use our own shared library libmdMatchup.so. However, when we use version 3.6.1 cg_annotate callgrind.out.22533 --auto = yes

Hi everyone, I'm trying to print a table justified to the left, but it doesn't work. Any hints? KennArt <- data.frame(NR=c(171,172,174,175,176,177,181,411,980), TYP=c("K?rnermais", ?"Corn Cob Mix", "Zuckermais", "Mischanbau (Silo)Mais/Sonnenblumen", ?"Mais mit Bejagungsschneise in gutem landwirtschaftlichen und ?kologischen Zustand",

Hello again R help. This is a simple question perhaps(laughing as I type) with a simple answer. Is the multcomp package appropriate for using with a lme built under the nlme package. I know i can get it to work (i.e report p-values for a tukey test) but i am not unsure if this is appropriate for a linear mixed effects model. Below is the code i used lmm1 <- lme(asinh(density_recruit) ~

Hi all: I have always used SPSS to estimate weekly covariance based on a linear regression model but have to hard code the model Std. Error and the Mean-Square and then execute one week a the time. I was wondering if someone could give me an idea on how to estimate weekly(WK) covariance using the summary and anova of "dfr"(lineal model below). I have to do this for 52

Hi everyone, I'm sorry if this turns out to be more a statistical question than one specifically about R - but would greatly appreciate your advice anyway. I've been using a logistic regression model to look at the relationship between a binary outcome (say, the odds of picking n white balls from a bag containing m balls in total) and a variety of other binary parameters:

Dear useRs, I am trying to model a time series with a transfer function. I think it can be put into the ARMA framework, and estimated with the 'arima' function (and others have made similar comments on this list). I have tried to do that, but the results have so far been disappointing. Maybe I am trying to make 'arima' do something it can't... The data are time series of