similar to: problem with abine(lm(...)) for plot(y~x, log='xy')

hello, I only want to get the slope of a linear regression of ca. 100 variables against time. I can do for each response (100 times) var1.lm <- lm(response~predictor) but I thought that there might be an easier way of doing this. If I am including more variables it is doing a multiple regression and the output (slope) differs. any idea? thanks Martin

Hi, All: Does 'nlme' require the nonlinear function to be differentiable? I'm fitting structural change models to many related time series, and I'd like to improve the estimates of the change point times through some type of pooling. Unfortunately, I've so far been unable to get 'nlme' to work for me. The following toy example is the closest I've come

** Delivery incomplete ** There was a temporary problem delivering your message to curtiskwong101 at gmail.com. Gmail will retry for 46 more hours. You'll be notified if the delivery fails permanently. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.centos.org/pipermail/centos-docs/attachments/20180919/4030c740/attachment-0002.html>

** Delivery incomplete ** There was a temporary problem delivering your message to curtiskwong101 at gmail.com. Gmail will retry for 47 more hours. You'll be notified if the delivery fails permanently. The response was: Receive rate too high -------------- next part -------------- An HTML attachment was scrubbed... URL:

** Delivery incomplete ** There was a temporary problem delivering your message to curtiskwong101 at gmail.com. Gmail will retry for 47 more hours. You'll be notified if the delivery fails permanently. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.centos.org/pipermail/centos-docs/attachments/20180919/7249d4a2/attachment-0002.html>

Hi, I'm trying to plot a map with the location of the stations that I am using require(MASS)require(maps) mapaBrasil <- read.table("BASE_BR.txt",sep="") latlonRS <- read.csv2("coordenadas.txt",h=F,sep="") png(filename="mapa_grupos.png",width=600,height=800) # criating a map for RSmap(xlim=c(-80,-35),ylim=c(-45,10),

I am currently trying to use R to construct a regression model to explain output based on temperature. I have combined my output and temp data into a notepad file. there is no problem with loading the data into R. > data.df output temp 1 850 17 2 849 17 3 905 17 4 925 17 5 1043 19 6 1104 20 7 1097 18 8 979 19 9 926 18 10 1133 18 ~~ ~~ 240 1124

Is there a way to teach xy.coords, when given log="xy", or just "x" or "y" separately, to do a decimal log10 instead of the natural log? Cheers, Alexy

I'm now trying to read.zoo in a rather long txt file with two columns: date/time and value in kW e.g. 432.2189 The read.zoo function finally ran without errors but not sure it is correct. I took the header off, and put in commas and added a " at the beginning and " at the end. z=read.zoo("Kevin-0-comma-ITPower.txt","%m/%d/%Y

Dear R-help, I am looking for a function that will plot error bars in x- or y-direction (or both), the same as the Gnuplot function 'plot' can achieve with: plot "file.dat" with xyerrorbars,... Rsite-searching led me to the functions 'errbar' and 'plotCI' in the Hmisc, gregmisc, and plotrix packages. As I understand the descriptions and examples, none of

Hi - I get the following complaint from R windows 632 on this command > plot.default(ppmZnSGS+1, ppmZnALS+1, main = "5 Lab Check ppmZn SGS vs ALS",log="xy") Error: missing value where logical needed The following doesn't complain > plot.default(ppmZnSGS+1, ppmZnALS+1, main = "5 Lab Check ppmZn SGS vs ALS",log="") log="xy" is included

In ?xy.coords it says: If 'y' is missing and 'x' is a formula: of the form 'yvar ~ xvar'. 'xvar' and 'yvar' are used as x and y variables. list: containing components 'x' and 'y', these are used to define plotting coordinates. time series: the x values are taken to be 'time(x)' and the y

Hi, With current R devel and R 2.13.1: > aa=1:5 > save(aa, file="aa.rda", compress="xy") Error in identical(compress, "gzip") || compress : invalid 'y' type in 'x || y' Thanks! H. > sessionInfo() R version 2.13.1 (2011-07-08) Platform: x86_64-unknown-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_CA.UTF-8 LC_NUMERIC=C [3]

Hi Trying to finish a simple xyplot, but the key colors don't match those specified on the plot. I have only been able to change the color of the key text, but not the points. Hours <- c (25.88,32.31,24.04,23.27,16.64,19.75,27.25,17.9,20.45,29.65,20.83,27.17, 28.42,28,18.85,17.4,19.75,25.34,23.23,29.85,28.96,25.6,14.67,20.55,25.77 ,

I would appreciate any help for the following question. How to increase the font size of labels for the x and y axes in the following code? --------------------------------- xyplot(lbxglu~lbxgh|eth,reg.dat.5, panel = p, col = 1:2, par.strip.text = list(cex = 1.25,col="brown"), xlab="A1C",ylab="Fasting Plasma Glucose", main="FPG") p <- function(..., col)

I have 2 vectors, x and y and have done an xy plot. I want to plot the label (name?) of the vector on the plot rather than the value. text(x,y, labels = x) gives me the value of x. text(x,y, labels = labels(x)) gives me something like c("text1","text2"..) plotted for each point text(x,y, labels = names(x)) gives nothing print(x) gives me [,1] text1

Hi. I'm trying to add some vertical lines to an xyplot in which x axis is a temporal variable (class=dates) and y axis a factor variable. This is my code: xyplot(factor(Abacus$Emisor)~Abacus$Dia,xlab="Date", ylab="Fish",pch=124, scales=list(x=list(format="%b %Y",tick.number=30)) , panel=function(...){ panel.xyplot(...) panel.abline(v=08/08/2008) }) I'd

Hi, Has anyone know about how to get the correct plot? I have use this R script (as below), so I expect the plot is based on x axis, but the result was opposite. Any suggestion will be great. library(IRanges) data <-read.table(file="~/q20snpref/illusmp454merbed",sep="\t",header=F) colnames(data)<-c("Scaffold","sca_position","coverage")

Hi, I am new to R. How can I get labels/different colours on the lines in my xy-plots. I plot "het" against "temp" conditional on "year" and "station". library(lattice) MyLines <- function(xi, yi, ...){ #Draws line in the panels while avoiding spaghetti-plots I <- order(xi) panel.lines(xi[I], yi[I], col = 1)} xyplot(het ~ temp | fstation,

Hi! Where is ozone.xy? I've found ozone in package fields, but cannot find ozone.xy Actually, is there any way to find in which package lives a particular dataset (if the package is not installed in the local system?) Thanks Dr. Agustin Lobo Instituto de Ciencias de la Tierra (CSIC) Lluis Sole Sabaris s/n 08028 Barcelona SPAIN tel 34 93409 5410 fax 34 93411 0012 alobo at ija.csic.es