similar to: How to add multiple ablines

I really want to do this: abline( a=tan(-kT*pi/180), b=kY-tan(-kT*pi/180)*kX ) where kX,kY and kT are vectors of equal length. But I can't do that with abline unless I use a loop, and I haven't figured out the least unelegant way of writing the loop yet. So is there a way to do this without a loop? Or if I am to resort to the loop, what's the best way of doing it considering that I

Hi folks, Win7 64bit R 1.12.0 I run following command on R:- > ToothGrowth > attach(ToothGrowth) > plot(dose,len) > matrics=lm(len~dose) > abline(metrics) Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : plot.new has not been called yet Only a grey diagram is displayed without content > plot(abline(metrics)) Error in int_abline(a = a, b = b, h = h, v

hi List and Manuel, I have encounter the following problem with the function "lineplot.CI".? I'm running R 2.10.1, sciplot 1.0-7 on Win XP.? It seems like it's a scoping issue, but I couldn't figure it out. Thanks! ...Tao > lineplot.CI(x.factor = dose, response = len, data = ToothGrowth)??? ## fine > lineplot.CI(x.factor = dose, response = len, data = ToothGrowth,

Dear all How does one resolve the variance between the actual availability of data sets in the default implementation with those mentioned in the documentation? I am unable to attach some of the datasets, even though help() is available for the same datasets. For example, "ToothGrowth" is a dataset that is widely mentioned in the documentation as being part of the base install.

Hi to all, Does anybody knows why this is giving an error? data(ToothGrowth) # Two-way design with options bargraph.CI(dose, len, group = supp, data = ToothGrowth, xlab = "Dose", ylab = "Growth", cex.lab = 1.5, x.leg = 1, col = "black", angle = 45, cex.names = 1.25, density = c(0,20), legend = TRUE) Error in dn.call[[1]] :

Thank you for your guys reply for my previous question. But I got one more question about the boxplot. With the code in the R-help: boxplot(len ~ dose, data = ToothGrowth, boxwex = 0.25, at = 1:3 - 0.2, subset = supp == "VC", col = "yellow", main = "Guinea Pigs' Tooth Growth", xlab = "Vitamin C dose mg", ylab =

Hi all, I I need to manipulate the x axis in a stripchart. I will use one of the data sets included in R to explain what I need to do. attach(ToothGrowth) stripchart(len[supp=='VC']~dose[supp=='VC'], vertical=TRUE, group.names=c('A','A','A')) stripchart(len[supp=='OJ']~dose[supp=='OJ'], add=TRUE, vertical=TRUE, at=c(1:3)+.1,

Hi, I followed the following sequence of commands after installing 'XGobi': >library(xgobi); >data(ToothGrowth); >xgobi(ToothGrowth); I got the following warning message: C:/R13 1.0/rw1030/library/xgobi/scripts/xgobi.bat -vtitle 'ToothGrowth' -std mmx C:/DOCUME~1/h1452992/LOCALS 1/Temp/xgobi-TthGr5705 Warning message: argument `pkg' is deprecated. Use

Dear List, ? I am creating a boxplot with two subsets, very similar to the example by Roger Bivand at ?boxplot (reproduced below).? I am trying to change the labels on the x-axis to have one number to cover both subsets.? I can do this in other plots by using axis=FALSE followed by a separate axis() command.? I have also tried variations in the names= argument but can't get it to work.?

I'm trying to plot a scatterplot of two variables using pch to plot different characters based on a third factor. Here is my example > data("ToothGrowth") > attach(ToothGrowth) > levels(supp) [1] "OJ" "VC" > plot(len ~ dose,pch=as.numeric(supp)) > legend(locator(1),pch=as.numeric(supp),legend=levels(supp)) The command as.numeric(supp) returns 2 2

Dear list, How can I add the mean and standard deviation to each of the boxplots using the example provided in the boxplot function? boxplot(len ~ dose, data = ToothGrowth, boxwex = 0.25, at = 1:3 - 0.2, subset = supp == "VC", col = "yellow", main = "Guinea Pigs' Tooth Growth", xlab = "Vitamin C dose mg",

Dear all I have written a function to perform a very simple but useful task which I do regularly. It is designed to show how many values are missing from each variable in a data.frame. In its current form it works but is slow because I have used several loops to achieve this simple task. Can anyone see a more efficient way to get the same results? Or is there existing function which does

Hi all, I'm attempting to make a quite-specific plot where the axes cross at the origin and with gridlines for guidance. I've been using ablines to create the reference lines because I want a lot of control as to where they are placed on the axis. This command works very well for such control. However... These ablines don't seem to work when I specify the origin as 0,0. They go

Hello, When using bargraph.CI in package sciplot can the bars for each group be different colors? How do I select the color for each group? When I use this instead of the default (SD vs SE): bargraph.CI(x.factor = dose, response = len, data = ToothGrowth, ci.fun= function(x) c(mean(x)-sd(x), mean(x) + sd(x)) ) Am I getting 95% CI bars? Thank you kindly, Michael Just

Hi all, I have a csv matrix "KT.csv" and it has Inf and NA I want to calculate the mean of each row so I use rowMeans(KT,na.rm = TRUE) but with this Inf cannot be omminted. I?m trying to use before running rowMeans(KT,na.rm = TRUE) KT<-range(KT,finite=TRUE) but it doesn?t works... Do you know a simple way to ommit Inf en the calculations? Many thanks. I have tried also to

https://bugzilla.mindrot.org/show_bug.cgi?id=2680 Bug ID: 2680 Summary: Regression in server-sig-algs offer in 7.4p1 (Deprecation of SHA1 is not being enforced) Product: Portable OpenSSH Version: 7.4p1 Hardware: Other OS: Linux Status: NEW Severity: enhancement Priority: P5

Hi, I'm trying to fit nonlinear mixed effects model using nlme function but getting an error message. Here is what I have: fitted_model = nlme(scores~spline(b1,b2,b3,kt,time), fixed = list(b1~1, b2~1, b3~1, kt~1), random = b1+b2+b3~1, groups= ~id, data = sdat, start = c(b1=3.5,b2=2,b3=.60,kt=3.5),verbose=T) Error: Error in

I'm running an LSODA to generate some graphs, but I need to stop at a certain point and use those values to generate another LSODA output. This is working fine, but when I try to run the second LSODA, I get the "Error in eval(expr, envir, enclos) : object 'N' not found". Any ideas what can be causing this? I have no object 'N' anywhere in the script. I made an

Hi, Thanks for dovecot. It's the best IMAP server ever. Okay, here's my problem: I use mbox INBOXes from /var/spool/mail/%u, and the INDEX lies in /var/spool/mail/dovecot/%d/%u. Users' other folders are found in /%h/mail. Now, I notice that when a user is accessing her INBOX, the imap process cwd's to the user's home, which I'd like to prevent as much as possible.

ComRades, How do you put a superscript on a scalable delimiter? I want to put 'b' as the power of the expression in the following plot: t <- 1:25 K <- 0.2 y <- ((1-exp(-K*t))/(1-exp(-K*t)*exp(K)))^3 plot(t,y,"l",main="K=0.2, b=3") text(15,5,expression(bgroup("(",frac(1-e^-Kt,1-e^-Kt*e^K),")"))) Plotmath examples in demo(plotmath) do not