similar to: Matrix function

I posted this earlier, but am not certain that it was in fact posted, so I will try again. Hi, I have an equation Xcp and would like to look at errors affecting it. I am applying errors of -.58 to .1 to R and -.45 to .47 to Xc. I would like to look at all combinations. I set up a matrix function as follows and it does not work. On the last line, when it run, it does not seem like the two X’s

Hi all and thanks in advance. I am regressing Time and Weight, and then predicting Weight at different Time. The format of the Time data is day/month/year. How can I get R to use time series data such as this? Keith -- M. Keith Cox, Ph.D. Alaska NOAA Fisheries, National Marine Fisheries Service Auke Bay Laboratories 17109 Pt. Lena Loop Rd. Juneau, AK 99801 Keith.Cox at noaa.gov marlinkcox at

Hi, I would like to add a color under a lm line and not the plotted line. Is this possible? In the example, I do not want the area under the curve red, but rather under the llm line. x=seq(0,5,len=100) y=-(x-5)^2 llm<-lm(y~x) plot(x,y) polygon(c(x,x[length(x)]), c(y, y[1]), col='red') abline(llm) keith -- M. Keith Cox, Ph.D. Alaska NOAA Fisheries, National Marine Fisheries Service

Hello, hopefully simple question, but I cannot find the answer. I need to change the color from the standard default. Still want the scaled colors, but need different colors for different graphs. Code is: wireframe(z ~ y*x, mat.df, drape = TRUE, zlab = list("Water mass error (%)",rot=92), zlim=c(-50,180), xlab = list("Resistance error

I would like to create a point plot with the following two sets of points: #1 plot(Day,Total) and #2 (Day,cons_hat). Total is the actual value seen and cons-hat is a predicted value. If possible, I do not want to stack them as they are quite long. (FYI, I did use the reshape on a previous post, thanks, but this one is different). Day Tank Tanks Total cons_hat 1 a a4 5.651017 5.59

I would like to calculate a first and second derivative and am having problems finding a simple solution. My syntax may be off as I am not a mathematician, so pardon ahead of time. data: t<-seq(0,4, by=.1) The function is: H(t) = t^3-6*t^2+5*t + 30 from here I plot the curve: plot(x,y ,xlab="x-values", ylab="f(x)", type="l") But would like to similarly plot the

I need a two sample t.test between M and F. The data are arranged in one column, x. Cant seem to figure how to run a two sample t.test. Not really sure what this output is giving me, but there should be no difference between M and F in the example, but summary p-value indicates this. How can I run a two sample t.test with data in one column. x=rep(c(1,2,3,4),2)

I need grayscale formatting for a wireframe. The only col.regions that I can find are color palettes are all colored: rainbow(n, s = 1, v = 1, start = 0, end = max(1,n - 1)/n, gamma = 1, alpha = 1) heat.colors(n, alpha = 1) terrain.colors(n, alpha = 1) topo.colors(n, alpha = 1) cm.colors(n, alpha = 1) The code follows: X11() library(lattice) par(family="serif", cex=1.2)

I have a data set with four columns and need to make boxplots from them. Data is as follows: tank Tanks Total cons_hat 1 a a4 5.651017 5.59 2 a a5 5.017499 5.29 3 a a6 4.894238 4.69 4 c c4 3.986347 3.40 5 c c5 4.099442 3.58 6 c c6 4.150522 3.64 7 h h4 5.187792 6.32 8 h h5 6.713422 6.44 9 h h6 5.168555

I need the parameters estimated for a non-linear equation, an example of the data is below. rm(list=ls()) Time<-c( 0, 0, 0, 0, 0, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 8, 8, 8, 8, 8) Level<-c( 100, 110, 90, 95, 87, 60, 65, 61, 55, 57, 40, 41, 50, 47, 44, 44, 42, 38, 40, 37, 37, 35, 40, 34, 32, 20, 22, 25, 27, 29) plot(Time,Level,pch=16) Keith

When using a wireframe, I need to move the colorkey from the "right" position (default0 towards the plot. I have also needed to adjust the height and used the code colorkey=list(T,space='right',height=.5) I have looked at documents (within levelplot) but cannot find a way to move the colorkey other than right, left, bottom and top. I do not understand corner interacts with x,

This data are kilojoules of energy that are consumed in starving fish over a time period (Days). The KJ reach a lower asymptote and level off and I would like to use a non-linear plot to show this leveling off. The data are noisy and the sample sizes not the largest. I have tried selfstarting weibull curves and tried the following, both end with errors. Days<-c(12, 12, 12, 12, 22, 22, 22,

OK, I need help plotting. I have column headings of Day, Wgt, Foodin, Rep, Grp and Tanks. Rep=c(1,2,3) and Tanks=c(a1,a2,a3,a4,a5,a6, c1,c2,c3,c4,c5,c6, h1,h2,h3,h4,h5,h6). I created a subset where I only would like Rep=2, and Tanks=c(a4,c4,h4) and would like to graph (points) of Wgt and Day. I would think that I only need 3 colors, but when I run with only 3, only 2 lines show up. When I add

Here is a runable program. When I plot Day and Wgt, it graphs all the data points. All I need is daily.sub1 plotted. I also need each "Tanks" to have its own col or pch. When I run it with the line with pch, it gives me nothing. rm(list=ls()) Trial<-rep(c(1,2),each=12) Tanks=rep(c("a3","a4","c4","h4"),each=3,2) Day=rep(c(1:12),2)

Hello, hopefully simple question, but I cannot find the answer. I need to change the color from the standard default. Still want the scaled colors, but need different colors for different graphs. Code is: wireframe(z ~ y*x, mat.df, drape = TRUE, zlab = list("Water mass error (%)",rot=92), zlim=c(-50,180), xlab = list("Resistance error

Greetings! I am trying to compare simulated environmental conditions from a model against a recruitment time series for a species of crab by first dropping 5 data points, and then using the remainder to attempt to simulate the missing data as a measure of best fit and using the following code: all.mat<-as.matrix(comb,ncol=ncol(comb),nrow=nrow(comb)) obs<-as.matrix(R2,24,1) mod<-all.mat

I have re-labeled tick marks on the x axis. The problem is that by using axes=FALSE, the axes disappears and when they are called back using axis(side=1)..etc. the axis on sides 1 and 2 do not meet at the bottom left corner of the graph. I would also like to have the 3rd and 4th axes in there as well, all meeting in their respective corners. par(mfrow=c(1,2)) gut<-c("Full",

I have long axis labels (8units) on the x axis and 2 unit number labels on the y axis. I need to move the x axis title without moving the y axis title. I have used mgp=c(title, label and line), but mgp moves both the x and y titles. The following puts the y axis title way way out to the left. Anyone know how to move just the x axis. X11() conif<-c("Control", "Switched

Dear R users, I am querying an Oracle database using sqlQuery() from the RODBC package. What I would like to do is embed my R Sql query in a for loop. For example suppose the data.frame "vessel" contains two columns: vessel[,1]= vessel id and vessel[,2]=permit year I am using "vessel" as an input in a SQL based function that require inputs of permit year and vessel id

Hello all, I have used describeBy to generate the following summary statistics. I simply need x and y error bars on a plot that has CQN (xaxis) and Price (yaxis). There should be four total points on the graph (one for each supplier). Using "error.crosses(desc$CQN, desc$Price)" does not work. group: a vars n mean sd median trimmed mad min max range skew CQN