Displaying 3 results from an estimated 3 matches for "_year".
Did you mean:
year
2005 Dec 15
2
Why is bubbles() creating empty png graphs?
...quot;_Y", year, ".png", sep="")
progress <- progress + 1
cat(progress, 'of', progress_end, " | ", fn, " | ")
flush.console()
if (!file.exists(fn))
{
q <- paste("select X,Y,RX from", tn, "where (_run=", run, ") AND
(_Year=", year, ")")
rs <- dbSendQuery(con, q)
positions <- fetch(rs, n = -1)
while (!dbHasCompleted(rs))
{
cat(".")
flush.console()
}
cat(length(positions$X), "points")
png(filename=fn)
bubble(positions, do.sqrt=FALSE, main=q)
#plot(positions)
dev.off()
pause()
ca...
2007 Feb 24
0
Multiple comparisons when interacction]
...bs=2;
input year plot block treat yield;
run;
proc mixed data=uno;
class treat year block;
model yield=block year treat treat*year;
lsmeans year treat /pdiff;
lsmeans treat*year /slice=year pdiff;
ods output diffs=dos;
run;
data tres;
set dos;
if year=_year;
proc print data=tres;
var year _year treat _treat estimate stderr df tvalue probt;
run;
Data are attached as a file: data.csv.
In fact, I do not know if this is the best approach to analyze the data:
- Should block be considered as random? We use the same file and
randomization every year...
2007 Jan 29
3
Multiple comparisons when interacction
In the model:
lm.1 <- lm(variable ~ BLOC + TIL * YEAR , data=selvanera)
I found TIL*YEAR interaction significant. Then I am trying to compare
means of the different levels of TIL inside every YEAR using:
mc.2 <- glht(lm.1, linfct = mcp(TIL*YEAR="Tukey"))
summary(mc.2, test = univariate())
but it does not work.
There is any way of doing this, like the SLICE option in